General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizationai source. It is being reieased in the interest of making avaiiabie as 
much information as possibie. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



I\I/\SA 

Technical Memorandum 83852 


A HIGH RESOLUTION 
SPECTROSCOPIC STUDY OF 
THE P2 band of hydrogen 
SULFIDE AND THE 1-0 BAND 
OF HYDROGEN IODIDE 

(MASA-T/l-d3d5^) k iilGii tlESOiUliCN 

SPECIiiOdCOPlC SlUti OF Th£ V SUB 2 bAhO of 

hlfOiiUJEil dOLFlDE AND IHt 1-0 EANL OF 

iiYDfiOGbd i.ODlD£ Pb.D. Ihesis - Maryidua 

Ucxv. (NASA) 155 p nC AOo/flF A01 CSCL 07D GJ/25 


Lawrence Larrabee Strow 


DECEMBER 1981 



National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt. Maryland 20771 


N 82- 1 
UUCldS 

1 b3b2 


A HIGH RESOLUTION SPECTROSCOPIC STUDY OF 
THE V2 BAND OF HYDROGEN SULFIDE 
AND THE 1-0 BAND OF HYDROGEN IODIDE 


by 

Lawrence Lambee Stiow 


Dissertation submitted to the Faculty of the Graduate School 
§, of the University of Maryland in partial fulfillment 

of the requirements for the degree of 

^ Doctor of Philosophy 

1981 


ABSTRACT 


Title of Dissertation; A High Resolution Spectroscopic Study of the »2 

Hydrogen Sulfide and the l~0 Band of Hydrogen Iodide 

Lawrence Larrabee Strow. Doctor of Philosophy, 1981 

Dissertation directed by: Dr. Thomas D. Wilkerson 

Research Professor 

Institute for Physical Science and Technology 

A tunable diode laser spectrometer has been constructed and used to study; (1) the 
effects of centrifugal distortion on the transition frequencies and strengths of the V2 
of H 2 S, and (2) nuclear quadrupole hyperfine structure in the 1-0 band of HI. Some 
additional transition frequencies in the band of H 2 S were measured with a high resolu- 
tion grating spectrometer. 

A total of 126 line frequencies and 94 line strengths in the vj band were measured. 
The average accuracy of the line frequency measurements was ±0.(XK)l6cm‘*. An error 
analysts of relative frequency measurements using a diode laser is given. The line strengths 
were measured to an average accuracy of about 3%. The eflect of the finite spectral 
width of the diode laser on the measurement of line strengths is discussed in detail. 

The observed H 2 S line frequencies were fit to Watson's AS and NS reduced Hamil- 
tonians in both the K and IIP coordinate representations in order to determine the best 
set of rotation-distortion constants for the upper state of the vi band. Successful fits 
were obtained for all the reduced Hamiltonians except the AS Hamiltonian in the HP 
representation. The characteristics of this Hamiltonian arc examined in detail to under 
stand its poor behavior relative to the other reduced Hamiltonians. Comparisons of 
observed line strengths in this band to calculated rigid rotor line strengths are also 
presented. 

Nuclear quadrupole hyperfine structur * in the l-O vibration-rotation band of HI. 
was observed with the diode laser .spectrometer. The R(0). P(l>. P(2). and P(3> lines of 


this band were recorded, all of which clearly exhibit nuclear quadrupole splittings. The 
upper vibrational state nuclear quadrupole coupling constant, determined from the ob- 
served splittings, was found to be ■ 1850MHz ± l2MHz, or 1.2% ± 0.7% larger than the 
ground state coupling constant. It is believed that this is the first observation of nuclear 
quadrupole structure in the vibration-rotation s.iec^rum of a diatomic molecule. 

The results of the H2S measurements described here should serve as a severe test of 
present and future theories of centrifugal distortion in asymmetric rotor spectra. This 
work also demonstrates the utility of diode lasers for highly accurate studies of infrared 


molecular spectra. 
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CHAPTER I 


INTRODUCTION 

Information about the geometry, potential function, and electric dipole moment func- 
tion of a molecule cm be obtained from a quantitative study of its vibration-rotation 
spectrum. As the experimental methods of infrared molecular spectroscopy improve, bet- 
ter molecular models can be constructed from the interpretation of the observed spectra. 
Accurate experimental studies of molecular vibration-rotation spectra are also important 
for a number of diverse applications such as laser isotope separation and the study of 
planetary atmospheres (including the Earth’s atmosphere). 

The development of tunable diode lasers has made it possible to obtain molecular 
spectra in the 500cm*' to 2500cm*' spectral region with essentially Doppler limited reso- 
lution. Very precise spectroscopic measurements can Ue made with diode lasers in a rela- 
tively straightforw.’rd manner without complex instrumentation. The high re!/)lution of 
diode lasers makes them very attractive for the measurement of absorption line strengths, 
an aspect of molecular vibration-rotation spectra that is much less studied than transition 
frequencies. Tunable diode lasers can also be used to measure transition frequencies if 
absorption frequency standards are available in the spectral region under study. 

In this work, a tunable diode laser spectrometer was constructed and used to study 
the effects of centrifugal distortion on the transition frequencies and strengths of the 
band of HiS. Some additional line frequencies in the V 2 band were measured with a high 
resolution grating spectrometer. A study of nuclear quadrupole hyperfine structure in the 
low J lines of the 1-0 band of HI was also made with the diode laser spectrometer and 
will be presented here. 

The main portion of this work is concerned with the r'T vibration-rotation band of 
H 2 S. centered at 1 IKOcm*'. While the spectra of light asymmetric rotors such as H 2 S and 
H 2 O have been understood qualitatively for many years, the quantitative characterization 


1 


of both their energy levels and transition strengths has remained a difficult problem. This 
is due in part to the relativ>‘ly large effects of centrifugal distortion on the spectra of light 
asymmetric rotors. Since the vibration of H 2 S is largely unperturbed by other vibrations of 
the molecule, this band is ideal for a study of the effects of centrifugal distortion without 
other complicating interactions. Also, the V| and vibrational bands of H 2 S exhibit what 
are perhaps the largest intensity perturbations due to centrifugal distortion ever observed 
in the spectrum of an asymmetric rotor. 

In a low resolution (»5cm"*) study of the 1^2 band, Emerson and Eggers^*) found 
that the P branch region of tb'c band is about five times less intense than the R branch 
region. They were able to qualitatively account for the observed intensity anomalies by 
including centrifugal distortion in a calculation of the line stiengths of this band. Since 
their work, no further studies of the strengths in this band have been made. In fact, no 
individual line strengths in the infrared spectrum of H 2 S have ever been measured. 

The only previous rotational analysis of the line frequencies in the ^^nd of H 2 S 
was performed by Allen and Plyler.^^) They assigned 55 transitions from 1080cm** to 
1257cm'* from a spectrum recorded with a resolution of 0.3 cm**. They were able to ob- 
tain values for the quadratic rotational constants A, B, and C in the upper vibrational 
state of the 1^2 band, but no centrifugal distortion constants were determined. 

In the present work, a total of 1 26 line frequencies and 94 line strengths in the 1^2 band 
were measured. The line frequencies were successfully fit to several of Watson’s reduced rota- 
tional Hamiltonians.*^) and a number of centrifugal distortion constants in the first excited 
state of the 1^2 vibration as well as the quadratic rotational constants were determined. 

Also, the measured line strengths were compared to calculated rigid rotor strengths in order 
to determine the magnitude of the centrifugal distortion perturbations. The line strengths 
measured here should serve as a stringent test of the existing theories of centrifugal distor- 
tion perturbations in light asymmetric rotor line strengths. 

In Chapter II. the derivation of Watson’s reduced rotational Hamiltonian is summarised 
as well as the basic theory of line strengths in asymmetric rotot vibration-rotation spectra. 
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Chapters 111 and IV deal with the experimental portion of this work. In Chapter 111, the 
diode laser spectrometer is discussed and several examples of spectra of H 2 S taken with 
this instrument are presented. The determination of the H 2 S line frequencies and 
strengths from the diode laser spectra is outlined in Chapter IV and a detailed error 
analysis given. 

' The least-squares analysis of the observed i '2 band frequencies using Watson's AS and 

NS reduced Hamiltonians is presented in Chapter V. Successful fits to the frequencies 
were obtained using the AS Hamiltonian in the I*” coordinate representation, and the NS 
Hamiltonian in both the F and IIF representations. The AS Hamiltonian in the Hi'' 
representation, though, was unexpectedly found to fit the frequencies very poorly. A 
large part of Chapter V deals with the characteristics of this particular Hamiltonian and 
why it fits the frequencies so poorly. 

In Chapter VI the results of the H 2 S line strength measurements are (X>mpared to cal- 
culated rigid rotor line strengths. The dependence of the observed line strength perturba- 
tions on J. Kg. and K^. is also examined. 

In Chapter VII, the study of nuclear quadrupole hyperfine structure in the 1-0 
vibration-rotation band of hydrogen iodide near 2200cm'i is presented. Using the diode 
laser spectrometer tlie R(0), P(l), P(2), and P(3) lines of this band were recorded, all of 
which exhibit nuclear quadrupole hyperfine structure. Nuclear quadrupole splittings, which 
arise from the interaction of the electric field gradient of the molecular electron cloud 
. with the nuclear electric quadrupole moment, are commonly observed in microwave pure 

rotational spectra. In infrared molecular spectra, these hyperfine splittings are usually 
masked by either the Doppler width of the absorption lines or by the low resolution of 
the spectrometer being used to record the spectrum. 

In the case of hydrogen iodide, its extremely large nuclear quadrupole moment and 
molecular weight combine to give relatively large hyperfine splittings with small Doppler 
widths. This allowed the resolution of some structure in the low J transitions of hydro- 
gen iodide using the diode laser spectrometer. From these spectra, the nuclear quadrupole 
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coupling constant of hydrogen iodide in both the ground and upper vibrational states was 
determined. It is believed that this is the first observation of nuclear quadrupole hype^ 
fine structure in the vibrational spectrum of a diatomic molecule. 


CHAPTER II 


THEORY OF ASYMMETRIC ROTOR VIBRATION-ROTATION SPECTRA 


II. A. Development of the Rotational Hamiltonian 

The complete molecular Hamiltonian for a polyatomic molecule can be written as 


P 2 p 2 

H = 2 ±L + ±L. + V__ + + V„.. (ll-l) 

"■* 2M 2M "" ” 


where P^^ and are the momentum operators for the nuclei and electrons. The poten- 
tial energy terms and come from the nucleus-nucleus, electron-electron, 

and nucleus-electron interactions respectively. A great simplification to H can be made 
using the Bom-Oppenheimer approximation, which is adopted in this woric. This approxi- 
mation assumes that since electrons are much lighter than nuclei, they move much more 
quickly. Therefore the electron’s motion at any time is the same as if the nuclei were fixed 
at their present instantaneous positions. This allows one to separate Schrodinger’s equation 
into an electronic part, and a vibration-rotation part, to be solved separately. In the elec- 
tronic wave equation the nuclear positions are taken to he adjustable parameters. In 
the vibration-rotation wave equation the electronic energy as a function of nuclear posi- 
tion becomes the potential energy term in the Hamiltonian. To solve the vibration- 
rotation problem we can treat the molecule as a set of point masses representing the 
nuclei lying in some potential field. With this model the classical Hamiltonian can be 
developed from which the quantum-mechanical Hamiltonian is obtained by proper replace- 
ment of coordinates and their conjugate momenta with the corresponding quantum- 
mechanical operators. 

For a molecule of N nuclei, there are 3N degrees of freedom. Three of these degrees 
of freedom describe the translational motion of the molecular center of mass thus leaving 
3N-3 degrees of freedom to describe the rotation and vibration of the molecule. The 
rotation can be described by introducing a coordinate system which rotates with the mole- 
cule, called the molecule-fixed axes, whose origin is at the molecule’s center of mass. Three 
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coordinates, such as the Euler angles, are required to give the orientation of this axis svstem 
with respect to space fixed axes. The remaining 3N>6 degrees of freedom are associated 
with the normal mode vibrations of the molecule, there being three of these for a three atom 
molecule such as hydrogen sulfide. The molecular vibrations are described relative to the 
rotating and translating molecule-fixed axes with the normal vibrational coordinates Q|. 

The ciassical expression for the molecular kinetic energy is obtained by transforming 
the kinetic energy in the space fixed frame, T = Vi 2 *£n* *he molecule-fixed axis 

frame. This was originally done in a rigoix>us fashion by Wilson and Howardl^^ who obtained 

T = - n„) (Jj - nj) + w E (II-:) 

where a. ^ = x, y, ?, in the molecule-fixed frame and i runs over the 3N-6 nonnal coordi- 
nates. In this expression = dT/dca^^ is the component of angular momentum, 
being the comjxinent of the angular velocity of the molecule-fixed axes, rijj is a quan- 
tity close to the vibrational angular momentum about the a axis but has some contributions 
from the rotational motion and is not conjugate to any part’cular angular coordinate. The 
term is a function of the nori;.al coordinates and is almost, but not quite, equal to the in- 
stantaneous inverse moment of inertia tensor. Pj = 3T/dOj is the momentum due to vibration 
and is conjugate to the normal coordinate Qj. The translational energy of the molecule has 
been left out of the kinetic enei-gy expression, ('omplete definitions of the terms used here 
can be found in several placesl-*'-^^. It is assumed that the potential V can be written solely 
as a function of the normal awrdinates Qj. In the limit of small vibrations only terms in 
Qj- will appear in the potential. 

The quantum-mechanical Hamiltonian has been derived from the classical expression 
by Wilson and Howardl^V Darling and Dennisonl^'f and more recently by Watsiml**^ Its 
derivation from the classical Hamiltonian is non-trivial because the momentum tenn is 
not conjugate to any cirordinate. The rather complicated original expres.sions for the rota- 
tional Hamiltonian in references(4) and (7) have been considerably simplified by Watson to 
give 
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H = >4 (j„ - n„) 0^ - n^) + I (p.2 + V (Q.)) + u (II-3) 

where J^. 11^^, Pj are the quantum mechanical momentum operators corresponding to the 
previously given classical quantities. The extra term U = 2 - 1/8 is a function only 
of coordinates and not momenta and therefore can be considered as part of the potential 
energy. 

As written, this Hamiltonian is too complicated for any practical calculations of mole- 
cular vibration-rotation energy levels. It is customary to simplify it by expanding the 
potential V and the effective inverse moment of inertia tensor about equilibrium values 
as a Taylor series in the normal coordinates Qj. For reasonably rigid molecules this series 
convei^es fairly rapidly. For light asymmetric top molecules such as H 2 S and HiO, though, 
this series converges relatively slowly because some components of diverge for large values 
of the bending mode normal coordinate. A perturbation calculation can then be made by per- 
forming a series of contact or unitary transformations on the expanded Hamiltonian to make it 
diagonal in vibration quantum number to the required degree of accuracyf'^’'^). This results in 

an effective rotational Hamiltonian for each vibrational state, so that we can write H = H^ + H\- 

V 

where the rotational parameters in Hj^ will depend on the vibrational state of the molecule. 

V 

From now on we will write Hj^ as just H|^. remembering that this Hj^ is for a particular 
vibrational state. The zero order term for Hy is 

= >/i 2 (Pj2 + X. Q.2), (11-4) 

which represents a set of uncoupled harmonic oscillators with energy eigenvalues 

Ey = h^jCnj + Vi). (11-5) 

1 

where the v. are the 3N-6 normal vibrational mode trequencies given by v.= — , the m = 

JJT 

0, 1 , 2, 3, . . . are the vibrational quantum numbers, and h is Planck’s constant. If there are no 
resonances between vibrations one then need only be concerned with diagonalizing the rota- 
tional Hamiltonian in order to calculate the rotational energies. In this work we are mainly 
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concerned with rotational structure within a given vibrational band and will not be further 
concerned with the details of the purely vibrational terms in the Hamiltonian. Approxima- 
tions to Hr will now be discussed for the specific case of an asymmetric rotor which has 
three unequal principal moments of inertia. 

The lowest order approximation for H|^ in a given vibrational state is 

Hr ** = Vi 2^ (II-6) 

where the are constant effective values for the inverse moment of inertia tensor. These 
come from integrals of the over the vibrational wavefunctions. If the molecular 
axes are rotated so they are identical to the principal axes of inertia then is reduced 
to diagonal form giving H^j^y = + Y'Jy^ + Z'}^ where X', Y', Z' = 1/21^, l/2Iy, 

1/21^ and ly, are the principal moments of inertia expressed in proper units. If the 
angular momentum operators J^, Jy, are written in units of h= h/2ir then X', Y', Z' will 
be defined to be in units of wavenumbers (cm“* ) so that all energies in this work can be ex- 
pressed in this convenient unit. Infrared transition frequencies and molecular constants 
will also be given in the cm~* unit, remembering that to obtain the true frequency one must 
multiply the quantity given in wavenumber units by the speed of light. Microwave transi- 
tions and molecular constants will be given in MHz, as is normally done. 

The next order approximation to Hr , derived by Wilson and Howard^”*), is 

where the are given to a very good approximation byi'O) 



The indices i and j run over some set of 3N-6 internal displacement coordinates Rj which are 
used to evaluate where the e stands for evaluating dp^^/dRi at equilib- 

rium. F"l is the inverse force field matrix written in terms of the coordinates Rj. The 
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equation for assumes a harmonic force field and that higher order terms in the ex- 
pansion of are negligible. The dependence of the on the inverse of the force 
constant matrix is physically reasonable since the smaller the force constants of the 
molecule th^ more it will distort at higher rotational states, and therefore the bigger the 
T^pyp . For a general asymmetric rotor there are 2 1 non-zero T^pyg which together take 
on only nine distinct values. The rest of the r^py^ are zero due to symmetry. The nine 
distinct r^py^ are Tg^aaa> ^aa00 ~ ~ ~ ~ 

x,y, z. 

It was later shown by Kivelson and WilsonO^) that by use of the angular momentum 
commutation relations the Hr given by Wilson and Howard can be written in the more 
simple form 


Hr ^ * yj/ + + y. ( 11 - 8 ) 

where ^yyyy - ^yyyy’ ^zzzz ~ ^ZZZZ 


t ^ 

^xxyy 

^xxyy 

+ 

^^xyxy 

t ^ 

^yyzz 

^yyzz 

+ 

2t 

'yzyz 

t _ 

^xxzz 

^xxzz 

+ 

It 

^’xzxz 


X 

Y 

Z 


X' -»■ */4(3rxyxy 
Y' + >4(3ry,y, 
Z' + */4(3 txzxz 


- 2t 

- 2t 

- 2t 


xzxz 

xyxy 

xyxy 


2t ) 
^'yzyz' 

^^xzxz^ 


2t ) 
yzyzf- 


By use of the commutation relations it has been shown that there are really only six inde- 
pendent distortion constants for a general asymmetric rotor and also that the quadratic 
coefficients depend upon the original quartic coefncKnts. 

For a planar asymmetric top molecule such as H2S, Dowling^ has shown that the six 
independent r' coefficients reduce to four independent coefficients. His relations between 
the quartic coefficients are only strictly valid for a harmonic approximation to the force 


9 


field and for the equilibrium coefficients, not vibrationally averaged ones that are observed 
in practice. If the observed energy levels of a molecule are extremely accurate then the use 
of these relations in the Hamiltonian can lead to a poor fit of the data because of these two 
approximations. The relations obtained by Dowling are given in Appendix A. 

It is important to know how many independent sets of constants there are in be- 
cause only the independent constants (which may be combinations of fundamental mole- 
cular constants) can be derived from a molecular spectrum. When accurate rotational energy 
levels of a non-planar asymmetric molecule were fit with Kivelson and Wilson’s Hr, the 
quartic distortion constants were found to be indeterminate, although the observed energy 
levels could be accurately reproduced with the constants obtained in the fit. Watson(l^) 
showed that the matrix elements of the six Kivelson and Wilson’s distortion operators were 
linearly dependent on each other, therefore leading to their indeterminant value in the least- 
squares fitting process. In doing so. Watson developed a general method for developing rota- 
tional Hamiltonians to high orders in angular momentum operators that insures independence 
of the matrix elements of the effective operators. This enables a unique set of distortion con- 
stants to be determined. The Hamiltonians developed by Watson were used in this work since 
they can be easily expanded to the high orders of angular momentum required to fit light 
asymmetric rotors. A summary of Watson’s development of the rotational Hamiltonian will 
now follow. This summary will go into some detail on those fwints important for the inter- 
pretation of the fits performed in this work on H7S rotational data. 

II. B. Watson’s Reduced Rotational Hamiltonian 

Watson developed a procedure to deal with the problem of indeterminate distortion co- 
efficients by systematically making use of the fact that the eigenvalues of Hr do not change 
when it is transformed with some arbitrary unitary operator. By choosing the unitary op- 
erator to be a power series in angular momentum operators, the transformed Hamiltonian 
becomes a power series like the original one but with different coefficients for the angular 
momentum operators. One can then choose the arbitrary parameters specifying the unitary 
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operator to eliminate as many terms as possible from the transformed Hr. resulting in what is 
filled a reduced Hamiltonian. Alternative methods of reducing Hr result in different individ- 
ual coefficients but not in independent combinations of coefficients for operators, which are 
always the same. The remaining terms are then the maximum number that can be determined 
from the rotational energy levels. It should be remembered that since an arbitrary unitary 
transformation of the Hamiltonian does not change the eigenvalues; only those combinations 
ol coefficients not dependent upon the arbitrary parameters in the unitary transformation 
are determinable. 

The steps required to obtain the reduced Hamiltonians used in this work will now be out- 
lined. The following discussion will concentrate only on the effects of the unitary transfor- 
mations on the quadratic and quartic coefficients for simplicity. The transformation of 
higher order coefficients and subsequent reduction cf the Hamiltonian is analogous to that 
of the lower order ones. For detailed discussions of this theory one can refer to papers by 
Watson^^’ 

If there are no vibrational degeneracies or resonances present for a particular vibrational 
state of a molecule, then the effective rotational Hamiltonian, Hr , for that state is a power 
series in components of the total angular momentum operators Jx, Jy, and J^. Watson has 
shown that if one takes the angular momentum commutation relations into account. Hr can 
always be written in the so-called standard form 

"R %!,. hpqr 

where p. q, r = 0. 1 . 2. . . . and the hpqr are constant coefficients. This expression is con- 
strained by the fact that Hr must be Hermitian and time reversal invariant. This requires the 
sum p + q + r to be equal to an even integer and the hpq^ to be real. For molecules with 
orthorhombic point groups a number of the h^^jr shown to equal zero. For a general 
(non-orthorhombic) molecule. Watson showed that the Hr of equation II-9 can be trans- 
formed to look like the Hr of an orthorhombic molecule, which when written up to sixth 
power in angular momentum, is given by 


II 


1 



< 


Hr = XJr! + YJ/ + t T., 

* J ♦.«.« •><.“ + *cal W.* k' * 1.*) 

' *xyz »x“ Jy’ >t * Jx' Jy’ Jx’>- (»-'<» 

The coefficients in this equation are related to Kivelson and Wilson’s parameters by = 
Vi Tg<yj 3 |j. Note that at this poim there are six quartic and ten sextic coefficients. This 
series can be written to higher powers in J if necessary. The transformed Hj^ is given by 

Hr = U-1 Hr U (II-ll) 

where Watson chose 

y = e'^5 giS7 

with 

S 3 ' '111 ('xJy'z Jz^y^x> 

S 5 = '311 <V JyJx * Jz 'y V) 

* '131 ^^x^y^ ^z V’ ^x^ 

+ S„3 JyJ^) (11-13) 

and so on so that U is unitary and time reversal invariant as required. If we write Hr = H 2 + 
H 4 + H 6 + .. . where Hp contains terms to the n*^ power of angular momentum and write 

Hr U = H 2 ' + H 4 ' + Hg' + . . . (II-I4) 

then the transformed Hr terms are given to first order by^-^) 
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0M5) 


Hj' = H2 

H4' = H4 + ilH2,S3l 

H^' = + ilH4,S3l ~ '/iUH2,S3l.S3l + ilH2,S5) 

where [ A, B) = AB - BA. The relations given here for the H only contain terms from 
the transformation that are the same order of magnitude as the untransformed terms, 
any smaller terms have been dropped. The transformation has been carried out using the ap- 
proximation ci^n % 1 + iS^ which is unitary to order S^. This is a relatively accurate proce- 
dure if the parameters snj. S3jj, sj3j, sn3, etc. are chosen to be small enough. When is 
re-written in the standard form previously introduced by using the commutation relations on 

each of the Hp' expressions, terms of lower degree than Hp' may be produced. The resulting 

^ ^ 

terms of similar magnitude can be collected together again to give Hr = H2 + H4 + H6 • 
In order for Hr^ to converge as well as Hj^ the unitary transformation parameters should be 
chosen so that Sjjj ~ T/B and *131 . S311 , Sjj3 ~ ^/B where B, T, and ^ are on the order 
of the untransformed quadratic, quartic, and sextic coefficients in Hr. The transformed 
quadratic and quartic coefficients are, from equation 11-15, 

X = X + 4(Z - Y) Sjii 

Y = Y + 4(X - Z) Sjji 

Z = Z + 4(Y - X) s,ij 

^xx “ ^xx’ ^yy “ ^yy ^zz “ ^zz 
fy, =Ty, + 2(Z-Y) Sj,, 

Txz = T,, + 2(X-Z) s,„ 

T^y = Txy + 2(Y - X)Sjij. (11-16) 

It should be noted that the relative changes in the quadratic coefficients are much smaller 
than those in the quartic coefficients. A typical transformed sextic coefficient is 
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^yyz “ ^yyz 2 (Z Y) Sjjj - ^ (Tyy - Ty^) Sj| j 


(11-17) 


+ 4(Z-Y)s,2, . 

All the transformed sextic coefficients are given by Watson and the higher order contri- 
butions for the quadratic and quartic coefficients are given by Typke(l^). 

A set of quartic coefficient combinations that do not depend, to a very good approxima- 
tion, upon the unitary transformation parameters are 

^xx’^yy’^zz’ ^1 “ ^yz ^xz ^xy’ 

^2 ■■ ^yz YTjj 2 ZYjjy . (11-18) 

These quantities can therefore be chosen as a set of determinable quartic coefficients. This 
choice for the determinable coefficients is only one out of an infinite number of possibilities, 
but any other choice for a set of determinable combinations can be expressed in terms of 
those given above. For the sextic part of Hp there are seven determinable combinations of 
coefficients. In general there are n+1 independent terms in Hp of degree n. 

We now wish to write Hp in a form convenient for calculation by using the operators 
~ lx * 'Jy instead of Jx- Jy. and Jy. This is normally done for the calculation of 
asymmetric rotor energy levels due to the algebraic complexity of the matrix elements of Hp 
if the operators Jx- Jy. and are used. In terms of these operators. Hp has the form 

“ ®200-i^ ®020 ^z^ ®002 '*■ "*■ Y400 A* + ^220 

* '^040>t‘ * ^202l’ * >-> * «To22 Uz' <V + JJ) 

+ (J+^ + + Tqq^ (J^.^ + J_^) + higher order terms. (11-19) 

The relationship between these coefficients and the previously used ones in Hp are given by 
Watson(3). The subscripts of the coefficients refer to the powers of 1 i, and in the 
operator for each coefficient. For example, the quadratic coefficients in this, form are re- 
lated to the original ones by 
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B200 = V2(X + Y) 


Bqo 2 = Vi(X-Y) 

®020 “ Z-‘/2(X + Y). (11-20) 

The transformed quadratic and quartic coefficients in Hr for this form are 


®200 


®200 

+ 

^®002 *111 

®020 

= 

®020 

- 

14Boo 2 Sill 

®002 

- 

®002 

+ 

2 ®020 *1 1 1 

^400 


^400 

- 

2 Bq 02 *111 

^202 

= 

^202 



^220 

= 

T22O 

+ 

12Boo2 *11 

Tq 22 

= 

'*'022 

- 

2 ^020 *1 1 1 

^040 

= 

^040 

- 

IOB002 *111 

^004 

= 

^004 

+ 

®002 * 111 - 


S|]jcan now be chosen to reduce the quartic terms in Hr from six to five. The most 
commonly used reduction, called the AS form, is obtained by setting 


^004 ^ ^^6 

®002 (Y-X) 


( 11 - 22 ) 


where 


= 1/16 ITxx + Tyy - 2TxyI . (11-23) 

This choice for Sj j j makes Tqq^ = 0 so that the operator (J.^^ + ) is dropped from 

Hr and is not included in fitting the observed energy levels to Hr. Similarly, three 
operators in the sextic p.irt of Hr can be dropped with proper choices of Sj3j , S31 1 , and 
S] 13. The resulting Hamiltonian, written in the most commonly used notation is 

Hr^ = H^ + hJ^ + H^ + . . . (11-24) 
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where 


Hj® = ‘/j + Y^) j,* + Y^®)] 


+ Va (X^^ - Y^^) (J+^ + jJ ) 


Hf = -Aj - Ajk h" - Ak j/ - 6j (J+ ■ + i.h 


>/ 6 £k ( 4 ^ 4^5 


Hf = Hj + Hjk J,' ^ Hkj j/ + Hk j/ 

+ hj (j+^ -*■ J.4 + ‘/ihj,^ (V + 4^)+ (4^ ■•• J- ) Jz^i 

+ Vih,^ Uz^< 4 ^ 4 ^) + ( 4 ^ + L^) 4 ^) 


nf ~ Lj xi* ■*■ 4 jk Jx 4 ■*■ 4k J- 4 ^kkj 4 ■*■ ^k 4 


|8 4 . I 1 ^ 4 - I ... 1 ^ I ^ + I j 6 ^ . ,8 


+ «J /{j+^ + ij) + ‘/^CjK i.^ 14^4^ + J.^)+(J+^ + 4^)4^ 


+ Vii^^j }} (4^ <4^ + J/) + (4^ + 4‘)4"1 


2 . . 2, , 4, 


+ Viif^ (4^(4^ + 4^) + (4^ 4^)Jz^l- 


This is the most commonly used reduced Hamiltonian because of its relative ease of calcula- 
tion, which will be discussed in Section II-C. 

Another possible reduction of Hr. that is sometimes used is called the NS form.(^) re- 
sults from setting 


022 


‘111 


2 B, 


020 


Z - (X Y) ’ 


(U-25) 


where 


Rj = 1/8 (Tjjjj Tyy 2 Tj^j, + ^yz)- 


(H-26) 
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? 2 2 

This choice for Sj jj sets Tq22 = 0 so that the term with the operator h [ 3 jT (J+ + ) 

+ ) is dropped fiom Hj^. The equations for Sj^i , s^n , and S||3 for both 

the AS and NS reductions are giver by Watson^^). The NS reduced Hamiltonian can be 
written 


hNS ^ j^NS ^ j^NS ^ j^NS ^ 


(11-27) 


where 



+ iz‘"" - + 


/NS 


NS 


,NS 


i, 


+ '/4 - Y^^) (J+^ + jJ) 


Hf = ^Dj I* ~ D;k !,* 

+ 6j J-(J+“ + J_‘) + 52(J/ + J_'*) 

H6* = h; J*' + HV I* + H|u r + Hjj J/ 

4 -> -> 1 4 4 

+ h, r (J+ + J-") + ^2 V 0 + + j_ ) 

+ h, (J+ ’ + J_^) 



Lj Ljjj( ^ Ljk I 

+ Lk + Cl r (J+‘ + 

+ ?3 i" (J+^ + + ?4 (J+^ 



+ c. r 


+ 


L' J-J^ 
^KKJ < •'z 


4 4 

(J+ + J_ ) 


These two forms for Hr were used in this work. The H2S pure rotational spectrum in 
the microwave spectral region was fit by Helniinger et al.f using the AS from written in 
terms of powers of the operator J - Jy” instead of J 4“, which results in some slightly 
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different operators in and The microwave data as well as rotational line fret]uencies 
in the 2»^2« *^3 hands of H 2 S have been fit by Gillis and Edwards^ to a reduced 

Hamiltonian developed by Typke^ which is very similar to the NS form. Typke also em- 
ployed pc)wers of the operator Jy“ instead of to write Hr , which again leads to 
slightly different higher order operators in Hr. 

II. ('. C'u culation of the Rotational Energy Levels 

In this section the calculation of energy levels from Hr will be discussed as well as the 
molecular symmetries which affect this calculation. Unlike a symmetric rotor, the energy 
levels of an asymmetric rvitor. given by Hr I = Ej I >. cannot be written in closed 
form, even in the rigid rotor appi\)ximation. To obtain the Ej. the Hamiltonian matrix is 
computed in some basis set and is then numerically diagonalized with a computer. 

A amvenient basis set consistsof the symmetric top wavefunclions U.K>. 1 he M de- 
pendence of these functions, where M represents the quantum number corresponding to the 
2J + I orientations of the total angular momentum in the space fixed coordinate system, is not 
shown because field free rotational energies do not depend upon it. Since the Hamiltonians 
we have discussed arc diagonal in J. the only matrix elements required for the calculation of 
the Ej am < J. K’ i Hr | J. K > . The Hamiltonian matrix is block di;igonal in J and can 
therefore be diagimali/ed separately for each J value. The only non-/en> matrix elements for 
a rigid rotor are 

< J. K I Hrigjjl J . K > = ‘ i (X + Y) J (J + I ) + |Z ' : (X + Y) | K" 

<J.K±:|Hrij,jj|J.K> = U(X Y)||J(J + 1) K(Kil)|* 

|J(J+ I) (K 1 l)(K ± :)I| - . (11-28) 

so that H^j^j^j is in tri-diagonal form. 

The numerical diagonalization of Hr is facilitated by taking into account the symme- 
tries of asymmetric rotor wavefunctions. which belong to the four-group V(a. b. c) defined 
by the three rotation operators t',*'. and the identity o|U'rator E. The labels a. b. 


L 
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c denote the molecular axes and are defined by the requirement that the principal mo- 
ments of inertia satisfy the relation As shown in Figure U-1 . the planar HiS 

molecule lies in the ah plane, with the b axis having 2 fold symmetry. C-,**. for example, 
represents a ISO® rotation of the molecule about its a axis. The a. b, c axes of the molecule 
can be identified with the x. y. / molecular axes used up to now in the six ways slmwn in 



0 ^ 92 " 

c axis 1 to page 

r„.s « I..UA 

Figure 11-i : HiS (kHMUcliy, C'cnlcr of Mass at Origin 

Table 11-1. The character table for V(a, b. c) is given in Table 11-2, with a + 1 representing 
a symmetric wavefunction for the given operation and a -1 an antisymmetric wavefunction. 


tabu: 11-1 

t'oonlinate Repix'si'ntations 



l*^ 

11^ 

llF 

1^’ 

11^’ 

Isl'^ 

X 

h 

L* 

a 

c 

a 

b 

y 

c 

a 

b 

b 


a 

/ 

a 

b 

c 

a 

b 

c 
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TABLE 11-2 

Character Table for V(a, b, c) 


Symnictry Species 


Operation 


+ 

1 

Type 

E 

r ^ 

c ^ 

r c 

ee 

A 

1 

1 

1 

1 

eo 

Ba 

1 

1 

-1 

-1 

oo 

Bb 

1 

-1 

1 

-1 

oe 

Be 

1 

-1 

-1 

1 


It was observed by Wangl that the I J. K > symmetric roto*- basis functions do not trans- 
form as a representation of the four-group V(a, b, c), while the linear combinations of these 
functions 

I J,0+> = I J,0> 

I J,K+> = 1/v^ {l J,K> + i J,-K>( K>0 

|J,K"> = 1/n/2 {|J,K> - 1J,-K>( K>0 (11-29) 

do have this property. For the 1 J, K > basis functions 1 J, K > = J(J + 1 ) I J. K > and 
I J. K > = K IJ, K >. where the z subscript refers to the molecule fixed z axis. If 
Hr is set up in this basis, it can be factorized for each J (except J = 0 or 1), into four inde- 
pendent tri-diagonal submatrices. This can be done because the matrix elements connecting 
I J, K"*” > and 1 J, K“ > vanish as well as those connecting even and odd K. For a given J 
we have 

h{^ = £■*■ + + E" + O' (11-30) 

where E and O refer to the evenness or oddness of the K values in the matrix elements and 
the + and - refer to the matrices with only | J, K'*’ > or | J, K" > matrix elements. This 
factorization of Hj^ is called the Wang transformation. Although the Wang transformation 
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is not necessary for the computation of the Ej, it does reduce the size of the matrices that 
must be diagonalized from 2J + 1 by 2J + I to approximately J/2 by J/2. which can be a signifi- 
cant reduction in computer memory, especially for high J values. The diagonalization of the 
four Wang submatrices gives 2J + 1 energy levels for each J value, which in general are non- 
degenerate. 

The behavior of the energy levels can be described using the parameter k, where k = 

(2B- A-C)/(A~C) and A, B, C = l/2Ia, 1/21^,, 1/21(. in the proper units. If the moments of 

inertia of an asymmetric top are varied from the limiting prolate symmetric top (B = C) to 

the limiting oblate symmetric top (B = A), k will vary from -I to +1 , the most asymmetric 

top having k= 0. For H 2 S k +0.5 and thus it is a fairly asymmetric rotor. The asymmetric 

rotor energy levels are generally designated by . where Kj is the K value of the limiting 

symmetric prolate top level (k = -1) with which the asymmetric top energy level connects 

and Kg is the K value of the limiting symmetric oblate top level (k = +1) with which this 

same asymmetric top energy level connects. For a given J the asymmetric top energy levels 

do not cross, except in extreme cases which were not encountered in this work,^^^) making 

this labeling procedure possible. The identification of Kj. Kg values with the eigenvalues of 

the Wang submatrices for a given J can be found in a number of references^ and is 

easily incorporated into computer programs that calculate asymmetric rotor energy levels. 

This notation also conveniently gives the symmetry of the rotational wavefunctiors 

c 

from the parity of K, and Kg, as shown in Table 11-2. 

For the previously discussed AS form of the Hamiltonian the only non-zero matrix 
elements are 

<JK|Hr^|JK> and < J K±2 1 I J K > 

so that the calculation of the energy levels is identical to that for a rigid rotor, the only 
modification being the addition of higher degree terms to each matrix element. For the NS 
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Hamiltonian there are non-zero matrix elements for AK = 0, ±2, ±4, and ±6. Still, the 
Hamiltonian matrix can be Wang factorized and line identifications assigned as for the rigid 
rotor. The required matrix elements for both the AS and NS Hamiltonians have been 
tabulated by KirshneK^O. 


11. D. Vibration-Rotation Line Strengths of an Asymmetric Top Molecule 

To interpret the vibration-rotation spectrum of an asymmetric top molecule a knowl- 
edge of the transition strengths (and selection rules) as well as the molecular energy levels 
is required. If the approximation is made that there is no interaction between vibration 
and rotation, the calculation of the relative line strengths within a vibration-rotation band 
is straightforward. The relative line strengths obtained using this approximation, which is 
called the rigid rotor approximation, often agree with experimentally observed line 
strengths to 30% or better. The band of H 2 S. when observed under low resolu- 
tion {= 5 cm’*)( 22 ) exhibits an intensity profile very different from a rigid rotor pro- 
file. The H 2 S intensity profile, which is characterized by the P branch region being about 
I/S as intense as the R branch region of the band, was studied by Emerson and Eggers,(l) 
who showed that the intensity anomaly is caused by centrifugal distortion of the nonrigid 
H 2 S molecule. In this section we will discuss the fundamental relations for the vibration- 
rotation line strengths and also the rigid rotor approximation for these strengths. The 
approach taken by Emerson and Eggers to include some of the effects of centrifugal dis- 
tortion in the line strengths will also be summarized. 

The line strength is defined to be the integral of the absorption coefficient, k(v), 
over the line. 


S = 


/ 


k(i») d»». 


(II-3I) 


over the line 

In terms of fundamental molecular quantities, the line strength is given by^^^) 


S = 


8ir3|; Ngj e’^'^‘‘’^ 


3hc 0 


(1 - e-bWkl) K^>|2 


(11-32) 
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where N is the number of molecules/cm^-atm, v is the frequency of the transition, E, is 
the lower state energy of the transition, gj is the statistical weight or degeneracy of the 
lower state energy level E,, Q is the vibration-rotation partition function, and K/i>l2 is 
the square of the transition dipole moment matrix element. The partition function is 
given by 

Q = Sgj (11-33) 

j ^ 

where Ej are the molecular vibration-rotation energy levels. 

The square of the dipole moment matrix element is given by 

I<M>12 = _E i<ilMFlf>|2 (H-34) 

F=X, Y , Z 

where /ip is the component of the dipole moment along the X, Y, or Z space-fixed axes, 
and ii> and If > are the initial and final states of the transition. In the absence of any 
external fields the X, Y, and Z components of <fi> all have the same value so that only 
one component (the Z component for example) needs to be evaluated, giving 

K/i>|2 = 31<ilMzlf>l^- UI-35) 

The component of the dipole moment along the space-fixed Z axis can be related to the 
dipole moment relative to the molecule-fixed axes with the direction cosines, </>pg by 

‘ 

where g is summed over the three molecule-fixed axes, a, b, and c. 

The »»2 vibration of H 2 S studied in this work is illustrated in Fig. 11-2. For this 
vibration and Mc not change, so that only needs to be included in equation 
11-36 for P2- permanent dipole moment of H 2 S is also along the b axis, and gives 
rise to the pure rotational ^ectrum of H 2 S. Since H 2 S is a planar molecule = 0. 

The only time fig not zero is during the asymmetric vibration can now 

write the dipole matrix element for the 1^2 vibration as 

l<pi>l2 = 31<i|^b (H-37) 

At this point, the approximation that there is no interaction between vibration and 
rotation can be made. This means that does not depend on rotation coordinates and 
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c axis i to page 

Figure 11-2: The V 2 Vibration of H 2 S 

that I i> and 1 f> can be separated into vibrational and rotational parts | V> and | R>. so 
that now |<;x>|2 can be written as 

\<H>C- = 3|<Ril^ZblRf>l^ l<ViliiblVf>|2 (11-38) 

where 

|i> = |Ri>|Vj> and | f>= ( Rf > | Vf >. 

This is called the rigid rotor approximation in this work, although the dependence of | Rj> 
and I Rf> on the average centrifugal distortion within each vibrational state is included 
when evaluating <Rj li^^b I v'omparison to the observed V 2 band line 

strengths. The dependence of <Rj|(^2;b !^f^ vibrationally averaged centrifugal distor- 

tion in each state is obtained by using one of Watson’s high order reduced Hamiltonians to 
generate | Rj> and | Rf>. In the usual rigid rotor approximation, | Rj> and | Rf> are generated 
from the rigid rotor Hamiltonian given in equation 11-6, using the same rotational constants 
for both the lower and upper stales. The direction cosine matrix elements in the |JKM> rep- 
resentation are given in many placesH^'-^^ and are used to evaluate <Rjlv?2b l*^f> s'"ce the 
asymmetric lop rotational eigenfunctions are calculated in terms of |JKM> basis functions. 

In the rigid rotor approximation the total band intensity of the vibration-rotation band 
is determined by the <Vj|/jH 5 iVf> matrix element. If the dipole moment is expanded about 
the equilibrium value as a Taylor series in the normal coordinates Qj, we can approximate 
Mb by 
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1 


Mb 


‘=l\3Qi/Qi=0 


( 11 - 39 ) 


where is the permanent dipole moment of H 2 S. For the 1^2 band the vibration quantum 
number n 2 given in equation 11-5 changes 0 to 1 . Only Q 2 is non-zero for this vibration, giving 


<VilMblVr> = — <VilQ2lVf>. 
3Q2 


( 11 - 40 ) 


since <V||ju^ 1 Vj-> = 0 for a rigid rotor. 

The selection rules for changes in *^c 1^2 vibration can be determined 

using Table 11-2. For the direction cosine matrix element to be non-zero, the product of 
the characters of the representations of the initial and final rotational wave functions and 
of the dipole moment change must be +1 for alt of the group operators. For the »'2 vibr 
tion, the change in the dipole moment belongs to the type Bj, symmetry species, so that 
the selection rules for Kg and are 

ee ► 00 , eo oe. 

This constrains the changes in Kg and to be ±1, ±3. ±5, . . . Generally transitions with 
changes of ±1 in Kg and K^. are the strongest. The dipole selection rule AJ = 0, ±1 
also holds. The vibration-rotation transitions will be denoted by JV' v* - Jif ic in this 
work, with the primed upper state values of J. Kg and K^. given first, followed by the 
lower state values of these numbers. It should be noted that gj, the statistical weight 
factor, is three times larger for the eo and oe levels than for the 00 and ee levels of the 
same J due to the nuclear spin statistics of the two identical protons n HiS. 

To take into account the effects of centrigual distortion on the line strengths, 
Fmerson and Fgger^*^ used vibrational wavefunctions having coordinates which include 
centrifugal distortion displacements of the atoms from their rest equilibrium positions. 
They calculated the shift in the Qi ^md Q 2 normal coordinates from their rest equilibrium 
positions in both the upper and lower rotational states of each transition. This was done us- 
ing classical expressions for the molecular distortion, which depend upon the rotational 
energy and the square of the angular momentum about the a. b. and c axes. The rotational 
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energies and angular momenta were determined by Emerson and Eggers using the rigid rotor 
Hamiltonian given in equation II-6. Using these new shifted coordinates in harmonic oscilla- 
tor vibrational wavefunctions. they evaluated <Vj|jat)| Vf> and obtained a first order correc- 
tion to the rigid rotor line strengths, 

S = Srigidd +76pg/OMb/dQ2))2- (11-41) 

where 

6 = AQ^ - AQij 

is the difference between the centrifugal distortion of the normal coordinate Q 2 in the 
final state and the initial state, and 7 is a factor that depends upon the 1^2 vibrational 
frequency. 

By choosing and bii^/dQ 2 to have the same sign Emerson and Eggers were able to 
obtain qualitative agreement between the observed spectrum and a spectrum calculated 
using their strength correction term, showing that centrifugal distortion was indeed respon- 
sible for the observed intensity anomalies in this band. Since 6 has the same magnitude 
but a different sign for the transition J'j^' j^' ^ compared to its reverse Jk K 

J'k' k' • branch of the band is enhanced while the corresponding branch on the other 
side of the band center loses intensity. The intensity perturbation is also dependent on 
the ratio /ag/(9/ijj/9Q2), which is very large for H 2 S compared to other molecules because 
dn^/dQ2 's so small. The small size of 9^b/^Q2 illustrated by the fact that the 1^2 
of H 2 S is roughly 130 times weaker than the V 2 band of H 2 O. 

The magnitude of the line strength perturbation is denoted by the F factor, which is 
the ratio of the observed line strength to the calculated rigid rotor line strength. For the 
H 2 S lines observed in this work, the F factors were experimentally tound to vary from 
about 0.06 in the branch to greater than 4 in the qRP 

branch. (The branches 

are designated using the convention AJ^*’^*^c.) A discussion of the observed F fac- 
tors can be found in Ch. VI. 

Perturbations of HiO line strengths by centrifugal distortion have also been observed 
and several theoretical formulations have been devised to take these effects into 
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account .( 24 - 26 ) Although the F factors for H2O lines are much closer to unity than those 
for H2S for the same transitions, lines in the 1^2 of H2O involving very high energy 
rotational states have been observed, with F factors ranging from 0.07 to 3 . Generally, 
the above referenced theoretical strength calculations were able, by least squares adjust- 
ments of free parameters, to fit the observed line strengths in the t^2 band of H2O to 
\07r-207r. which is approximately the accuracy of the H2O line strength measurements. 

It is hoped that the H2S vj band line strengths measured in this work will be useful 
for a stringent test of existing theories developed for H2O. The H2S line strengths should 
be ideal for this purpose due to the wide range of F factors observed and to the high 
accuracy of the line strength measurements made here. 
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CHAPTER III 


DIODE LASER SPECTROMETER AND DATA ACQUISITION 
III. A. Introduction 

The tunable diode laser spectrometer, constructed as part of this work, was used to 
measure all of the line strengths and the majority of the line frequencies of H 2 S reported 
here. The use of semiconductor diode lasers allows almost Doppler limited spectra to be 
taken in continuous segments about 0.5 cm*^ to 1 cm'^ wide. Diode laser fabrication 
technology has been developed to the point where spectral studies extending over 1 00 an** 
are practical using a single diode laser, so that rather complete studies of even light asym- 
metric top absorption bands can be made with these lasers. 

The laser spectrometer was constructed so to meet several experimental requirements. 
These include: monitoring the relative frequency tuning of the laser simultaneously with 
the gas absorption signal; having a relative frequency calibration capability close to the 
accuracy of the best frequency standards (®»0.0004cm** ); being able to produce sample 
pathlengths from several meters to greater than SO meters in order to detect the rather 
weak and variable strength H 2 S absorption lines at low pressure; and recording the data in 
digital form because of its large volume and information content. In this chapter the laser 
spectrometer will be described as well as procedures used to obtain a spectrum with this 
instrument. Several examples of spectra will be shown and discussed to illustrate the 
capabilities of this technique and to show typical characteristics of the spectrum of H 2 S. 

A number of the line frequencies in the R branch region of the 1^2 ^2^ 

were measured with a high resolution (0.04 cm** ) grating spectrometer at the National 
Bureau of Standards (NBS).(27) ^ short summary of the experimental work performed 
with this instrument will also be given. 
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III. B. Experimental Apparatus 

The Pb-salt (Pb Si_x Se^. Pb|.x SHx Se) tunable diode lasers used in this work were 
fabricated by Laser Analytics Inc. Laser action in these semiconductors results from 
stimulated emission across the energy gap between the conduction and valence band. A 
population inversion is achieved by applying a forward bips current to the diode and 
thereby injecting charge carriers across the p-n junction. The recombination of these carriers 
provides the gain mechanism for laser action. The end faces of the laser crystal, made by 
cleaving along natural crystal planes, form the laser resonator. Pb-salt diode laser crystals 
are typically 400pm long and have cross sections of 200pm by 200pm. Pb S].x SCx 
lasers can be made to operate from about 2S00cm*l to 1200cm~^ and Pbj.x Snx Sc lasers 
from 1200cm** to 300cm** by controlling their chemical composition (x) during 
fabrication. 

The output frequency of any one diode laser can be tuned 20cm** to 200cm** by 
changing the laser temperature from lOK to greater than 60 K. The laser can be fine-tuned by 
varying its bias current, which changes the amount of l^R heating of the laser and its electrical 
contacts, and therefore the operating temperature. While the overall tuning of the laser comes 
from the variation of the band gap energy with temperature, the fine-tuning within one mode 
is dominated by cavity mode shifts due to the decrease in the index of refraction of the laser 
material with increasing temperature. 

Most diode lasers oscillate in several simultaneous modes separated in frequency by the 
cavity mode spacing of several wavenumbers. A single mode will usually continuously tune 
0.5cm'* to I cm**. Diode laser linewidths have been measured by a number of worker^ 28, 33) 
and can vary from less than I MHz to greater than 20 MHz FWHH. The diode laser linewidth 
will be discussed in more detail in Ch. IV. Single mode power levels range from several tens of 
microwatts to greater than one milliwatt for diode lasers. A PbSSe laser was used for the 
H2S R branch region measurements and a PbSnSe laser for the P branch region measurements. 

A second PbSSe laser was used for the HI study. Several reviews of Pb-salt tunable diode 
lasers can be consulted for details of their fabrication and operation.^ ^-^5) 
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A layout of the diode laser spectrometer is given in Fig. Ill* I. The laser is kept in a 
Cryogenics Technology Inc. closed cycle gaseous helium refrigerator that has an ultimate 
temperature of about lOK. The laser is mounted on a cold ; nger with a large thermal 
inertia that fs vibration isolated from the rest of the refrigerator. The heat sink tempera- 
ture can be varied with a small heater that is feedback controlled using a silicon diode 
temperature sensor. The heat sink can be varied in temperature from about lOK to 
greater than lOOK and can be held stable to less than 0.3 mK. The refrigerator and tem- 
perature control system are available commercially from Laser Analytics Inc. 

The laser radiation is collected with a f/1 KRS-5 lens. The beam is then focused 
onto the entrance slit of a 50cm focal length grating monochromator to select a single laser 
mode. A 400 Hz tuning fork chopper at the entrance slit modulates the beam for phase- 
sensitive AC detection. The peak resolution of the grating monochromator is about 0.25 cm’* 
and the wavenumber accuracy is about 0.5 to 1 cm’*. 

The radiation emerging from the monochromator is collimated with a 5cm diameter, 
30cm focal length off-axis parabola. This results in a maximum beam width of about 2.5 cm. 
The collimated beam is split approximately fifty-fifty with an uncoated piece of germanium, 
wedged slightly so not to produce interference fringes. One beam then passes through a 7.b5 
cm long solid germanium etalon for relative frequency calibration and is focused onto a I mm 
square HgCdTe liquid nitrogen cooled detector with a 5 cm diameter. 15 cm focal length off- 
axis parabola. The second beam traverses a 1 meter base path multiple pass absorption cell 
(White cell^^^*) and is focused onto a second 1 mm square HgCdTe detector with a l.'icm 
focal length lens. The two lenses marked L2 and L3 in Fig. lll-l , as well as the White cell 
windows, are made of Bap2. Lens L2 has a 27 cm focal length and focuses the collimated 
beam to a point coincident with the front surface of mirror M| . 

The optical design of a White cell has been discussed in many referenced and will not 
be repeated here in detail. Spherical mirrors M| ( 10 X 6cm), M2 (9 X 7cm), and M3 (9 X 7 
cm) all have a 100cm radius of curvature. M] is separated from M2 and M3 by their common 
radius of curvature so that the beam is always focused when it hits M], irrespective of the 
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Figure III-l . Diode Laser Spectrometer 





number of passes, and no radiation is tost dvie to v^netting. The top half of M| is notched in 
the usual manner for passage of the entrance and exit beams. By adjusting M 2 about a vertical 
axis, the beam can be made to traverse the celt from 4 to approximately 60 times in multiples 
of 4 passes, the only power loss ct)ming from mirror reflection losses or from absorption by 
the gas in the cell. All three mirrors are gold coaled for high reflectivity throughout the infra- 
red. The White cell housing is a inch ID glass pipe. Aluminum end caps sealed to the pipe 
with "O’* rings hold the three mirrors. A rotary vacuum feedthrough near mirror Mt abows 
adjustment of the number of passes while the cell is under vacuum. The White cell path- 
length is given by 

L * 90.855 • N 28,840cm (lU-1) 

where N is the number of pas.ses the radiation makes in the cell. This equation for the path- 
length was obtained from careful measurements of the cell dimensions and is accurate to 
about 0.3%. The cell tcmixrature is monitored with a thermoavupic attached to the cell wall. 

The cell is evacuated with a 21 cfm mechanical pump via a liquid nitixigen trap. This 
pump can evacuate the cell to several millitorr in a few minutes. A pure gas sample is achieved 
by tilling the cell once with HiS. pumping the cell to several millitorr, and then refilling it a 
second time with H 2 S. Three MKS Baratron capacitance pressure gai^es with full .scale ranges 
of 1000 torr, 10 torr, and I torr monitored the cell pres.sure. Over the operating range of each 
gauge the accuracy is stated to be 0.08% of the pressure reading, if the gauge is zeroed properly. 
To get an accurate zero setting the gauges can be valvevi off fnun the large White cell volume 
and punqx'd with a diffusion pump to obtain a high vacuum. Intercompariams of the three 
gauges indicate that their accuracy is at least a factor of two wt>rse than the stated 0.08% of 
reading figure. Care was taken to not allow the I torr gauge pressure to go above about 10 
torr without it being re-zeroed again, a procedure that takes siweral hours. If this gauge is let 
up to atmospheric pres.sure its zero reading usually changes a significant amount. In any case, 
leaks in the White cell (or outgassing of cell components) suggests u maximum uncertainty of 
1% for the H 2 S gas pressure, well above the inherent uncertainty of the pressure gauges. The 
I torr full scale gauge, used for almost all the pressure measurements here, was calibrated 


to a quoted accuracy of less than 0.08% by the manufacturer just before the start of 
this work. 

The H 2 S gas was supplied by Scientific Gas Products and had a stated purity of 
99.90%. Two lecture bottles ot H 2 S were used in the course of the experiment. One H 2 S 
line strength was measured using samples of H 2 S gas from each lecture bottle. Both meas- 
urements gave the same line strength within the precision of the experiment. 

The solid germanium etalon has a free spectral range of about 0.01 6 cm~^ and pro- 
vides relative frequency calibration. Unfortunately, the index of refraction of germanium 
is very sensitive to temperature, resulting in laige shifts of the etalon fringe peaks as the 
room temperature flucuates. M. El-Sherbiny et al.(^^) r«luced the magnitude of these 
shifts by placing their etalon in a copper cylinder whose temperature was stabilized by 
flowing water from a temperature controller through the cylinder walls. They achieved 
a stabilization of better than 0.008 K. for time periods up to one hour, which corresponds to 
a fringe drift of O.OOOScm*^ since for a germanium etalon 

dvpeak 

1.05 X IO-^v(cm*VK). (lll--2i 

dT 

di/peak 

The dependent of on the thermal expansion of the etalon is negligible. 

dT 

We have followed this approach by placing the etalon in a copper cylinder as shown 
in Fig. III-2, where the top half of the cylinder has been removed for viewing. Water 
from a Thermomix 1480 temperature controlled water bath is circulated at a rate of 
about 10 liters per minute through four channels drilled in the copper. The inner surface 
of the cylinder was hand lapped to achieve good thermal contact between the copper and 
the etalon. The water bath temperature was set at 303. 25K. about 7K above room tcmpe^ 
ature and was kept constant to about ±0.004 K. At 1 200cm** " » 0.1 2cm**/K so the 

±0,004K variation in temperature translates into ±0.0005cnr*. The copper cylinder is rela- 
tively massive, providing thermal inertia to help keep the frequency drift of the fringes 
below ±0.0005 enr* over the time scale of a diode laser scan (2 to 20 minutes). An extremely 
slow tuning, and therefore very stable, diode laser was available to test the etalon 
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stability. Over a peiiod ot one naif hour the Irequeney separation ol a fringe peak and a gas 
absorption line next to it was measured six times. The average fringe peak drift between meas- 
urements was found to be only 0.0001 cnr'. and the maximum drift over the whole one-half 
hour time period was O.OOOZcnr*. an acceptable amount considering other error sources in 
the measurements. 

A block diagram of the diode laser spectrometer electronics is shown in Fig. 111-3. 

The laser is tuned by varying the injection current with a Laser Analytics Model LCM cur- 
rent controller. A 50Hz to 1 kHz AC current ramp can be applied to the laser for real 
time display of the detector signals on an oscilloscope. Digitized data is recorded by ex- 
citing the laser with a slowly ramped DC current, mechanically chopping the radiation, 
and using lock-in amplifiers for AC phase sensitive detection. A switching unit allows 
the outputs of the lock-in amplifiers to be digitized sequentially at a combined rate of 
up to 70 samples per second. The readings from the 3-1/2 digit voltmeter are stored in a 
Hewlett Packard 9825 computer during the spectral scan. The digitized absorption 
spectrum and the etalon fringes can then be examined on a CRT display and stored on 
magnetic tape. Housekeeping information such as gas pressure, temperature, cell path- 
length, laser temperature, etc. are input by hand to the computer and are also stored on 
magnetic tape. 

Each spectrum consisted of 4000 points for both the absorption spectrum ar.'^ etalon 
fringes. Although this many points were not required for most of the shorter scans, keeping 
this number constant made operation of the data recording and analysis programs simple. 

Only the scan time desired (always 1 10 seconds or greater) and housekeeping information 
(if It changed) need be entered to the computer to start a scan. The outputs of the lock-in 
amplifiers could also be recorded on an analog X-Y-Y recorder when optimizing the experi- 
ment parameters before digitization. 

111. C. Operation of the Diode Laser Spectrometer and Sample Spectra 

The laser spectrometer optics were aligned with a He-Ne laser and a white light source. 
The etalon was purpcisely mis-aligned to prevent hack rellections of the laser radiation off 
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the etalon from re-ontcring the exit slit of the grating monochromator and feeding back to 
the la:.er. Both detector dewars were also tilted slightly so rellections of the laser radiation 
from the detectors cx'mld also no* travel back to the diode laser. Without these mis-align’ients 
composite optical cavities are formed that can set up laser modes that may broaden the ap- 
parent linewidtb of the laser and increase the noise level.^’*^' 

Final alignment of the laser spectrometer is performed by observing the detector signals 
from the diode laser radiation on a dual trace oscillosa^pe. By reducing the grating mono- 
chromator slitwidths and adjusting lens L|. the image from a single laser mode is focused at 
the center of the exit slit. The entrance slit is then opened very wide to prevent reflections 
off the slit jaws from returning to the laser. The flat mirror before the etalon and the para- 
bola [*2 ‘'djusted to obtain etalon fringes with the highest finesse possible, ('are was taken 
si> that these adjustments did not a*low back reflections from the etalon to re-enter the grat- 
ing monochromator or tilt the etalon loo much relative to the a>llimaled beam. The finesse 
of the etalon fringes was also improved by masking off about 1 3 of the 1 inch diameter 
etalon. This improvement anild be a consrt|iience of poor optical quality of part of the eta- 
lon or may be from the effects of unwanted reflected radiation in the optical system. Also, 
a .^'4 inch diameter iris diaphragm placed just before lens L> increaseil the absorption beam 
amplitude stability and prevented, a slight amount of radiation reflected from the etalon beam 
from reaching the absorption beam detector. This iris diaphragm decreased the absorption 
beam signal very little. 

The previous alignment procedures did not need to be repealed very often, l ach new 
laser mode could usu.illy be aligned by adjust mg lens I |. the grating angle, lens I the White 
cell mirror M>. and the parabola IN With caie. amplitude jiller of the laser signal caused by 
the refrigerator vibrations could be kept below O.IT. Before recording the spectrum the grat- 
ing was scanned to be sure that the other laser modes were sufficiently far away in frequency 
from the mode being used The ability to observe the absorption signal and the etalon 
fringes together m real lime during the alignment process was invaluable. Very often the 
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gas absorption spectrum would appear to be of high quality while the etalon fringes would 
indicate a inode hop or frequency instability. 

After these adjustmentswere made, the chopper was turned on ;>nd the laser slowly 
ramped with a IK' current for detection with the lock-in amplifiers. A time constant of 0.1 
second on the outputs of the lock-in amplifiers was always used. To be certain tha* the ab- 
sorption lines exhibited no distortion, the laser scan rate was set so that the full width at half 
height of a line was scanned in a time period approximately equal to twenty time constants 
(2 seconds). For the strength measurements the gas in the White cell was allowed to come to 
equilibrium for a period of at least 1 5 minutes after filling. Whenever (xissible the gas pres- 
sure was kept below 0.3 torr and the pathlength chosen to give a peak absorption ot to 
70' (. The optimum peak absorption for determining line strengths fuiin the peak absorption 
coefficient, as was done here, is Slit widths equivalent to about 0.4cnr* were used 

to separate the laser modes as completelv as jxissible. Usually 3 digitized scans each ot an ab- 
sorption line at 3 different pathlengths were made. A shutter was placed in the laser beam lor 
a short period during each scan to determine the 100'.' absorption level. Also, in most cases 
a scan of each mode was made with at least one abs».>rption line saturated (i.e. absorbing lOO'i 
of the incident radiation at line center) to test for mode purity. To check for non-linearities, etc. 
Ill the apparatus a set of scans was taken with a constant pressure and 13 different pathlengths. 
.\nother set of scans with a constant pathlength and It different gas pressures was also re- 
corded lor the same reason. Fhe line s' -ngths obtained from these scans using the methods 
described in t'h. IV showed no evidence of non-linearities in the intensity scale or in the pres- 
sure readings greater than the measurement precision of alx>ut 1 '.?. 

lo measure the H-»S line frequencies, a calibration gas was mixed with the H 2 S in the 
White cell, ( alibralion gases included t'C) 2 . N 2 O, and NH 3 (see t'h. IV toi ..etails on the cali- 
bration lines). Usuallv three, and sometimes up to seven, scans were made for each I'cq'ieiKN 
measurement. For these scans the grating monochromator slits were opened much wider 
than in the line strength scans to allow the I cm*' wide modes to be recorded w ithout large 
changes in laser |x>wer across the >pectrum 
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The general procedure in the search for good laser modes was to start with low laser tem- 
peratures and currents and then slowly increase both of these parameters to find new modes. 
This was done since las.r failure is much more likely at high temperatures and currents. By 
systematically varying the laser temperature, current, and the grating angle all useable modes 
were found, so as to achieve as cx>mplete a spectral coverage of the M 2 S spectrum as possible. 
Also the laser (heat sink) temperature was kept as high as ^lossible and the injection current as 
low as possible when working with a given mode. This usually resulted in more stable (but 
weaker) laser emission and also reduced the number of modes at other frequencies. 

A number of diode laser spectra of HsS are shown in Figs. 111-4 to 111-18. Most of the 
spectra shown are calibration scans since they are longer than the strength scans and better 
illustrate the spectrum of IhS. Most of the strength scans include only a single line over a 
spectral region of five or six fringes. Lines not identified in these spectra are due to calibra- 
tion gas absorption lines other than the calibration line itself. 

Figs. 111-4 through 111-10 show spectra containing two or more calibration lines per 
mode. These spectra and others were used to test the accuracy of the relative frequency cal- 
ibrations performed here, which are discussed in C'h. IV. High signal-to-noise ratios were 
obtained for most of the spectra, even at pathlengths as long as 48 meters as can be seen in 
Fig. 111-4. The baseline variations seen in the spectra are due to laser pciwer variations, the 
slit function of the grating monochromator, and other causes. It should be remembered that 
as the laser is scanned its image is translated in the focal plane of the monochromator, thereby 
producing baseline variations. 

These figures show that many modes of the laser used for the P branch measure- 
ments were at least 1 cm** long. The mode m Fig. 111-10 is greater than 1 ..‘'cnr* long Fig. 
111-1 5 IS a scan more than 2cnr* long that is part of a .‘'cnr* mode, the longest one for this 
pariicular laser linfortunately this mode contained very few H 2 S lines. The scan shown here 
was used to transfer a calibration from the II 2 S b;4~^25 

calibrated with a N 2 O line in the lower frequency [XHtion of this mode. The 
isi>topic line was then used to calibrate the blended doublet 7 q 7 - 8 |j^. ^| 7 - 8 q}< at 1101 OoOcin'* 
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Figure III-5. Diode Laser Spectrum of H 2 S and N 2 O Near 1 14'4.7cm 
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Figure III-7. Diode Laser Spectrum of H 2 S and N 2 O Near I J 23 cm"* 
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Figure III-9. Diode Laser Spectrum of H 2 S and N 2 O Near 1 1 29 cm' 
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Figure Ill-I I. Diode Laser Spectrum of H 2 S and CO 7 Near 1095.7cm" 
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Figure 111-14. Diode Laser Spectrum of H 2 S and N 2 O Near 1 1 20.5 cm 
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Figure III-15. Diode Laser Spectrum of H 2 S Near 1099 cm*l 
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Figure IFI-17. Diode Laser Spectrum of H 2 S Near 1229.34cm" 
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Figure IH-1 8. Diode Laser Spectrum of H2S Near 1257cm"* 


in the high frequency portion of this same mode. This was the only time a H2S line was cali- 
brated in an indirect manner and over such a long frequency interval. 

Whenever possible, CO2 lines in the very accurately measured 9,6ftm band were used as 
calibration standards. Some of the spectra using CO2 standards are shown in Figs. III-l 1 
through 111-13. Since the CO2 lines used here are at the far end of the CO2 R branch, 
relatively large optical depths of CO2 were required. The weakest CO2 line used was the 
R(58) transition for which a 52 meter path and 10 torr of COi was required. Laig;e opti- 
cal depths of H2S were sometimes used to bring out weaker lines in the spectrum while 
saturating other lines, as can be seen in Figs. 111-4 and 111-6. 

Another feature of the H2S spectrum is absorption lines due to isotopic H2S. The 
natural abundances of the sulphur 34 and sulphur 33 isotopes are 4.29f and 0.75% re- 
spectively. Isotopic lines are indicated in the spectra shown in Figs. Ill-band 111-9. 

Fig. 111-14 is an example of a spectrum not used for any measurements due to insta- 
bilities in the frequency scale. These frequency non-uniformities are evident in the 
etalon fringes and illustrate the importance of recording the etalon fringes at the same 
time as the absorption signal. Subsequent scans of the same laser mode showed no evi- 
dence of frequency instabilities. 

Some examples of spectra used for line strength measurements are given in Figs. 111-16 
through 111-18. A number of scans at different pathlengths of the same transitions are 
superimposed in Fig. 111-16. Also shown is a scan with the gas pressure and pathlength 
increased enough to saturate one transition for the determination of the 100% absorption 
level. The spectra shown here have been normalized to correct for a slightly sloping back- 
ground, These spectra, as well as those in many of the other figures, show H2S doublets char- 
acterized by 3 to 1 intensity ratios which result from the effects of the nuclear spin 
statistics on the lower state rotational energy level populations. 

Fig. 111-18 is an example of the lowest signal to noise ratio obtained for any of the 
spectra recorded. This is also the most dense spectrum recorded for H2S. This spectrum 
has the only lAK ±11=3 line observed in the i >2 band, the 541 - 5)4 transition. 
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Fig. 111-17 shows the shortest laser mode used in this work. This spectrum and those in 
Fig. 11-16 are representative of the signal-to -noise ratio for the majority ol the strength meas- 
urements. The appearance of two H 2 S lines in the mode shown in Fig. III-17 made the line 
assignments possible. 

III. D. Grating Spectrometer Measurements 

About forty percent of the H 2 S line frequencies were measured with a high resolution 
grating spectrometer at the National Bureau of Standards. This work was performed before 
completing the construction of the diode laser spectrometer. Only a portion of the more 
strongly absorbing R and branches from 1 180cm'^ to 1250cm'^ was recorded. High 
resolution measurements of the weaker P branch region could not be made since only a 
60cm long absorption cell was available. With a continuous spectrum 70 cm"^ long in the 
R branch the first quantum number assignments in the band could be easily made and an 
initial band fit obtained before the diode laser study. The molecular constants from this 
fit were used to generate a calculated spectrum of the 1^2 band of H 2 S. This spectrum en- 
abled the line assignments in the diode laser spectra to be made quickly and easily. The 
grating spectra were also taken since the frequency measurements made with this instru- 
ment are more accurate than the calibration line frequencies available for use with the diode 
laser spectrometer from 1220cm"* to 1235 cm**. 

The spectrometer, which has been described in detail elsewhere^^^), has an Ebert opti- 
cal design with a 3.81 meter focal length and a 40.8 cm wid.* echelle diffraction grating. It 
was not practical to measure any line strengths with the grating spectrometer due to poor 
signal-to-noise, low resolution, and non-linear 100% transmission levels. 

The H 2 S spectra were taken at a resolution of 0.04cm"* and a signal-to-noise ratio of 
about 70 with a 1 second time constant on the output of the lock-in amplifier. The 
60cm long gas cell contained about 20 torr of H 2 S. The spectrum from 1 180cm'* to 
1250cm'* was scanned and recorded on magnetic tape three times. It was calibrated by 
the standard method of inserting well known CO absorption lines appearing in other orders 
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of the grating into the H 2 S spectrum every 4cm'^. A computer peak finder program found 
the line centers and computed a quadratic or cubic polynominal calibration curve from 
which the H 2 S line center frequencies were calculated. The reported frequencies are aver- 
ages of the three trials. The calibration lines could generally be fit to better than ±0,001 
cm** and the overall error in the measured frequencies is believed to be about ±0.002 cm**. 
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CHAPTER IV 


EXPERIMENTAL RESULTS - DETERMINATION OF THE H^S LINE 
FREQUENCIES AND STRENGTHS 


IV. A. Introduction 

The determination of the H 2 S frequencies and strengths from the diode laser spectra 
will be discussed in this chapter. This will include a brief description of the computer 
analysis of the raw digital data. Errors in the line center frequency measurements from 
using solid germanium etalons to determine the relative frequency scale will be analyzed 
and tests of the overall experimental accuracy presented. The method of measuring the 
line strengths from the line center absorption at low pressure and the errors associated 
with this method will be discussed in detail. These errors include the effects of laser line- 
shape and pressure broadening on the line center absorption. Corrections to the measured 
line strengths for either a Lorentzian or Gaussian laser lineshape will be presented. 

A total of about 1200 spectra were recorded with the diode laser in the course of 
this work. A total of 1 26 line frequencies were measured, distributed as 32 P branch. 25 R 
branch, 14 branch, and 55 Q*^ branch lines. Of these lines. 77 were measured with the 
diode laser spectrometer and '*9 with the grating spectrometer. All 49 lines measured with 
the grating spectrometer were in the R and Q*^ branches between 1081 cm'* and 1 250cm**. 
A total of 94 line strengths were measured (all with the diode laser), distributed as 32 P 
branch. 16 R branch, 13 Q***^ branch, and 33 0*^ branch lines. All of the lines measured 
here fell between 1081 cm** and 1260cnr*. 

The H 2 S line assignments were made using a spectrum calculated with the molecular 
constants obtained from a fit of the grating spectrometer data. This spectrum, along with 
the known frequency setting of the diode laser mode selecting grating monochromator and 
with the H 2 S and calibration gas line patterns in the spectra allowed the assignments to be 
made with ease. 
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IV. B. Frequency Mcasaiements and Error Analysi s 


The H2S line center frequencies were determined by measuring the frequency difference 
between a well known reference absorption line and the H2S line. The relative frequency 
scale was derived from the etalon fringe peaks using a calculated free spectral range for the 
etalon. A total of 66 lines was calibrated relative to N2O lines. 7 relative to CO2 lines, and 
4 relative to NH3 lines. The N2O calibration stand .rds^^^) were either from the band 
centered at 1 284.9 cm'^ or from the 2i>2 band centered at 1 168.1 cm"*. The errors quoted 
for the N2O line frequencies range from ±0.00! 8cm"* to ±0.00098 cm"* for the v\ band and 
from ±0.00 18 cm"* to ±0.00035 cm"* for the 2v2 band. The CO2 calibration line frequencies^^®) 
are known to better than 0.00001 cm"*, well within the precision of this experiment. The NH3 
lines have been measured with a high resolution grating spectrometer^^* ) to an accuracy of 
about ±0.005 cm"*. The average uncertainty in the calibration line frequencies used here was 
±.0.00 10 cm**. 

The digitized absorption spectra and etalon fringes were analyzed with a computer pro- 
gram to obtain relative line center frequencies, line center absorptions, linewidths, and etalon 
fringe peak positions. A least squares fit of the \0Q9c transmission region surrounding each 
absorption line to a cubic polynominal was used to determine the 100% transmission level 
across the line. The line center frequency, line center absorption, and linewidth were deter- 
mined from least squares fits of portions of the line profile to a cubic polynominal. 

Cubic polynominals were chosen for the fits to the line profile rather than a single fit to a 
Doppler lineshape for reasons of computing speed (a desktop computer was used and over 
1 100 spectra were analyzed) and because of possible deviations of the true lineshape from 
the Doppler, especially at highei pressures. The etalon peaks were also determined from 
fitting the region around each peak to a cubic polynominal. Since most of the line 
strength spectra were only 0.1 cm"* to 0.2 enr* long, the line profile fits usually included 
30 to 60 points. For the longer frequency calibration scans the number of points per line 
fit ranged from 6 to 12. 
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The fringe peak positions were fit to cubic spline functions for interpolation between 
fringe peaks. In most cases the fringe separations were constant enough that linear interpola- 
tion between peaks would have sufficed. The H 2 S line frequencies and widths were then 
converted to wavenumber units using the free spectral range of the etalon and the calibra- 
tion line frequency. 

It is of utmost importance to use an accurate value for the etalon free spectral range. 
For example, a 0.5 cm"* separation between the calibration line and unknown line trans- 
lates into about 30 etalon fringes. If we wish to measure this frequency difference to 
0.0002 cm"* the etalon free spectral range must be known to 6.7 x 10"^ cm"* or 0.04%. 

The free spectral range of an etalon is given to a high degree of accuracy by 


At'fsR - 


1 


/ dn 
2t n + — 
\ dX 




(lV-1) 


where n is the index of refraction of germanium at the wavelength X and t is the etalon 
thickness. To calculate Aj^fSR accurately very good values for n(X) and t are required. A 
careful measurement of the etalon length gave t = 7.65048 ± 0.00038 cm, accurate 
enough so to not introduce large errors into the calculated free spectral range. Fortunately, 
the index of refraction of germanium has been measured at eleven wavelengths between 
2.55pm and 12.36ium to an accuracy of 0 . 01 5 %.^'*^) a cubic spline fit to this data was 
used to interpolate between the measured values of the index of refraction and to deter- 
mine dn/dX. The cubic spline coefficients for this fit are given in Table IV-1. Before cal- 
culating the free spectral range, the room temperature values for the index of refraction 
obtained using Table IV-1 were corrected to T = 303.3 K, the temperature at which the 
etalon was stabilized. 

There are several other errors associated with the use of etalons in relative frequency 
calibration. One of these, fringe pattern fluctuations from etalon temperature drifts, was 
shown in Chapter III to have a maximum value of ±0.0001 cm"*. Another possible cali- 
bration error comes from angular misalignment of the etalon relative to the collimated 
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Table IV-1 


Spline Fit Coefficients for Index of Refraction of Germanium 


(a) Wavelength, Xj 

(Min) 

^**)lndex, nj 

B 


C 


D 

2.554 

4.06230 

-2.6322 x 

10-2 

-3.3061 X 

10-* 

1.0568 

2.652 

4.05''54 

-6.0672 X 

10-2 

-1.9902 X 

10-2 

1.0568 

2.732 

4.05310 

-4.3565 X 

10-2 

2.3374 X 

10-* 

-9.5556 X 10-1 

2.856 

4.04947 

-2.9676 X 

10-2 

-1.2173 X 

10-* 

7.2885 X 10-* 

2.958 

4.04595 

-3.1760 X 

10-2 

1.0130 X 

10-* 

-2.6202 X 10-1 

3.090 

4.04292 

-1.8714 X 

10-2 

-2.4617 X 

10-3 

3.2369 X 10-3 

4.120 

4.02457 

-1.3483 X 

10-2 

7.5404 X 

10-3 

-2.1275 X 10-3 

5.190 

4.01617 

-4.6538 X 

10-3 

7.1115 X 

10^ 

-4.1458 X 10-5 

8.230 

4.00743 

-1.4794 X 

10-3 

3.3305 X 

10^ 

-1.0377 X 10-5 

10.270 

4.00571 

-2.5008 X 

10*4 

2.6955 X 

10^ 

-1.0377 X 10-5 

12.360 

4.00627 







and nj are experimental values from ref. (28) for T = 297 K. 
(b)n(X) = nj + B(X - Xj) + C(X - Xj)2 + D(X - Xj)3, Xj < X < Xj+j . 


laser beam. Since the mode selecting grating monochromator is before the etalon in this 
experiment, the laser beam angle will shift as the laser frequency is scanned (the grating 
angle is kept fixed during the scan). This changing angular misalignment will produce an 
observed free spectral range different from the one calculated using equation lV-1. 

Flicker et iiave shown that this can result in large calibration errors when an air- 

gap etalon is used in conjunction with a high dispersion monochromator and a short focal 
length collimati.ig element before the etalon. This error, defined to be the actual separa- 
tion between two lines minus the observed frequency separation, was calculated using our 
experimental parameters and was found to equal +0.00006 cm'* for a 0.8 cm** measure- 
ment. It should be noted that the large index of refraction of germanium permits rela- 
tively large angular mis-alignments of the etalon to be made without introducing 
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significant changes in the observed free spectral range. The etalon errors are summarized 
in Table IV-2. The total error due to the etalon in the relative frequency measurements 
should range from approximately +0.0002 cm‘^ to -0.0001 cm’*. 


Table IV-2 
Etalon Errors 


Error Source 

•Error (cm’* ) 

Index of Refraction 

±0.(X)012 

Etalon Thickness 

±0.00004 

Scanning Error 

+0.00006 

Temperature Drift 

±0.00010 


*Assuming a 0.8cm‘* measurement, except for the Temperature Drift, which is independent 
of the frequency range measured. 


Fortunately, a total of seven laser modes were recorded that contained two calibrr- 
tion standards, separated by about 0.8 cm**, which could be used to test the accuracy of 
the frequency difference measurements. The frequency separation of the calibration 
standards was always very well known compared to the precision of our measurements. 

The results of these tests are given in Table iV-3. The measured frequency separations 
are averages of three to five trials and the assigned error is the standard deviation of 
these trials. 

The average deviation between the actual frequency difference and the measured dif- 
ference was -0.00004cm’*, well below the precision of any one measurement. This indi- 
cates that there are no significant systematic errors in the calculated values for the etalon 
free spectral range or in the data. Also, the average assigned expierimental error of 
0.00047 cm’* is close to the standard deviation of the measured frequency separations from 
the correct values. Because of this agreement, and since any systematic errors appear to be 
very small, the experimental error in the line frequencies was taken to be the standard 
deviation of the trials used in the measurement of each line. This error is probably domi- 
nated by random laser frequency fluctuations and etalon temperature drifts. 
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Table lV-3 


Diode Laser Calibration Accuracy Tests 


Lines Used 

Approx. Frequency 
Separation - A 
(cm-^) 

Calc A - Meas A 
(cm’l) 

Assigned 

Experimental Error 
in A (cm‘^) 

CO 2 - CO 2 

1.00 

-0.00078 

0.00037 

N 2 O - N 2 O 

0.77 

+0.00049 

0.00033 

N 2 O - N 2 O 

0.77 

-0.00036 

0.00072 

N 2 O - N 2 O 

0.77 

+0.00001 

0.00015 

N 2 O - N 2 O 

0.80 

+0.00049 

0.00064 

NjO - N 2 O 

0.80 

-0.00018 

0.00057 

N 2 O - N 2 O 

0.80 

+0.00004 

0.00050 



Average = -0.00004 cm'^ 

0.00047 



*(±0.00042 cm-1) 



*This is the standard deviation of the measured frequency separations from the true 
separations. 


The measured line frequencies are given in Table lV-4. The calculated error given in this 
table is the root mean square of the measurement error discussed above and the uncertainty in 
the calibration line frequency. The average calibration line error was ±0.0010cm‘* and the av- 
erage total error for the diode laser frequency measurements was ±0.001 3 cm'^ The lines de- 
noted with an asterisk in Table IV-4 were measured with the NBS grating spectrometer. The 
average total error for all the frequencies given in Table IV-4 is ±0.001 6 cm'T 

IV. C. Strength Measurements and Error Analysis 

Only a small fraction of the studies made with diode lasers have been concerned with the 
measurement of line strengths. Usually, the method of equivalent widths('^>'^5) has been 
used to determine the line strength from the observed absorption.('^) The line strength can also 
be determined from the line center absorption if the laser linewidth is negligible compared to 
the width of the absorption line.(^^) Or, if the laser linewidth is not too large, the line 


63 







Table IV -4 


Measured Transition Frequencies in »»2 Band of H2S 


Upper Lower 

J'(Ka'. K^.'H(Ka. K,.) 

Observed 
Frequency 
(cnr* ) 

Experimental 
Error 
(cm** ) 

(a)Obs-Calc 
Frequency 
(cm** ) 

6 < 2 , 5 )- 5 ( l , 4 > 

1257.0689 

0.0040 

- 0.0004 

6 ( l , S )- 5 ( 2 > 4 > 

1257.0100 

0.0034 

- 0.0016 

S < 4 , i )- S < l , 4 > 

1256.9713 

0.0015 

- 0.0005 

5 < 2 , 3 )- 4 ( 3 , 2 > 

1256.8985 

0.0025 

- 0.0017 

?< 4 , 6 >- 9 < 3 , 7 ) 

1256.7848 

0.0014 

- 0.0014 

5 < 3 , 2 )- 4 < 4 , l ) 

1256.7639 

0.0016 

- 0.0007 

9 ( 3 , 6 )- 9 ( 2 , 7 > 

1256.7313 

0.0014 

- 0.0012 

4 < 3 , 2 )- 3 ( 2 , i > 

1253.4679 

0.0010 

0.0000 

7 < i , 6 )“ 7 < 0 , 7 ) 

* 1249.2231 

0.0020 

0.0008 

3 < 3 , 0 )- 2 < 2 > i > 

* 1248.3825 

0.0020 

0.0018 

5 < 2 , 4 )- 4 < l , 3 ) 

* 1247.5366 

0.0020 

- 0.0006 

5 < l , 4 )- 4 < 2 , 3 ) 

* 1247.1994 

0.0020 

- 0.0001 

8 ( 4 , 5 )“ 8 < 3 , 6 > 

1246.0787 

0.0012 

- 0.0010 

8 < 3 , S )' 8 ( 2 , 6 ) 

1245.8463 

0.0011 

- 0.0010 

6 ( 1 , 6 )- S ( 0 , 5 > 

1243.6170 

0.0013 

0.0001 

9 < 6 , 4 >- 9 < S , 5 > 

1243.5780 

0.0011 

0.0010 

3 ( 3 , l )- 2 < 2 , 0 > 

* 1243.0679 

0.0020 

0.0010 

7 ( 3 , S )- 7 ( 2 , 6 ) 

1242.5596 

0.0011 

- 0.0013 

7 < 2 , 5 )- 7 < i , 6 ) 

1242.5022 

0.0013 

- 0.0012 

4 ( 2 , 2 >- 3 < 3 , l ) 

* 1241.3277 

0.0020 

0.0070 

6 < 2 > 5 )- 6 < i ^ 6 ) 

• 1239.3377 

0.0020 

- 0.0038 

8 ( S , 4 )- 8 ( 4 , 5 ) 

* 1239.1324 

0.0020 

- 0 . 0008 

4 < 2 , 3 )- 3 < i , 2 ) 

* 1238.3986 

0.0020 

- 0.0008 

9 ( S , 4 )- 9 < 4 , S ) 

* 1238.2410 

0.0020 

0.0021 

4 ( l , 3 )- 3 ( 2 , 2 ) 

1236.7275 

0.0015 

0.0004 

8 ( 4 , 4 )- e < 3 , S > 

1236.6533 

0.0030 

- 0 . 00 Q 3 

8 < 7 , 2 )- 8 < 6 , 3 ) 

* 1235.7061 

0.0020 

0.0004 

7 ( 4 , 4 )- 7 ( 3 , S ) 

* 1235.3761 

0.0020 

- 0.0022 

8 < 6 > 3 )- 8 < 5 , 4 ) 

1234.6136 

0.0019 

- 0.0004 

5 ( l , S >- 4 < 0 , 4 ) 

1234.5850 

0.0019 

0.0003 

5 < 0 , 5 )- 4 < i , 4 ) 

1234.5788 

0.0019 

0.0004 

6 ( 3 , 4 )- 6 < 2 > 5 > 

* 1232.1437 

0.0020 

- 0 . 0020 

7 ( 7 , 0 )- 7 < 6 , l ) 

* 1230.3327 

0.0020 

- 0 . 0003 

3 < 2 , 2 )- 2 < i , l ) 

* 1229.8479 

0.0020 

0.0002 

5 ( 2 , 4 >- 5 < i , S ) 

* 1229.3386 

0.0020 

0.0049 

5 < l , 4 )- 5 ( 0 , S ) 

* 1229.2731 

0.0020 

- 0.0013 

7 < 6 , 2 >- 7 < S , 3 ) 

* 1228.1318 

0.0020 

- 0.0034 

2 < 2 , 0 >- l ( i , i ) 

* 1227.5982 

0.0020 

0.0028 

9 < 7 , 2 )~ 9 < 6 , 3 ) 

* 1226.3977 

0.0020 

- 0.0021 

9 < 6 , 3 >- 9 < S , 4 ) 

* 1225.9818 

0.0020 

- 0.0023 

4 < i , 4 )- 3 ( 0 , 3 ) 

* 1225.4417 

0.0020 

- 0.0015 

6 ( 4 , 3 >' 6 ( 3 , 4 > 

* 1225.1607 

0.0020 

- 0 . 0012 

6 < 6 , i >- 6 ( 5 , 2 > 

1224.3336 

0.0048 

- 0.0005 

3 ( l , 2 )- 2 < 2 , l ) 

* 1224.2811 

0.0020 

0.0019 

8 ( 5 , 3 )- B < 4 , 4 > 

* 1223.9543 

0.0020 

0.1006 
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Table IV- 4 (Continued) 


Upper Lower 

J'(Ka'. Kc')-J(Ka. 

Observed 

Frequency 

(cm’*) 

Experimental 
Error 
(cm** ) 

(*)Obs-Calc 
Frequency 
(cm** ) 

7<4,3)-7<3>4) 

*1222.9320 

0.0020 

>0.0022 

6(3,3>-6(2,4> 

*1221.9920 

0.0020 

-0.0037 

S(3,3)>5<2,4> 

*1221.8686 

0.0020 

0.0044 

4<S,2)-6<4,3> 

*1221.6799 

0.0020 

0.0035 

S<2,3)-S(l>4) 

1220.6479 

0.0020 

-0.0023 

2<2,l)-i(l,0) 

*2219.9501 

0.0020 

0.0015 

8<6»2>-8(S>3) 

*1218.6091 

0.0020 

0.0050 

7<6,l)-7<5,2> 

*1217.9880 

0.0020 

0.0011 

5(5,l)-S(4,2) 

*1216.8640 

0.0020 

0.0029 

5<4,2)-5<3,3) 

1216.4365 

0.0013 

-0.0005 

3<i,3>-2(0>2> 

1216.2968 

0.0010 

-0.0004 

3(0,3)>2(1,2> 

*1215.9665 

0.0020 

0.0007 

7(S,2>-7<4,3> 

*1213.3871 

0.0020 

-0.0006 

4(3,2)-4<2,3) 

*1212.4180 

0.0020 

0.0002 

5<5,0>-S(4,l> 

*1211.1931 

0.0020 

0.0026 

6(4,2>-6(3,3> 

*1210.4286 

0.0020 

-0.0005 

4<4,l)-4<3,2> 

*1210.3639 

0.0020 

0.0013 

6(5,i)-6(4,2> 

*1209.5985 

0.0020 

-0.0004 

3(2,2)~3(i,3) 

*1209.5252 

0.0020 

0.0026 

5<3,2)-5<2,3> 

1209.0831 

0.0006 

0.0003 

4<2,2)-4(l,3) 

1208.5210 

0.0005 

0.0002 

3<i,2)-3<0,3> 

1208.0164 

0.0005 

-0.0002 

2(l,2)-i<0,i) 

*1207.4508 

0.0020 

0.0012 

2<0,2>-l(i,l> 

*1205.7807 

0.0020 

0.0026 

3<3,i)-3(2,2> 

*1205.0435 

0.0020 

0.0003 

S<4,l)-5(3,2> 

*1203.4497 

0.0020 

-0.0004 

4(4,0)-4(3,i> 

1202.8244 

0.0005 

-0.0006 

2<2,i>-2(l,2) 

*1201.0267 

0.0020 

0.0002 

4(3,i)-4<2>2> 

1199.5292 

0.0005 

-0.0003 

i(l>i)-0(0,0> 

*1197.9709 

0.0020 

0.0037 

3<2,i)-3<l,2) 

*1197.4967 

0.0020 

0.0007 

2<l,l)-2(0,2> 

1196.8217 

0.0004 

0.0001 

3<3,0)-3<2,i> 

*1196.1755 

0.0020 

0.0014 

2(2,0)-2<l,i> 

1191.5453 

0.0006 

-0.0001 

l<l,0)-i(0,l) 

1188.7735 

0.0004 

-0.0002 

6<3,3)-6(4,2) 

1171.4822 

0.0007 

-0.0006 

I<0,i>-2<1,2) 

1158.1694 

0.0005 

-0.0001 

2(i ,2)-3<«,3) 

1149.7713 

0.0005 

-0.0003 

2<0,2)-3<i,3) 

1149.4026 

0.0008 

-0.0004 

4<0,4)-4(i,3) 

1148.7198 

0.0004 

-0.0003 

4<l,4>-4(2,3> 

1148.4491 

0.0004 

-0.0001 

Ki,0)-2<2,i> 

1147.3581 

0.0006 

-0.0003 

2<2,l)-3(i,2> 

1144.2681 

0.0008 

0.0001 

3<l,3)-4<0,4) 

1140.1403 

0.0006 

-0.0007 

3<0,3)-4(l,4) 

1140.0855 

0.0008 

-0.0004 

Kl,i>-2<2,0) 

1139.5982 

0.0005 

-0.0001 

2<l,l)-3<2,2) 

1138.4452 

0.0004 

-0.0001 

3(2.2)-4<i,3> 

1132.8476 

0.0004 

0.0004 

7<i,6>-7<2,5> 

1132.7146 

0.0007 

-0.0000 
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Table IV-4 (Continued) 


Upper Lower 

J'(Ka', Kc')-J(Ka, K^) 

Observed 
Frequency 
(cm"' ) 

Experimental 
Error 
(cm-* ) 

(a)Obs-Calc 
Frequency 
(cm** ) 

7<2,6)-7<3,S) 

1132.6721 

0.0007 

0.0005 

4<1,4>-S(0,S> 

1130.5247 

0.0005 

0.0006 

4(0,4)-5(i,S) 

1130.5172 

0.0005 

0.0007 

S(5,0)-6(4 3) 

1129.1619 

0.0003 

0.0000 

9(2,7>-9(3>A) 

1128.5689 

0.0011 

0.0007 

9(3,7)-9(4,6> 

1128.5302 

0.0013 

0.0004 

6(G^6>-6(i,5> 

1128.3660 

0.0012 

0.0002 

6<l,6)-6(2,5) 

1128.3583 

0.0013 

0.0001 

2(2,0)-3(3,l) 

1 127 . 3454 

0.0007 

0.0005 

4(2,3)-5(i,4) 

1123.2380 

0.0010 

-0 . 0005 

8(i,7)-8(2,6) 

1123.0877 

0.0013 

-0.0007 

8<2,7)-8(3,6> 

1123.0794 

0.0014 

-0 . 0008 

4<i,3)-5(2,4) 

1122.8542 

0.0011 

-0.0007 

S<0,S>-6<1,6) 

1120.8093 

0.0016 

'0.0016 

4<4,i>-5(3,2> 

1120.5884 

0.0016 

-0.0025 

10(2,8)-10(3,7> 

1119.4978 

0.0018 

0 . 0002 

iOt3,8)-iO(4,7) 

1119.4895 

0.0019 

0 . 0003 

6<l,6)-7(0,7) 

1110.9883 

0.0050 

-0.0022 

6<2,5>-7(l,6> 

1104.9311 

0.0050 

-0.0039 

6(l,S)-7(2>6) 

1104.9201 

0.0050 

-0.0035 

5<4,2)-6(3,3> 

1104.3112 

0.0050 

-0.0032 

7(0,7)-8(l,B> 

1101.0598 

0.0010 

-0 . 0001 

7(l,7)-8<0,8) 

1101.0598 

0.0010 

-0.0001 

6(3,4>-7<2,5> 

1098.2698 

0.0004 

-0.0006 

6<2,4)-7(3,5> 

1097.9380 

0.0005 

-0.0003 

7(l,6)-8(2,7> 

1095.6973 

0.0003 

0.0000 

8(l,8)-9(0,9> 

1091.0226 

0.0001 

0 . 0000 

8(0,8>-9(l,9> 

1091.0226 

0.0001 

0.0000 

7<3>5)-8(2,6> 

1089.6241 

0.0002 

0.0002 

7<2,S>-8(3,6> 

1089.5570 

0.0004 

0.0003 

8<2,7>-9<l,8> 

1086.3629 

0 . 0002 

-0.0000 

7(3,4)-8(4,S) 

1082.0329 

0.0003 

0.0002 


♦Indicates that the transition frequency was measured with the NBS grating spectrometer. 
(a)See Ch. V for details on the calculated frequencies. 
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strength can be determined from the line center absorption and then corrected for the distor- 
tion of the absorption line profile by the laser. For the reasons outlined below, this second 
method (called the direct method) is expected to yield more accurate line strengths and was 
used in this work. 

One of the most difficult problems associated with any line strength measurement is 
the determination of the 100% transmission level in the region of the line. This is especially 
difficult for high pressure absorption when the line wing, described by a Lorentzian lineshape, 
absorbs an appreciable distance away from the line center. To avoid this problem the gas 
pressure was kept as low as possible in this experiment, usually between 0.15 and 0.25 torr. 
Under these conditions, the lineshape is very nearly Doppler, so that the absorption in the 
line wing goes to zero very rapidly. 

In the method of equivalent widths, the line strength is determined from the int^rated 
absorption of the line, and is therefore independent of the laser lineshape. This method has 
been used for most diode laser line strength measurements since the observed low pressure 
linewidths in these studies were 8% to 40% wider than the Doppler width, indicating that the 
observed absorption line profiles were significantly distorted. The accuracy of the equivalent 
width method, though, is more sensitive than the direct method to errors in the 100% trans- 
mission baseline^'*^) and to slight frequency non-linearities in the laser tuning.^^^) 

In our case, the observed low pressure linewidths were only about 1% larger than the 
Doppler width. Therefore, it was possible to measure the line strengths using the direct 
method, with corrections for distortion of the absorption line by the laser, and to avoid some 
of the errors associated with the use of equivalent widths. Also, the H2S data contained a num- 
ber of blended lines with 3 to 1 strength ratios. The strengths cf these blended lines can be 
measured much more easily using the direct method than the method of equivalent widths. 

The intensity, I(i^), of a spectrum containing a single absorption line is given by Beer’s 

law, 

I(i') = IqC*') (lV-2) 
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where Io(v) is the intensity if there is no absorption, k(i^) is the absorption coefficient, 
and X is the optical depth, which is equal to the pressure-pathlength product. For a 
Doppler linesliope the absorption coefficient has the form 


S 

k(p) = — 

'J'd 



(IV-3) 


where S is the line strength, Vq is the frequency of the line center, and 7 q is the Doppler 
width and is equal to the half widtii at half height of the absorption coefficient. The 
Doppler width is given by 


= »'o 


2kT ln2 
Mc2 


(IV^) 


where M is the molecular mass, k is Boltzmann’s constant, T is the temperature, and c is 
the speed of light. From the line center transmission, T^, = 
the line strength since 


Sd - 




PL 


In (1/To) 



(IV-5) 


The line strength calculated in this manner will be denoted by Sj) since a Doppler line- 
shape was assumed. This quantity was determined for each absorption line and was sub- 
sequently corrected for the effects of pressure broadening and the laser lineshape on the 
absorption line profile. 

If a saturated absorption scan was available, the 100% absorption level in each spec- 
trum was compared to the signal level recorded when the shutter was blocking the laser 
beam. If the difference between the snutter signal and the 100% absorption level was 
greater than a few tenths of one percent, the saturated absorption level was used instead 
of the shutter zero as the IgCVo) = 0 level. Saturated absorption scans were made for all 
but seventeen of the lines measured. These seventeen lines are marked in Table lV-5. 
which is a listing of the measured line strengths. 
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Fifteen of these seventeen lines were in the P branch region since the H 2 S absorption 
is so much weaker on the low frequency side of the band. The saturated absorption 
levels of all but two of the other P branch region lines were less than 1% higher than the 
shutter zero, and averaged only 0.3% higher. The other two saturated absorption levels 
were 1.9% and 1.6% higher than the shutter zero. These two lines were recorded with the 
highest currents used with the P branch r^ion laser, for which poorer mode quality and 
separation can be expected. Since the lines for which saturated absorption scans were not 
made occurred at lower laser currents, it is likely that the shutter zero is an accurate meas- 
ure of the 100% absorption for these scans. Rather high error bars were assigned to most 
of the strengths measured for these lines anyway since their spectra had to be recorded at 
relatively high pressures. 

Only two lines recorded with the R branch region laser could not be saturated to 
check the shutter zero level. The shutter zero and the saturated absorption level were 
within 3% of each other for all but four of the other transitions in this region. The satu- 
rated absorption level was from 3% to 9% higher than the shutter zero for these four transi- 
tions. The two lines that could not be saturated in the R branch region were recorded 
with low laser injection currents for which the agreement between the shutter zero and 
the 100% absorption level was always better than 1%, when it could be measured. 

Corrections to Sp for the effects of any pressure broadening can be determined using 
the Voigt lineshape, which is a convolution of the Doppler and pressure broadened 
Lorentzian absorption line profiles. For 80% of the lines measured here the pressure was 
less than 0.4 torr, so that these corrections are small, although not negligible when com- 
pared to other error sources. The corrections for these low pressure lines were determined 
in conjunction with corrections for distortion of the line profile by the laser and will be 
discussed later. Eighteen lines were recorded with pressures ranging from 0.4 to 4 torr. 

The pressure broadening corrections to these lines will now be discussed. 
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The Voigt lineshape is given 




y2 + (X - t)2 


(IV-6) 


where 



and 

7 l = TtP. 

kg is the absorption coefficient at Vq, the line center, for a Doppler lineshape. 7 l is the half 
width at half height of the pressure broadened Lorentzian lineshape, 7^ is the pressure broad- 
ening coefficient in cm'Vatm, and P is the pressure in atm. At the line center k''(v = Vq) can 
be considerably simplified(‘^8) to give 

= kocV^l -erf(y)). (IV-7) 

The line strength corrected for the effects of pressure broadening, denoted by Sp, is there- 
fore given by 

Sd 

S = . (IV-8) 

ey^{l -erf(y)) 


For the pressures used here y is always less than 1 so that erf(y) can be accurately calcu- 
lated from the power series 


erf(y) = - ^ f e*^^ dt = 

Jck 


2 


(IV-9) 


2 ^ (-l)ny2nti 

^ ho n!(2n+ !)■ 

The correction to Sd for pressure broadening as a function of y is plotted in Fig. lV-1 . 
For example, if a line with a Doppler width of 0.001 2 cm** and a pressure broadening coeffi- 
cient of 0.1 cm**/atm was recorded at a pressure of 1 torr. Fig. IV- 1 shows that the value of 
Sq found for this line must be increased 109f to account for the decrease in the line center 
absorption due to pressure broadening. 

If 7 l is known, the correction to Sq for pressure broadening can be calculated. 
Olivero and Longbothum^^^) have developed an empirical equation for the half width of the 
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Voigt lineshape, 7y, as a function of and 7£). This equation can be invert ed^**^) to 
obtain 

\ = T'V ^7.7254 ~ 6.7254 \/(l +0.3195 (7 d/7v> 2) ^ (lV-10) 

This relation was used to determine 7 l and therefore 7^ from the measured halfwid^hs of 
nine of the lines, all recorded at pressures greater than 1 torr. The residual line broadening 
of about 1% due to the diode laser linewidth was neglected since the observed linewidths 
were 20% to 50% wider than the Doppler width. 

For these nine lines, the measured pressure broadening coefficient was used for the 
correction of each line strength. The error in the 7^ values was estimated to be about 
± 0.02 cm"* /atm. The uncertainty in the line strength correction from the uncertainty of 
±0.02cm"Vatm in 7^ was included in the total error for these line strength measurements. 
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The measured values of 7^ ranged from 0.1 1 to 0.l9cm"*/atm, the average being 0.15 
cm’* /atm. The lower state J values of these lines varied from J = 2 to J = 10, 

The average value for of 0.1 5 cm’* /atm was used to determine the pressure broad- 
ening corrections for the lemainder of lines measured at pressures between 0.4 and 1 torr. 
Errors corresponding to an uncertainty of ±0.05 cm’* /atm in 7^ were included in the error 
bars for these line strengths. This uncertainty of ±0.05 cm" * in the average value of 7^ is large 
enough to encompass all nine of the values of 7^ measured here. 

The corrections to Sp for these 18 lines ranged from 7% to 60% and the resulting 
errors in the line strengths from the uncertainty in these corrections ranged from 2% to 
10 %. 

A good estimate of the diode laser lineshape is required to determine its effect on an 
absorption line profile. Studic^^®>^* * of line profile distortions in diode laser spectra have 
assumed that the laser lineshape is Gaussian. Also, previous line strength measuremci.ts(^^) 
using the direct method were not corrected for laser lineshape effects since the assumption 
of a Gaussian shape led to a negligible correction. 

In a review of the literature, six direct measurements of the lineshape and linewidth 
of a diode laser were found. These measurements were all performed hy heterodyn- 
ing the output of a diode laser with a highly stable and very narrow CO 2 or CO laser line. 

In all six cases the observed diode laser lineshape was clearly best described by a Lorentzian 
lineshape. The diode laser linewidths in these measurements varied from less than 2 MHz to 
greater than 20MHz full width at half height. We therefore determined our corrections to 
the measured line strengths using a Lorentzian lineshape for the laser emission p 'ofile. 

To determine the effect of the laser lineshape on the line prcTile and to find correc- 
tions for the measured line strengths a procedure similar to that used in grating spectroscopy 
was followed The laser lineshape was convolved with a scries of Doppler rbsorption 
line profiles for various ratios of the laser linewidth to the Doppler width. Using these 
convolutions, corrections to the line strength can be determined either from the laser line- 
width, it is known, or from the observed broadening of the line above the Doppler 
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width, if this broadening can be measured accurately enough. Convolutions were per- 
formed for both a Gaussian and Lorentzian lineshape. The Gaussian lineshape convolutions 
were made for comparison purposes and to check our convolution routine since the results 
of other convolutions of a Gaussian with a Doppler line profile are available in the 

literature.(^2) 

The convolved line profile, for a laser lineshape described by the function f(v - v') is 
given by 

“ f ^TRUE(^ ) (IV-11) 

-oo 

where Aqbs and AjpuE ^^e observed and true absorption of the line. For a Lorentzian 
laser lineshape 


Kp - p') 

and for a Gaussian laser lineshape 


J ODL 

n (p - p')^ + 


(lV-12) 


f(p - p') 

where apL is the laser linewidth. 
so that 



(IV- 13) 


«DL V JT 

The lineshape functions have beer propeiiy normalized 



f(i' - i^') dt'' = 1 . 


(IV- 14) 


The convolutions were performed numciically with a desktop computer. The ratio of 
the laser linewidth to the Doppler width was varied from 2% to 50% and the true absorp- 
tion of the Doppler profile was varied from about 3% to 97%. To save computation time 
the convolved profile ' as only evaluated at line center and in the vicinity of its halfwidth 
point, defined as the halfwidth of the effective absorption coefficient giving rise to this 
convolved profile. 

The results of these convolutions, in terms of corrections to the measured line 
strengths, are plotted in Fig. IV-2 for a Lorentzian laser lineshape and in Fig. lV-3 for a 
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Gaussian lineshape. In each of these figures, lines of constant observed absorption are 
plotted as well as correction curves for a constant ratio of the laser width to the Doppler 
width. 

As one would expect, for a given laser linewidth the line center absorption is decreased 
much more by the Lorentzian lineshape since its wings extend so far compared to the wings of 
the Gaussian lineshape. For example, if the laser width is 20% of the Doppler width, the line 
strength correction for a Guassian laser lineshape is 2% while the correction for a Lorentzian 
laser lineshape is about 23%. If the line strength cxirrection is determined from the observed 
increase in the linewidth the difference between the two laser lineshapes is smaller. If the laser 
lineshape is Gaussian and the absorption line is broadened by X%. the line strength correction 
will also be about X%. If the laser lineshape is Lorentzian though, a broadening of X% in the 
observed line indicates that ? line strength correction of about 2X% is required. 

Since the diode laser linevidth could not be directly measured, the strength correc- 
lions were dete.^mined from the increase of the observed H 2 S linewidths over the Doppler 
width. Since the Doppler vhdth of H 2 S in the band is only about 0.0012cm"*, the 
halfw'idth meas:i-*einent of a single line is not nearly accurate enough to use for the deter- 
mination of the strength correction. Also, even at pressures below 0.4 torr, pressure 
broadening can distort ine line profile as muta a-; the l.i.^er lineshape. Therefore, the aver- 
age of a large number of low pressure lin«.'.vid£'i a;?asuremcms was used to determine the 
combined line strength correction due to pr-'ssure broadening and the laser iineshape. 

A to* 0 low pressure linewidths from 66 different transitions were used to 

measure the average broadening. Tlie .average pressure for these linewidth measafernents 
was 0.24 torr and the average absorption was 45%. Care w.as taken to exclude any 
blended lines from the average. Several transitions recorded using laser modes with poor 
frequency stability were also excluded. The mean ratio of the observed halfwidth to the 
Doppler width foi these transitions was 1.013. The standard deviation of the average line- 
width ratio for ei-ch of the 66 transitions about the mean ratio was ±0.021 . This number is 
believed to be a conservative estimate of »he amount of variation in the line broadening 
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of these transitions since it may also contain a contribution from the rather large random 
error in any one individual linewidth measurement. The probability that this estimate is 
on the high side is supported by the fact that the mean ratio minus one standard deviation 
is less than one, implying an observed width smaller than the Doppler width. The test re- 
sults given in Table IV-3 show that there is no possibility of a systematic error in the line- 
width measurements due to a poor relative frequency scale. The average linewidth ratio 
obtained here should statistically be very accurate, since it was determined from S13 line- 
width measurements. 

The average linewidth ratio plus one standard deviation, or 1 .034, was taken as the 
upper limit of the observed line broadening. If due totally to the laser, this excess broad- 
ening of 3.4% implies a L “r linewidth of 2 MHz. This maximum ratio of the observed 
halfwidth to the Doppler width implies a maximum correction to the measured line strengths 
of 7%, for any combination of broadening due to the laser lineshape and the gas pressure. 
The minimum possible correction to the line strengths was determined to be 3%. This 
minimum is the combined correction for a laser linewidth of 0.5 MHz, the smallest 
linewidth usually observed for diode lasers, and for a pressure broadening coefficient of 
0.07cm*Vatm, which is one half the value of the average pressure broadening coefficient 
measured for H2S. Our best estimate for the low pressure line strength corrections is therefore 
5% ± 2%. If the laser was assumed to have a Gaussian rather than Lorentzian lineshape, 
the line strength coirection was found to be 3% ± 2%, 2% lower than the correction for a 
Lorentzian lineshape. 

As previously mentioned, the strength corrections were determined from a set of ab- 
sorption lines whose average peak absorption was 45%. If the observed excess line broad- 
ening was due solely to the diode laser linewidth. the correction to the line strength would 
depend upon the peak absorption of the line being measuted. In this case, a laser line- 
width that results in a 5% correction at an absorption of 45% should be used for the calcu- 
lation of the strength corrections. If the excess broadening was due completely to pressure 
broadening, the line strength correction is not a function of the line center absorption. 


77 



However, since the correction curves shown in Fig. lV-2 for different laser linewidths are 
almost parallel as the peak absorption is varied, it is not necessary to know the relative 
amounts of laser and pressure broadening very accurately. Also, because the vast majority 
of the strength measurements were made using lines with peak absorptions between 30% 
and 70%, the variation with absorption of any one of the possible laser correction curves 
is at most only ±0.6%. We could therefore use a constant correction factor of 5% for the 
line strengths. 

Eight of the transitions measured here were blended doublets whose lines have a 
knowii strength ratio of 3 to 1. The frequency separation of the two lines in each doublet 
can be calculated from the V 2 band constants, which are discussed in Ch. V. With a knowl- 
edge of the frequency separation and the relative strengths, the strength of the stronger 
line, S|, can be determined from the relation 

Td In O/Tp) 

1 

PL 



where Tg is the minimum transmission of the blended doublet. The factor a is given by 

-ln2 

\ / + . 

3 


a = e 


emission oi ine oienuea aouoiei. i 

\ / + — e V ■>'d / 


(IV-16) 


where u is the frequency at which the blended doublet profile has its peak absorption, and 
and the unblended line center frequencies of the stronger and weaker tines in the 

doublet. The frequency differences v - and v - vj were determined from a theoretical 
calculation of the blended line profiles using a Doppler lineshape and the known strength 
ratio and unblended frequency separation. Alt but one of the blended doublets were meas- 
ured at low pressures. The doublet 122 10 “ ^^3,9* *“3,10 ” *-4,9 measured at about 
?..Z torr, but the error in using a Doppler instead of a Voigt lineshape in this calculation 
should be small compared to other errors in the measurement of its line strength. Sj for 
each doublet was then increased by 5% to correct for the effects of the laser linewidth and 
pressure broadening. Two of the doublets were measured at line separations of less than 
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0.00001 cnr*, while the other separations varied from 0.0001 Scm*^ to 0.00196cm‘*. Since 
the line blending was very severe for all of these doublets, a always had a value close to 4/3. 

The line strength measurements will now be summarized and the use of the direct 
method in making these measurements discussed. The first step was to determine, if possi- 
ble, the 100% absorption level fr^m saturated absorption scans. Then the quantity S[>, 
given in equation IV-5, was calculated for all the unblended lines with the assumption that 
the line profile was Doppler. For all the lines measured at pressures less than 0.4 torr 
was increased b> 5% to correct for the effects of pressure broadening and the laser line- 
shape on the observed line profile. Three sources of error were considered in determining 
the total error for these line strengths. The total error was taken to be the root mean 
square of the 2% error in the combined correction for pressure broadening and the laser 
lineshape, the 1% estimated error in the gas pressure measurement, and the standard devia- 
tion of the individual strengths used to obtain an average line strength for each transition. 
The uncertainty of ±0.3% in the gas cell pathlength was small enough to be neglected. The 
average total error for these line strengths was 2.7%. 

Eighteen transitions were measured at pressures greater than 0.4 torr. Their strengths, 
S[), were corrected for pressure broadening using equation IV-8. To calculate this correc- 
tion one must know the pressure broadening coefficient, 7^. for the line being measured. 
Nine of these eighteen line strengths were corrected using the measured value of 7^ for 
each line. The average value of 7^ for these nine lines was used to calculate the correc- 
tion for the other nine lines, which were recoided at too low a pressure for their pressure 
broadening coefficient to be measured. The resulting line strengths were then increased 
another 3% for the effects of the laser lineshape on the line pro Tile. The 5% correction 
used for the tow pressure measurements was for both pressure broadening and for distor- 
tion of the line profile by the laser so that increasing these eighteen lines by the full 5% 
would result in a slight over-correction for pressure broadening. Therefore, an estimated 
correction for the laser distortion alone of 3% was used, this being the lower bound of 
full 5% ± 27r correction. An uncertainty of ±2%' was also assigned to this 3% correction. 
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In general, though, the total error in these higher pressure measurements was dominated 
by the uncertainty in 7^- The errors in these eighteen strengths were determined the same 
way as the errors in the low pressure strength measurements, except for an extra term 
added to account for the error from the uncertainty in the values of 7^ used to calculate 
the pressure broadening corrections. The average total error for these line strengths was 
6 . 6 %. 

The averaged strengths and their experimental errors are listed in Table lV-5. Their 
values have been normalized to 300 K using the relation 

/ T \ 5/2 - ' t) 0 - e-hi'/kTo) 

S(To- 300K) = S(T) — ) e “ (lV-17) 

\tJ (1 - e-hWkT) 

which can be derived from equation 11-32 and the dependence of Q and N on temperature. 
The partition function Q is primarily dependent on temperature through its rotational 
part, given in equation Vl-2. N, the number of molecules/cm^-atm, is inversely propor- 
tional to temperature. The listed strengths have also been divided by 0.9502 to account 
for the naturally occurring abundance(53) of the sulphur 32 isotope. The average total 
error for all the line strengths in this table is 3.4%, 

We have found the direct method of measuring line strengths to have an accuracy of 
about 3%. However, the relatively high accuracy of the direct method can only be re- 
tained if the line broadening due to the diode laser is less then several percent. Also, the 
gas pressure must be kept very low, or else the pressure broadening coefficient must be 
measured for each transition if the highest possible accuracy is to be achieved. In other 
line strength studies using diode lasers,^'^) low pressure linewidths 8% to 40% greater 
than the Doppler widths have been observed. To keep the error in the strength correction 
for these lines below 5% for example, their linewidths must be measured to an accuracy of 
about 2%, which is difficult to achieve without a large number of measurements, especially 
for heavy molecules with small Doppler widths. Also, the line strength correction factor 
becomes more sensitive to the details of the actual laser lineshape and to the observed ab- 
sorption level as the distortion of the line profile by the laser increases. 
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Table IV -5 


Measured Line Strengths in V2 Band of H2S 


Upper Lower 

K^HlKa. K,,) 

Observed 

Strength 

(cnr*/atm-tin) 

Experimental 

Error 

(%) 

(a)p Factor 

Calculated 

Frequency 

(cm-*) 

6(2,S)-S(t,4) 

i .62C-02 

2.7 

1 44 

1257.068 

6(l,5>-5(2,4) 

5.44E-03 

3.i 

1 .45 

1257.010 

5»4,1>-S(i,4) 

3.11E-04 

3.9 

4.71 

1256.972 

S(2,3>-4(3,2) 

1.13E-02 

3.8 

1 .84 

1256.899 

9<4,6)-9(3,7) 

1.41E-03 

2.8 

3.95 

1256.783 

S(3,2)-4(4,l) 

2.34E-03 

2.4 

1.81 

1256.764 

9(3,6>-9<2,7> 

4.31E-03 

2.7 

4.01 

1256.729 

4(3,2)-3(2,l) 

l,85E-02 

2.5 

2.23 

1253.468 

8(3,6>-8(2,7> 

7.S0E-03 

2.3 

4.43 

1252.874 

8(2,6)-B(l,7) 

2.45E-03 

2.4 

4.34 

1252.862 

8(4,5>-8(3,6) 

6.S4E-03 

2.5 

3.72 

1246.076 

8(3,5)-8(2,6) 

2.29E-03 

3.1 

3.89 

1245.843 

D 4(i,6>*-S(0,5) 

1.66E-02 

2.7 

0.96 

1243.617 

« 6(0,6>-5(i,5> 

5.54E-03 

2.7 

0.96 

1243.616 

7<3,5>-7(2,6) 

3.09E-03 

2.3 

3.64 

1242.557 

7(2,5)-7<l,6) 

9.11E-03 

2.6 

3.57 

1242.500 

4(l,3)-3<2,2) 

5.03E-03 

2.4 

1.52 

1236.727 

8<4,4)-8<3,S) 

1.70E-03 

3.4 

2.92 

1236.650 

8(7,2>-8(6,3> 

2.92E-03 

2.8 

2.61 

1235.705 

8(6,3)-8(S,4) 

4.38E-03 

2.6 

2.84 

1234.613 

5<l,S)-4(0,4) 

6.80E-03 

2.4 

1.14 

1234.585 

5(0,5)-4(i,4) 

2.04E-02 

2.4 

1.14 

1234.579 

7(S,3>-7(4,4> 

2.44E-03 

2.8 

2.70 

1229.354 

5<2,4)-S(i,5) 

3.S3E-03 

2.6 

3.14 

1229.333 

5<l,4>-5<0,5> 

1.07E-02 

2.4 

3.18 

1229.273 

4(i,4)-3(0,3> 

1.91E-02 

2.8 

1.13 

1225.443 

4<0,4)-3(l,3) 

6.40E-03 

3.0 

1.14 

1225.395 

6(4,3>-6(3,4) 

1 . 08E-02 

2.4 

2.65 

1225.161 

6(6,i)-6(5,2) 

5.06E-03 

3.6 

2.45 

1224.333 

3(i,2)-2(2,i) 

7.84E-03 

2.5 

1.46 

1224.279 

S<2,3)-5<1,4) 

1.25E-02 

3.2 

2.49 

1220.649 

6(6>0)-6(S,i> 

1.58E-03 

2 4 

2.22 

1220.589 

2<2,i)-Kl,0> 

1.38E'02 

3.6 

1.46 

1219.949 

4(i,3)-4(0,4) 

3.6SE>I)3 

2.4 

2.59 

1218.947 

S(S,l)-5<4,2> 

2.56E-03 

2.4 

2.26 

1216. >'60 

5<4,2)-5(3,3) 

4.17E-03 

3.3 

2.38 

1216.437 

3<i,3)*2(0,2> 

S.b8E-03 

2.6 

1.23 

1216.297 

3(0,3)>2<i,2> 

1 .59E-02 

2.9 

1.16 

1215.965 

7(S,2)-7<4,3) 

7.30E-03 

2.6 

1.96 

1213.387 

4<3^2)-4<2»3) 

1.37E-02 

3.0 

2.24 

1212.418 

S(’S,*)-S<4,i> 

8.04E-03 

2.7 

2.11 

1211.190 

6(S,l)-6(4,2) 

3.25E-03 

2.5 

1.94 

1209.599 

3(2,2>-'3(l,3> 

3.76E-93 

3.0 

2.28 

1209.523 

5(3»2)-5(2,3) 

1.40E-02 

2.7 

2.08 

1209.082 

3<i,2)-3<l,3) 

l.lSE-92 

2.3 

2.25 

1208.016 
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Table IV-5 (Continued) 


Upper Lower 

J'(Ka'. Kc')-J(Ka. Kc> 

Observed 

Strength 

(cm'Vatm-cm) 

Experimental 

Error 

(%) 

(a)p Factor 

Calculated 
Frequency 
(cm‘* ) 

(b) 

2(i,2>-i<0,i> 

i.2Se-02 

2.3 

1.29 

1207.449 


3(3,i>-3(2,2> 

3.S0E-03 

3.3 

1.84 

1205.043 


4(4,0)-4(3,i) 

3.79E-03 

2.4 

1.78 

1202.825 


4(3,l>-4(2,2> 

5.86E-03 

2.3 

1.66 

1199.530 


2(1,1)*2(0,2) 

3.76E-03 

2.4 

1 .68 

1196.822 


2<2>0)-2<l,l) 

4.9SE-03 

2.S 

1.31 

1191.546 


3(2;l)-3(3,0) 

S.OSE-03 

z.z 

0.56 

1175.161 


5<3,2>-S<4,l) 

3.69E-03 

2.6 

0.47 

1175.088 


l(0,i)-2(i,2) 

6.40E-03 

2.3 

0.78 

1158.169 


2(l,2)-3(0,3) 

9.05E-03 

2.6 

0.81 

1149.771 


2<0,2)-3(i,3> 

2.9iE-03 

2.5 

0.78 

1149.403 


4(0,4)-4ti,3) 

7.12E-0S 

4.4 

0.063 

1148.720 


4<l,4)-4<2,3> 

2.01E-04 

5.0 

0.059 

1148.449 


l(i,0>-2(2,l> 

3.06E-03 

2.5 

0.41 

1147.359 


2<2,1)-3<1,2> 

2.39E-03 

2,6 

0.56 

1144.268 


3(i,3)>4<0,4) 

3.42E-03 

2.5 

0.80 

1140.141 


3(0,3)-4(i»4) 

1.02E-02 

2.4 

0.80 

1140.086 


2<l,i>-3<2,2) 

9.28E-04 

2.5 

0.41 

1138.446 


3(2,2)-4(l,3) 

1.27E-03 

2.4 

0.53 

1132.847 


4<1,4)-5<0,S> 

1.92E-02 

2.3 

0.80 

1130.524 


4<0,4>-5<i,5) 

3.39E-03 

2.4 

0.80 

1130.517 


5<5,0)-&(4,3) 

1.32E-04 

4.2 

0.87 

1129.162 

X 

9(2,7>”9<3,6) 

1.22E-04 

4.8 

0.17 

1128.570 

X 

9<3,7)“9<4,6) 

4.34E-05 

8.3 

0.18 

1128.531 

X 

6<0,6)>6(i,5) 

3.76E-0S 

9.4 

0.062 

1128.366 

X 

6<l,6>-6<2,5) 

i.04E-04 

9.1 

0.057 

1128.359 

X 

2(2,0)-3(3,l> 

4.32E-04 

5,5 

0.16 

1127.345 

X 

4<2,3>'5(l,4> 

3.9SE-03 

2,6 

0.51 

1123.238 


e<i ,7>-8(2»6> 

e.S2E-05 

10,2 

0.22 

1123.085 


B<2*7)-8<3,6) 

2.S1E-04 

10.1 

0.22 

1123.077 


4<l,3>-5<2,4) 

i.29E-03 

2.4 

0.50 

1122.854 


« 5<l,5>-6<0,6> 

3.06E-03 

2.5 

0.80 

1120.813 


I> 5<0,5)>6<i»6) 

9.19E-03 

2.5 

0.80 

1120.812 


4<4,l)-S(3,2> 

3.42E-04 

6.3 

0.39 

1120.591 


i0<2,8)-i0(3^7> 

6.08E-0S 

6.9 

0.45 

1119.500 

X 

10<3,8)-i0<4,7) 

1.72E-04 

6.6 

0.43 

1119.491 

X 

5<2,4)-6<l,5> 

1.19E-03 

2.6 

0.51 

1114.083 


5<l,4)-6<2,5> 

3.63E-03 

2.4 

0.51 

1114.015 


D 6(i,6>-7<0,7> 

7.94E-03 

2.9 

0.82 

1110.992 


♦ 6<0,6)-7<l,7> 

2.65E-03 

2.9 

0.82 

1110.992 


5(2,3>-6<3,4) 

1 .OSL-03 

2.6 

0.24 

1105.720 


6<2,5>-7<l,6) 

3.01E-03 

2.7 

0.51 

1104.935 


6(i,5)-7(2,6) 

9.92E-04 

2.8 

0.51 

1104.923 


S(4,2)>6(3,3> 

i.lSE<04 

3.4 

0.18 

1104.316 

X 

« 12(2,18}*i2(3>9> 

5.20E-05 

9.6 

1.50 

1101.248 

X 

D 12(3,i0>-&2(4,9> 

i.56E-04 

9.6 

1.50 

1101.248 

X 

D 7<0,7>-8(l,8) 

6.02E-03 

2.7 

0.80 

1101.062 


* 7<1,7>“B<0,8> 

2.01E-03 

2.- 

0.80 

1101.062 



(c) 


X 

X 


X 

X 

X 

X 

X 

X 


X 

X 

X 


X 

X 

X 
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Table IV-5 (Continued) 


Upper Lower 

Kc )-J(Ka, K,) 

Observed 

Strength 

(cm"*/atm-cm) 

Experimental 

Error 

(%) 

Factor 

' aiculated 
Frequency 
(cm**) 

(b) 

(c) 

6(3,4>-7<2,S> 

i.0iE'03 

3.2 

0.28 

1098.270 

X 

“x 

* 7<2,6>-8(l,7> 

7.8SE-04 

2.3 

0.53 

1095.699 



D 7(l,6>-8<2,7> 

2.35E*03 

2.3 

0.53 

1095.697 



« 8(0^8>-9(t,?> 

i.46E*03 

2.4 

0.80 

1091.023 



D 8(1,8>>9(0,9> 

4.39E-03 

2.4 

0 . PO 

1091.023 



7(3,5>-8(2,6> 

2.60E-04 

3.7 

0.29 

1089.624 

X 


7<2,5>-8(3,6) 

7.61E-04 

2.4 

0.28 

1089.557 

X 


0 8(2,7>-9(l,8> 

1.67E-03 

2.6 

0.52 

1086.363 



* B<l,7>-9<2,8) 

5.55E-04 

2.6 

0.52 

1086.363 



7(3,4>-8<4,5) 

1.74E-04 

S.O 

0.10 

1082.034 

X 

X 


(a)Observed line strength divided by calculated rigid rotor line strength - see Ch. VI. 
tb>An X is placed in this column if a saturated absorption scan was not available to deter- 
mine the lOO'/r absorption level for the strength measurement. 

(c) An X is placed in this column if the strength was measured at a gas pressure greater than 
0.4 torr. 

D- Denotes stronger line in blended doublet. 

• - Denotes weaker line in blended doublet. This line strength is obtained by di iding the 
stronger line strength by three. 
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In conclusion, this woik has ^own the feasibility of using diode lasers for wide rang- 
ing (>150cnr*) studies of infrared vibrational bands. Relative line frequencies were meas- 
ured to an accuracy of about ±0.0004 cm"* with the diode laser spectrometer. Any sys- 
tematic errors in the frequency measurements were shown to be very small. It was 
demonstrated that line strengths can be measured to an accuracy of about 3% using the 
direct method, if the diode laser linewidth is narrow enough. Corrections *o the measured 
line strengths to account for the distortion of the line profile by the diode laser were cal- 
culated and the importance of using the proper laser lineshape to determine these correc- 
tions was shown. 
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CHAPTER V 


ROTATIONAL ANALYSIS OF THE ^2 BAND OF H 2 S 

V. A. Rotational Structure of the t >2 Band and Line Assignments 

In this chapter, a least-squares analysis of the observed V 2 frequencies using 
Watson’s AS and NS reduced Hamiltonians will be presented. A targe part of this discus- 
sion will center on the characteristics of the AS Hamiltonian in the III'^ representation, 
which was found to fit the transition frequencies about 100 times worse than the other 
reduced Hamiltonians used in this work. First, the tine assignments and some of the gen- 
eral features of the i >2 band of H 2 S will be discussed. 

In order to assign the observed transitions a trial spectrum of the V 2 band was calcu- 
lated with a computer program written by MakK^A)^ which uses Watson’s AS reduced 
Hamiltonian in the I*" representation. This program also provides relative rigid rotor line 
strengths as well as the transition frequencies. Initial values for the band center and for 
A, B, and C in the upper vibrational state were taken from the work of Allen and Plyler,(2) 
who performed the only previous rotational analysis of this band. They were able to 
assign 55 transitions from 1080cm'* to 1257cm** from a spectrum recorded with a resolu- 
tion of 0.3cm**. The ground state constants for this trial calculation were obtained from a 
least-squares fit of the H 2 S pure rotational frequencies^*^) to the AS-F Hamiltonian. 

Terms up to sixth order in J were included in the ground state Hamiltonian. Although a 
better fit of the pure rotational frequencies could be obtained if higher order terms are 
included in the Hamiltonian, these terms are not required for the purpose of making line 
assignments. 

With this first trial calculation about 20 line in the R branch spectrum taken with the 
NBS grating spectrometer could be positively assigned. The measured frequencies for these 
20 transitions were then fit to the AS-F Hamiltonian to get approximate upper state values 
for A. B, C and the five quartic constants, which were in turn used in a second trial 
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calculation of the V 2 band. This calculated spectrum was accurate enough to allow the 
rest of the R branch grating spectrum to be assigned. Another fit to the AS-F Hamilton- 
ian was then made using all of the lines assigned in the R branch grating spectrum so to 
determine as many upper state constants as possible. These constants were used to calcu- 
late one last trial spectrum of the V 2 band in order to assign the H 2 S lines observed in the 
spectra taken with the diode laser spectrometer. Since the J values of the transitions ob- 
served with the diode laser spectrometer were no larger than those in the R branch grating 
spectrum, this calculated spectrum was sufficiently accurate for the assignment of all the 
lines in the diode laser spectra. 

Although the spectrum of a highly asymmetric rotor appears to be a somewhat random 
collection of lines at first glance, there are some regularities in the spectrum which help with 
line assignments. A number of sub-bands with fairly regularly space lines were identified in 
the spectrum of H 2 S. In the R and 0*^^ branches, the six sub-bands 

(J)j j - (J - ' (J)o.J " ’ 

(J)2,j-i '(J - 

were observed with J usually running from one to .>out nine. As J increases, the two sub- 
bands shown together on each line above converge. For intermediate values of J , this gives 
rise to doublets in the spectrum from the two converging sub-bands. The assignment of the 
lines in these doublets to the proper sub-band is easy because the spin statistics of H 2 S require 
the ratio of their strengths to be equal to either 3 to 1 or 1 to 3, depending on the value 
of J. By J = 9 some of these converging sub-bands are separated by less than 0.0001 cm • 
so that only a single line is obscived in the spectrum. Detailed discussions of the general 
structure of asymmetric rotor vibration bands can be lound in the book by Allen and 
Cross.^ 

A stick spectrum of only those transitions whose frequencies were measured in this 
work is shown in Fig. V-1. Below the plot of these transitions is a calculated stick 
spectrum for the band. In thetic plots the length of the stick is proportional to 




HlDN3diS 3AllV13a 
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Figure V-1 ; Stick Spectrum of V2 B^nd of H 2 S with Experimentally Observed Lines Indicated, 

Theoretical Rigid Rotor Strengths Used in this Plot 


tht relative rigid rotor strength for that transition. This spectrum was calculated using the 
AS-1*^ Hamiltonian to sixth order in angular momentum. The rotational constants for this 
Hamiltonian were determined from a least-squares fit of the hand and microwave pure 
rotational frequencies of H2S. Though higher order terms in the ground state Hamiltonian 
are required to fit the microwave rotational data to its experimental accuracy, any changes 
to the stick spectrum from these terms would not be visible. These least-squares fits will 
be discussed in detail later in this chapter. In Fig. V-2 the calculated spectrum of is 
again .shown, but with the P. 0. and R branch lines displaced vertically from each other. 

In this figure the (.Dij - (J - ' ~ sub-bands are indicated to 

illustrate two converging sub-bands. By J = b these two sub-bands arc separated by only 
0.0008cnr*. less than the Doppler widths of the lines. The frequency difference between 
adjacent lines in each of these s ib-bands is about 9.1 cnr*. and is equal to the line separa- 
tions in the RlJ, K = J) series of an oblate symmetric top molecule with 
given by 2(6^^ - = 9.1 cnr*. The gap :n the calculated spectrum at the band 

origin (Vjj 1 182.bcm*l) and the relatively spread out Q branch are the two main features 
that ciiaiacteri/e this band as a type B asymmetric rotor band. When viewing the stick 
spectra in Fig. V-1 and V-2. it should be remembered that the line strengths shown are 
for a rigid rotor while in the experimentally observed spectrum the R branch is much 
stronger than the P branch. 

As discus.sed in Chapter IV. the frequencies of 49 R and Q br. nch lines were meas- 
ured with the NBS gratiirg spectrometer and 77 P, Q. and R branch lines with the diode 
laser spectrometer. The maximum values of J, Kj,, and K^. involved in these transitions 
were 10. 7. and 9 respectively. These line frequencies wore fit to the AS and NS forms 
of Watson’s rotational Hamiltonian in both the F and lll*^ representations using a computer 
program written by Kirshner.l-* ^ Conventionally the (or I**) representation is used for a 
prolate a.syminetric top (-1 < k < Ot and the III*^ (or lll^) representation tor an oblate 
asymmetric top (0 < k < I ). The representation is usually chosen in this way because as k 
approaches -1 the otT-diagonal matrix elements of the rigid rotor Hamiltonian become 
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FREQUENCY (cm ') 

Figure V-2. Stick Spectrum of ^2 Band of HjS, Theoretical Rigid Rotor Strengths Uied in this Plot 


very small for the K representation, and as k approaches -t-l they become very small for 
the III*’ representation. Since k 0.5 for H 2 S. the two most recent studies of HiS rota- 
tional spectra used either the Hl’tn) or the representations. For a molecule as 

asymmetric as H 2 S. though, the choice of representation would appear *0 be somewhat 
arbitrary. Since a program was available that could fit both reduced Hatniltonians in 
either representation, fits wen? made using all four of the Hamiltonian fonn AS -•/, 
AS-mr. NS-K, and NS-llir 

To predict the t »2 band transition frequencies one needs both the ground and upper 
vibrational state Hamiltonian cx^nstants. In this study the ground state constants were de- 
termined from the very accurate (=»0. 1 MHz) pure rotational transition frequencies reported 
by Helminger et al.(*^) and from the ground state combination differences (GSCD) ob- 
tained from the 1^2 band transitions. A ground state combination difference is the magni- 
tude of the frequency difference between two lines in a vibrational band with common 
upper states. Tb's frequency difference is equivalent to the frequency of a pure rotational 
transition between the ganind state levels of these two lines. The pure rotational mica> 
wave transitions had maximum values of J. K^. uid K^. of ‘J. 7. and 7 respectively. Best fit 
ground state constants were determined for all four Hamiltonian forms. These ground 
state constants were then used to calculate the lower state energy levels of the band 
transitions in order to obtain the upper state eneigy levels required for the least-squares 
fit of the upper state Hamiltonian constants. 

V. B. Least-Squares Fitting of the Rotational Hamiltonian 

Least-squares fits of the rotational Hamiltonian to the transition frequencies were per- 
formed with a computer program written by Kirschner^’l )and extended by Hillman.^^^^ The 
methods used in this program are discussed in detail by Kirshner and will only be summarized 
here, stressing those points important for the interpretation of our results. 

We wish to find the set of statistically significant constants, C j. that can most accurately 
predict the measurvHi frequencies from the Hamiltonian Hr = 1 Cj Op*, where Op* is an 
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angular momentum operator. Since the calculated frequencies, are not linear functions 
of the Cj, an iterative approach is used in the least-squares fitting routine. This is done by 
making initial guesses for the C, and then finding the set of corrections, 6Cj, that best solve 
the set of linear equations 

yi = 2xki6Ck (V-1) 

where yj = v?*’® - Xy = and is the observed frequency of the i*^ transi- 

tion. In this equation i runs over the N transitions and k over the M Hamiltonian constants. 
The best 5Cj, in a least-squares sense, are those that minimize the goodness of fit expression, 

X2 = I Wi(yi-2xy5Ck)2 (V-2) 

1=1 k 

where Wj = 1 /a? is the weight of the i**' transition, and Oj is the experimental error assigned 
to that transition. When the best SCj are found, the original values of the Cj are corrected and 
another set of corrections to these new Cj are found. This process is continued until the 6Cj 
satisfy the relation l6CjK 10“2 a(Cj), where a(Cj) is the standard deviation for the constant 

q- 

For a vibrational transition where the ground state constants are kept fixed during the 
fit we have, for example, 

(V-3) 

where | i > is the i^ upper state asymmetric top eigenfunction which satisfies 

E[‘PP" = <i|HR li>= (ilOp*"!!). (V-4) 

For a pure rotational transition we have 

x,^j= — - — = <i| Op*^ li>- <i' lOp*^ li'> (V-5) 

dCk 

where the upper state energy is 

Eupper = ( i I Hr I i > (V-6) 
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and the lower state energy is 


glower = iHj^ I i'). (V-7) 

To find the set of 5C| that minimize one usually generates the normal equations from 
the condition By^/b(SC^) = 0. The normal equations are 

Y = A 6C (V-8) 

w^ere 

'^k = ? Wj yj 

and 

Aki=2WjXj^Xji. 

The corrections 5C^ are found by inverll.ig the normal matrix A to give 

8C = A-‘ Y. (V-9) 

In practice, difficulties can arise in the inversion of A when a large number of 5Cj are 
being determined, because the determinant of A can become very small. This can be due to 
near linear dependencies in the independent variables (the expectation values of the angular 
momentum operators) or from a lack of sufficient experimental data to determine some ot 
the 6Cj independently of the others. Even in cases where near linear dependencies do not 
exist round-off errors incurred when inverting the normal matrix can be significant. 

To overcome this difficulty Kirshner used an algorithm called ORTHO^^^) that avoids 
the need for matrix inversions in the least squares calculation. We originally wanted to find 
the 5Cj that best solve the set of equations V-1 , which can be written 

l = XBC, (V-10) 

but problems in the inversion of the normal matrix A arise from near linear dependences in 
the columns of X, or between the vectors^. To treat this problem, the set of M vectors ^ 
can be transformed to an orthonormal set of vectors using the Gram-Schmidt orthonormal- 
ization procedure. The transformation of X is given by E = XD. We now must solve the set 
of equations 
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)^ = e 6 C' (V-ii) 

to get 

6 C' = eT^, (V-12) 

which can be written this way since the£j are orthonormal and therefore . The 

original corrections, 6 C, can easily be shown to be given by 

6 C = AeT^. (V-13) 

The 5C' obtained tf is way are the best least-squares estimates for the 6 C'. Because the f ; 
vectors are orthonormal, the 6 C| are their Fourier coefficients, which, it can be shown,C^^) 
are identical to the best least-squares estimates for the 5C|. Therefore a matrix inversion is 
avoided and only the Gram-Schmidt orthonormalization matrix A and the vectors e: need be 
calculated. The weighting of the yj can be accomplished by scaling the yj and xj^ byy/~^. 

A natural product of the orthonormalization routine in ORTHO is information on the 
orthogonality of the x: vectors, which indicates how independently the 6 C; can be deter- 
mined from the measured transition frequencies. The quantity found by ORTHO is the ratio 
GFj/ENFj where ENFj is the magnitude of the vector Xj and GFj. called the Gram Factor of 
the orthonormal vector ej, is the magnitude of the component of Xj orthogonal to x i . x ■> , 

X 7 . . . . Xii . In our case the vector Xi is associated with the constant X. ,':t with Y, xi with 
Z, X 4 with Aj or Dj, and so on in the order the exjnstants are given in equations 11-24 for the 
AS form or 11-27 for the NS form. A small GFj/ENFj ratio indicates a near linear dependence 
of X: on the previous i - 1 vectors while a GFj/ENF; ratio close to one means that X; is nearly 
linearly independent from the previous vectors. The over-all degree of linear dependence in 
the fit is given by the Gram Determinant. 

GMDT = n (GFi/ENFi)2. (V-14) 

Quantities determined for each fit were the best fit values for the Hamiltonian coeffi- 
cients, Cj; the standard deviation of the coefficients, o(Cj); the correlation matrix. the 
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calculated frequencies, ; the frequency residuals, the standard deviation of the 

residuals, Or ; and the weighted standard deviation of the fit, Oj,- . The standard deviation o'" 
the residuals, hereafter often referred to as the standard deviation of the fit, is given by 

1 N \l/2 

2 (ycalc_^pbs)2 ] (V-15) 

(N-M) i=l ' ‘ / 

Oj:, which will always be referred to as the weighted standard deviation of the fit, is given by 

and for a perfect fit should equal one. For more rigorous definitions of the above quantities 
see Kirshnei^21) or the book by Hamilton.^^^) 

In order to obtain the best least -squares fit for a given data set and Hamiltonian, the set 
of constants (and their associated operators) included in the Hamiltonian were varied until all 
statistically significant constants (Cj > 2o(C^)) were used, and hopefully the standard devia- 
tion of the residuals became approximately equal to the average experimental error. Gener- 
ally, quite a number of fits (>10), using different combinations of constants in the Hamil- 
tonian, had to be made until it was clear that the standard deviation of the residuals could 
not be lowered by adding another term in the Hamiltonian that would in turn be statistically 
significant. 

V. C. Results of Least-Squares Fits for the Rotational Constants of the V 2 Band 
The standard deviations of the residuals for the best ground and upper state fits are 
given in Table V-1 . Also included in this table are the standard deviations for fits that only 
included terms in the Hamiltonian up to second, fourth, and sixth power in angular momen- 
tum. The fitting program used here contains all Hamiltonian terms up to eighth power for 
the AS form but does not contain the off-diagonal eighth power terms for the NS Hamil- 
tonian. The program was not modified to include the eighth order off diagonal terms in the 
NS Hamiltonian because good fits of the data were obtained with the AS-F Hamiltonian. 

No eighth order terms were required for any of the upper state fits. 
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Table V-1 

Standard Deviation of Ground and Upper State Fits 
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In Tables V-2 through V-9 the best fit molecular constants for all four Hamiltonian 
forms for both the ground state and the first excited state of the ^2 vibration are given. The 
microwave transitions were weighted using the uncertainties in their frequencies given by Hel- 
minger et. al.(58) Xhe weights for the GSCD and the infrared transitions were determined from 
the experimental errors given in Table lV-4 for the infrared frequencies. The errors quoted for 
the molecular constants equal twice their standard deviations, which corresponds to 95% confi- 
dence intervals. From Table V-1 we see that fits with standard deviations comparable to the 
experimental errors were obtained for the ground state with the AS-F Hamiltonian and for 
the upper vibrational state with the AS-F, NS-F, and NS-IIF Hamilto.aians. The weighted 
standard deviations of the fit for these three upper state Hamiltonian forms were all close to 
one. It is not surprising that the microwave frequencies could not be fit to the experimental 
accuracy of about 0.2 MHz with the NS-F and NS-IIF Hamiltonian forms since they did not 
contain any eighth order off-diagonal terms. These terms had to be included in the AS-F 
Hamiltonian to get a fit comparable to the accuracy of the data. If the Hamiltonian is trun- 
cated to sixth order, the NS-I^ and NS-llF forms fit the microwave data as well as the AS-F 
form, so there appears to be no fundamental inaccuracy in the NS form Hamiltonians. 

It was very surprising though, that the fits using the AS-IIF form had such large stand- 
ard deviations. When all four forms are truncated to fourth power, the standard deviation of 
the residuals is about 1 3 times greater for the ground state AS-IIF Hamiltonian and about 21 
times greater for the upper state AS-IIF Hamiltonian than for the other three ground and 
upper state Hamiltonian forms. When the Hamiltonians are extended to sixth power the 
standard deviation of the AS-llF form fit is greater than the standard deviations for the other 
three forms by a factor of 68 for the ground state and 94 for the upper state. The best eighth 
order ground state fit using the AS-IIF form has a standard deviation about 100 times larger 
than the experimental accuracy of the microwave data. For the upper state fit, where only 
terms up to sixth order were found to be statistically significant, the standard deviation of 
AS-IIF form fit was again about 100 times larger than the accuracy of the infrared frequen- 
cies. As would be expected, the standard deviations of the AS-llF Hamiltonian constants 
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Table V-2 


Molecular Constants for Ground State, AS Form, F Representation 




Constant (MHz) 

2a (MHz) 

< 

< 


310590.038 

0.032 

fiAS 


270357.131 

0.038 

CAS 


141825.442 

0.030 



19.5708 

0.0025 

Ajk 


- 68.3542 

0.0089 



111.0269 

0.0086 



8.8627 

0.0015 

«K 


- 3.9769 

0.0052 

Hj 

X 10^ 

8.038 

0.061 

Hjk 

tf 

- 45.17 

0.25 

Hkj 

9t 

35.55 

0.84 

Hr 

tf 

43.37 

0.88 

hj 

tt 

- 4.159 

0.051 


ft 

- 14.71 

0.13 

hR 

ft 

36.96 

0.34 

^JR 

X 10 ^ 

71 

10 

^RJ 

tt 

-108 

13 


tt 

- 3.29 

0..‘’4 

«JR 

tt 

35.35 

0.83 

®RJ 

tt 

- 93.2 

4.1 


tt 

31.2 

2.1 


Standard Deviation of Microwave Transitions = 0.1 5 MHz 
Standard Deviation ofGSCD = 84 MHz 
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Table V-3 


Molecular Constants for Ground State, AS Form, Representation 




Constant (MHz) 

2o (MHz) 

aAS 


310192.7 

2.4 

bAS 


270875.0 

2.2 

(-.4S 


141706.2 

2.1 



50.51 

0.30 



-160.79 

0.30 



113.07 

0.41 



6.056 

0.024 



-256.63 

0.40 

Hj 

X 10^ 

99 

13 

Hjk 

tt 

-351 

24 

Hrj 

tf 

675 

48 

Hr 

ft 

-399 

30 


ft 

4.31 

0.94 

^JR 

•t 

-170 

16 


ft 

-2003 

60 


X 1Q6 

-1558 

182 

Ljjr 

ft 

1260 

132 


ft 

24.1 

8.0 

*^RJ 

tf 

-6622 

1600 

h 

It 

24948 

2400 


Standard Deviation of Microwave T ransitions = 1 4.6 MHz 
Standard Deviation of GSCD = 2572 MHz 
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Table V-4 


Molecular Constants for Ground State, NS Form, Representation 




Constant (MHz) 

2o (MHz) 

aNS 


310583.521 

0.084 

bNS 


270367.70 

0.12 

CNS 


141820.05 

0.077 

Dj 


20.8654 

0.0073 



-76.239 

0.011 

Dk 


117.717 

0.016 



-8.8658 

0.0044 



0.6433 

0.0013 

H'j 

X 10^ 

10.33 

0.13 

H'jk 

ft 

-90.01 

0.33 

h'kj 

ft 

154.0 

1.1 

Hr 

ft 

-33.1 

1.4 

hi 

ft 

2.83 

0.13 

hj 

tf 

-1.014 

0.032 

hj 

tf 

1.187 

0.026 

L'j 

X 10^ 

-6.9 

2.1 

LjJK 

ft 

62.9 

5.9 


ft 

-196 

22 

l'kkj 

tf 

282 

21 

l'k 

ft 

-175.0 

9.1 


Standard Deviation of Microwave Transitions = 0 49 MHz 
Standard Deviation of GSCD = 84 MHz 
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Molecular Constant! fo?^rpund ^a^|fi£i@flh, 




‘ — 

Constant ^Hz) 

2o (MHz) 

aNS 


310667.57 

ms 

bNS 


270331.55 

QA2 

cNS 


Ul 789.59 

o.n 



33.130 

OJ031 

^JK 



0.021 

Dk 


27.831 

0.032 

5l 


-^.2647 

0.004Q 

*2 


-8.3348 

0.0018 



13.92 

0.5S 

B'jk 

H 

-48 

QM 

h'kj 

« 

57.9 

2.S 

M'k 

It 

-23.3 

2j0 

hi 


5.185 

OM 

h2 

tt 

6.819 

0M6 

1-3 

ft 

0.928 

omo 


Standard Deviatiofi of M krowaffro T 2J3BifiM 
Standard Deviatioa of GSCD = 77 ISHz 
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Table V--6 


Molecular C'oiistanis for Upper State of Band. AS Form. 
F Representation 



Constant (enr* > 

2o icnr* ) 

aas 

io.7::i: 

0.00022 


^.::43b 

0.00014 

(AS 

4.668S6h 

0.(XHK)62 

% xuy» 

7.550 

0.033 


-:7,35 

0.14 

^ ” 

45.56 

0.10 

S,1 

3.457 

0.017 

% ” 

-0.144 

0.06 1 

H| X llF’ 

0.377 

0.034 

h,k ' 

-1.00 

0. i 8 


l.lo 

0.34 

% 

2.36 

0.26 


0.158 

0.016 


-0.428 

0.082 

% 

1.756 

0.08C 


1 182.57006 

0.tXX)43 


SteadaM Itevwtion ol Residuals - 0.001 8 enr ’ 
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Table V-7 


Molecnlar Constants for Upper State of v. Band, AS Form, lll'^ 

Representation 



Constant (cm”* ) 

2o (cm‘* ) 

^AS 

10.696 

0.012 

rAS 

9.268 

0.011 

^AS 

4.6445 

0.0090 

Aj X 10^ 

25.9 

3.2 

A ** 

^JK 

-89 

11 


49.6 

8.2 

«5j 

0 


«K 

-136 

24 

Hj X 10^ 

6.9 

2.1 

Hjk 

0 


Hkj 

-119 

22 

Hk 

97 

21 

■hj 

0 

- 

hjK 

-76 

21 

hK 

-I7I 

41 

"0 

1 182.583 

0.057 


*These constants were not statistically significant when included in the fit 
and were set equal to zero in the fit for the constants given above. 

Standard Deviation of Residuals = 0. 1 70 cm‘^ 
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Table V-8 


Molecular Constants for Upper State of ^2, Band, NS Form. 1*^ Representation 




Constant (cm'* ) 

2o (cm“^) 

aNS 


10.72208 

0.00019 

rNS 


9.22444 

0.00013 

(-NS 


4.668764 

0.000062 

Dj X 

lO** 

7.633 

0.030 

Djk X 

tf 

-27.96 

0.17 


tt 

46.17 

0.20 

6j X 

ft 

-3.480 

0.016 

6-, X 

If 

0.0524 

0.0084 

Hj X 

10^ 

0.368 

0.030 

Hjj, X 

It 

-3.35 

0.22 

^KJ ^ 

tf 

5.67 

0.42 

Hk 

It 

-0.62 

0.34 

hj X 

tl 

0.138 

0.015 

h, X 

ft 

-0.028 

0.011 

X 

tl 

0.0539 

0.0020 



1 182.57703 

0.00042 


Standard Deviation of Residuals = 0.0018 cm~* 
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Table V-9 


Molecular Constants for Upper State of Band, NS Form, IlF Representation 




Constant (cm"^) 

2a (cm“b 

aNS 


10.72538 

0.00019 

bNS 


9.22311 

0.00013 

CNS 


4.667511 

0.000058 

Dj X 

10^ 

13.594 

0.036 

^JK ^ 

u 

-25.805 

0.072 


u 

12.814 

0.058 

6j X 

i4 

-2.842 

0.040 

82 X 

It 

-3.292 

0.015 

X 

10^ 

0.725 

0.040 

^JK ^ 

It 

-2.65 

0.10 

H|,jx 

It 

3.26 

0.12 

Hk 

It 

-1.345 

0.048 

hj X 

tl 

0.291 

0.052 

hi X 

tl 

0.326 

0.028 

X 

It 

0.0477 

0.0056 

>^0 


1182.57694 

0.00041 


Standard Deviation of Residuals = 0.0018 cm‘* 
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i 

are much lai^ier than those for the other Hamiltonians. Also, the off-diagonal quartic con- 
stants and all higher order •'onstants in the AS-IIP form are. on the average, larger than the 
corresponding constants in the other forms. These AS-IIH fits are so poor that errors or 
numerical problems in the fitting program were at first suspected to be responsible for these 
results. The inability to obtain a good upper state fit is also accentuated by using such poor 
ground state constants to calculate the lower energy levels of the vibrational transitions. 

Later we will show that these results can be related directly to the unitary transformation used 
to obtain the AS-IIF Hamiltonian. Before discussing causes for the poor AS-HF form fits 
some observations on the fitting results in general will be made. 

The rigid rotor Hamiltonian was, of course, unable to fit the data even close to the ex- 
perimental accuracy because of the large centrifugal distortion corrections required for the 
higher J transitions. The rigid rotor fits were identical for all four Hamiltonian forms, as they 
should be. since is symmetric in X. Y. and Z and has the same form for both the AS and 
NS reduced Hamiltonians. When centrifugal distortion terms are included in Hr the accuracy 
of the fits varies among the Hamiltonian forms. For a given reduction, variations with repre- 
sentation are expected since the set of operators in each order is not symmetric in X, Y. and 
Z. Also, the reduction of the Hamiltonian transforms varying amounts of the original coeffi- 
cients to the new higher order coefficients depending upon the value of sj \ \ . whose value in 
turn depends ujx)n both the representation and the particular reduction being used. There- 
fore. if the Hamiltonian is truncated to a given power of angular momentum, it will not fit the 
data identically for different reductions and representations. If the Hamiltonian included 
enough higher order terms, the fits for different Hamiltonian forms would be identical 
since all the effects of lower order coefficients transformed to higher order coefficients 
would then be contained in the calculation. In our case, except for tfce AS-UF form, the 
standard deviations cf the fits for different forms with the Hamiltonian truncated to a 
given order are about the same. This indicates that the effects of the unitary transformations 
performed on the original Hamiltonian for these three reduced forms are very similar. 
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The ground state Hamiltonians were fit using both the microwave transitions alone and 
this data combined with the ground state combination difl'erences from the band. The 
GSCD had very little effect on the eighth order fits, usually changing constants obtained from 
a fit to only the microwave data by less than their standard deviations. This was due to the 
tow weights assigned to the GSCD compared to the microwave transitions and because the J. 
K;,. and values of the GSCD were about the same as those for the microwave transitions. 
In fact, the 51 GSCD derived from the I'l band contained only eleven energy levels not al- 
ready present in the microwave data. The microwave pure nitational frequencies are about 
500 times more accurate than the GSCD and were therefore weighted a factor of (500)^ = 

2.5 X 10'* times heavier than the GSCD. As seen in Table V-1 the standard deviations of the 
residuals for the GSCD are about equal to their experimental uncertainty, except for the 
AS-lD'’ Hamiltonian. This means that the microwave measurements and the band GSCD 
are consistent with each other to at least the accuracy of the GSCD. 

The weighted standard deviations of the fit are greater than one, as expected, for all the 
best ground state fits except the AS-F form. When both the microwave transitions and the 
GSCD are used in the fit of this Hamiltonian form the weighted standard deviation of the fit 
is 0.84. This figure is somewhat misleading since it is dominated by the residuals of the GSCD, 
which contribute very little to the determination of the ground stale constants. If a fit is per- 
formed with the AS-F form using only the microwave transitions the weighted standaRl devi- 
ation of the fit is reduced to 0.45. Therefore, the fit is “better” than it really should be con- 
sidering the quality of the data. Since the microwave data set contains only 39 transitions, 
which are used to determine 21 statistically significant rotational constants, there arc only 
N - M = 18 degrees of freedom m the fit. This small weighted standard deviation of the fit 
coupled with so few degrees of freedom suggests that the fit for the higlier order constants 
is somewhat parametne in nature. Tliis is probably also true when the GSCD are added to 
the data set. since they arc weiglited so lightly and affect the constants so little. 

Still, including the higlier order statistically significant constants in the fit is believed to 
be reasonable since they provide correc t'ons to the energy levels that cannot be handled by 


the lower order terms in the Hamiltonian. If these corrections are left out of the Hamiltonian, 
the fitting process will minimize the standard deviation by trying to fit the higher order be- 
havior of the Hamiltonian with lower order terms. This will lead to less accurate lower order 
constants than if the higher order behavior of the Hamiltonian is handled by these higher 
order constants, even if the result is a parametric type fit yielding values for the higher order 
constants that do not have physical significance. These comments should be kept in mind 
when interpreting the standard deviations of the constants obtained in these fits. To improve 
the molecular constants given here any further, transitions with higher values of J and Kg 
must be included in the fits. 

The correlation tables for the ground and upper state fits are given in Appendix B. As 
is usually the case for fits to a large number of Hamiltonian constants, there are a number of 
high correlations between constants. The best ground state fit. which used the AS-H Hamil- 
tonian. had a total of 15 correlations with lp|j l> 0.9. ranging in magnitude from 0.902 to 
0.986, with an average value of 0.94. There were fewer high correlations among the 
upper state constants. The AS-F form had only 3 correlations between upper state constants 
with |pjj|>0.9, and none of the forms had more than 5 of these high correlations. The AS-F 
form grounu state had 1 1 high correlations among the sixth and lower order constants com- 
pared to 3 high correlations among all the upper state constants. The upper state probably has 
fewer high correlations because more transitions were used in the fit for this state than for the 
ground state, making it easier to determine all the upper state c-onstants independently. Some 
high correlations are also expected due to the approximate validity of the planarity relations 
among the quartic and sextic constants. 

Comparisons of our AS-IIF form ground state fits to others will be made in the next 
section. The only other fit of the H^S microwave data not using the AS-llF form was by 
Gillis and Kdwardsf who used Typke’s Hamiltonian to sixth order in the IlF representation. 
Tliey included both the microwave data and GSCD from several vibrational bands in their fit 
to obtain the constants given in Table V-10. As previously discussed. Typke’s Hamiltonian is 
very similar to Watson s NS reduced form. ITie main differences between them will be 


107 



Table V-10 


Molecular Constants for Ground State, Typke’s Hamiltonian, IlF Representation - 

From Gillis and Edwards^ 




Constant (MHz) 

6a (MHz) 


A' 


310630.0 

9.3 

B' 


270292.9 

8.7 

C' 


141837.8 

3.9 

D'j 


32.98 

0.24 

D'jk 


-58.92 

0.39 

d'k 


27.79 

0.24 



-6.25 

0.29 



8.33 

0.21 

H'j 

X 10^ 

12.7 

2.8 

H'jk 

n 

^6.0 

7.8 

h'kj 

n 

56 

10 

Hr 

ft 

-22.6 

4.5 

H's 

ft 

9.7 

3.9 

H'e 

ft 

26.5 

4.8 

H'i6 

ft 

7.2 

2.1 


Standard Deviation of Microwave Transitions = 3 MHz 
Standard Deviation of GSCD = 90 MHz 
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noticeable in higli order (sixth and greater) off-diagonal constants, as can be seen by compar- 
ing our NS-llF form ground state constants in Table V-5 to their constants. The standard 
deviations of 3 MHz for the microwave data and 90 MHz for the GSCD obtained by Gillis 
and Edwards are very similar to our values of 2.1 MHz and 77 MHz for these quantities using 
the NS-IIF form Hamiltonian. 

We conclude that the best rotational constants for both the ground and upper state of 
the Pj band of H^S are those associated with the AS-H Hamiltonians. The observed Pj band 
transition frequencies minus the frequencies calculated with this Hamiltonian are given in 
Table IV-4. The AS-H form ground and upper state constants given in Tables V-2 and V-6 
were used to calculate these transition frequencies. Even though three of the upper state Ham- 
iltonian forms could fit the infrared data to experimental accuracies, the AS-H form constants 
should be more accurate since this form had the best ground state constants. The AS-IIH Ham- 
iltonian, though, was found to produce fits 10 to 100 times worse than the other three forms. 
Since this form was used by Helminger et. al. in the original fit of the microwave frequencies 
and because the IlH (or III®) representation is considered to be the natural representation to 
use for H,S, it is important to understand why this Hamiltonian produced such poor results. 
The following section will discuss the AS-lIH fits in detail and the cause for their poor 
behavior. 

V. D. Characteristics of the AS-llE Hamiltonian 

In this section, our results for the AS-llF form ground state fits will be compared to 
those of Helminger et and Gillis and Edwards.!*^) to see if they are consistent in order 
to help rule out the possibility that program errors are responsible for our poor AS-llE fits. 
Then the determinable quadratic and quartic constants for all four Hamiltonian forms will 
be compared to search for any obvious errors in the AS-IIF results. We will also show that 
the AS-HF determinable constants are converging to their correct values (taken to be the 
determinable constants for the AS-F Hamiltonian) as higher order terms are added to the 
Hamiltonian. Further evidence will then be given to support the conclusion that the poor 
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performance of the AS-IIK Hamiltonian is due to its very slow convergence. This will in- 
clude a quantitative study of the unitary transformation that produces the AS-IIH 
Hamiltonian. 

Helminger et. al. fit the AS-III® Hamiltonian to the 39 microwave transitions with a 
standard deviation of 0.38 MHz. They used a total of 24 constants in the Hamiltonian, in- 
cluding 5 tenth order ones, which are given in Table V-1 1 . See reference ( 1 6) for the defini- 
tions of the tenth order constants given in this table. As expected, their eighth and lower 
order constants had standard deviations somewhat lower than ours for the AS-llF 
form, but much larger than the standard deviations of our AS-F form constants. Since the 
standard deviations of our fits to the AS-IIF form improve rapidly as higher order terms 
are added, our results appear consistent with those of Helminger et. al. 

Gillis and Edwards also fit the 39 microwave transitions along with a large number of 
GSCD from several vibrational bands to the AS-IIF Hamiltonian. The constants from their 
fit are given in Table V-1 2. Their fit had a standard deviation of 1 80 MHz for the microwave 
data and 330 MHz for the GSCD when using sixth order terms and lower in the Hamiltonian. 
This can be compared to our standard deviations of 136 MHz and 3241 MHz for the micro- 
wave data and GSCD respectively for the AS-IIF Hamiltonian truncated at sixth order. Gillis 
and Edward’s GSCD had a much lower standard deviation because they were weighted only 
60 times lighter than the microwave transitions, while ours were weighted 2.5 X 10^ times 
lighter than the microwave data. Gillis and Edwards weiglited their GSCD about 1000 times 
heavier than the experimental accu»«cy would suggest because they fit so poorly when 
weighted properly. These heavier weights resulted in a better fit to the CSCD but at the 
same time led to a poorer fit to the microwave data. Therefore it is not surprising that their 
GSCD had a lower standard deviation than ours and their microwave transitions a higher 
standard deviation. This arbitral^' weighting scheme, while improving the fit to the GSCD, 
still results in a standard deviation for the microwave data 100 times worse than the experi- 
mental error. This highly unsatisfactory situation, resulting from their weighting scheme, 
indicates that a different approach should be taken to either improve the fit or understand 
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Table V-1 1 


Molecular Constants for Ground State, AS Form, III® Representation - From 

Helminger et. al.(*^) 




Constant (MHz) 

la (MHz) 

aAS 


310182.24 

0.60 

rAS 


270844.05 

0.51 

CAS 


141705.88 

0.51 



49.851 

0.038 

MU 


-159.696 

0.069 



111.851 

0.068 



-6.0191 

0.0050 

5k 


262.17 

0.21 

Hj 

X 10^ 

28.13 

1.05 


ft 

-228.3 

2.7 

Hkj 

ft 

459 

16 

Hk 

ff 

-276 

14 

hj 

ft 

-5.841 

0.081 

hjK 

ft 

242.8 

5.2 

hK 

ft 

287.0 

39 

^KKJ 

X 10^ 

2199 

1190 

^K 

ft 

-2343 

1260 

^KJ 

a 

-5327 

990 


ft 

-41648 

1900 

Pkj 

X 1Q6 

-304 

70 

**KKJ 

ff 

803 

180 

Pk 

ft 

-515 

114 

PjK 

ft 

-12.3 

5.2 

pkkj 

ft 

130 

13 


Standard Deviation of Microwave Transitions = 0.38 MHz 
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Table V-12 


Molecular Constants for Ground State, AS Form, III*^ Representation - From 

Gillis and Edward^ 




Constant (MHz) 

6o (MHz) 

aAS 


310172 

48 

fiAS 


270826 

42 

CAS 


141672 

23 



48.6 

1.1 



-161.6 

3.3 



115.8 

2.5 



5.6 

1.0 



-241.3 

1.5 

Hj 

X 10^ 

17 

14 

Hjk 


-708 

78 

Hrj 

t1 

2042 

130 

Hk 

rr 

-1357 

66 

hj 

tt 

7.3 

6.0 

^JK 

ff 

-298 

57 


tt 

305 

63 


Standard Deviation of Microwave Transitions = 180 MHz 
Standard Deviation of GSCD = 330 MHz 
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its behavior. Since the other Hamiltonian forms can fit the GSCD so much better than the 
AS-III'' form while using the proper weights, a different weighting scheme for this Hamil- 
tonian cannot be justified. 

A direct comparison of the quadratic and quartic constants obtained for each Hamilton- 
ian form can be made by transforming them into the determinable constants which are. to a 
good approximation, independent of Hamiltonian form. The relations for the determinable 
constants in terms of the AS and NS form constants have been derived by Watsonl^l and 
are given in Appendix A. The quartic determinable constants were introduced in equation 
11-18, where they are written in terns of the ^a/3 coefficients. The calculated determinable 
constants for the best ground and upper state fits are given in Tables V-1 3 to V-l 5. Table 
V-1-.* gives the upper state determinable constants in cm‘l, while Table V-l 5 gives these same 
constants in MH/. so they can be directly compared to the ground state determmablc constants 
given in Table V-l 3. As expected, the agreement among the determinable constants for the 
AS-l'". NS-f, and NS-llF forms is much better than between these a>nstants and the AS-lll’’ 
determinable constants. For the more precise ground state fits the AS-lll’’ determinable con- 
stants are in much better agreement with the others than foi the upper state fits. Since a 
higher order fit to more accuiate data leads to better deter:iiinable constants in the AS-llP^ 
form, it is unlikely that the relatively poor fits using this form are a result of program errors. 
If the AS-llH Hamiltonian converges more slowly than the other Hamiltonian forms, though, 
this result is understandable. 

Also supportive of our conclusion that the AS-111'^ Hamiltonian converges more slowly 
tlian the other Hamiltonian forms are the (IF/FNF values calculated for each fit performed 
m this work. Their values, along with the Ciram Determinants, are given in Table V-lb for 
the ground state and in Table V-i / for the upper stale fits. The Gram Determinant, which 
IS an over-all measure of the linear dependence o! the ilainillonian terms, is about eleven 
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Table V-13 


Determinable Constants for Best Ground State Fits 


Constant 


Hamiltonian Form 


AS, F 

AS, IlF 

NS, F 

NS, IlF 

aD 

310629.18 

310635.46 

310629.11 

310628.83 


270318.15 

270315.98 

270318.02 

270317.86 

c’^ 

141806.00 

141808.52 

141805.85 

141805.84 

^aa 

-62.24 

-62.80 

-62.34 

-62.32 

Tbb 

-37.30 

-38.59 

-37.31 

-37.28 

Tec 

-1.845 

-2.894 

-1.847 

-1.826 

Ti 

9.642 

8.662 

9.784 

9.753 

*T2 

5.683 X 10'^ 

-0.4733 X 10^ 

5.701 X 10^ 

5.761 X 10^ 


*Uniis forX 2 are (MHz)^, all others are in MHz. 
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Table V-14 


Detenninable Constants for Best Upper State Fits 


Constant 


Hamiltonian Form 


AS, 1^ 

<*>A... 111^ 

NS. 1^ 

NS, r.i'^ 

aD 

10.72364 

10.71948 

10.72364 

10.72363 


9.22247 

9.23708 

9.22250 

9.22250 

C^^ 

4.66831 

4.63948 

4.66821 

4.66825 

T,.(x !03) 

- 2.511 

-2.590 

-2.585 

-2.586 

Tbb 

-1.447 

-2.590 

-1.449 

-1.449 

T^.,. (X 

-6.450 

135.0 

-5.687 

-6.033 

o 

X 

— 

4.668 

11.3 

4.743 

4.772 

^^»>T2(x 10-^) 

1.8'>5 

4.916 

1.876 

1.921 


^“^Note that for the AS-lll*^ form the quartic constant 5j = 0. 
^^^Units for T •, are (cm" * , all others are in cm" * . 
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U)TabIc V-15 


Determinable Constants for Best Upper State Fits 


Constant 


Hamiltonian Form 


AS, K 

(b>AS. 111^ 

NS. r 

NS, IIF 

Al> 

321486.7 

321362.0 

321486.7 

321486.4 

bd 

276482.7 

276920.7 

276483.6 

276483.6 

CD 

13 ^) 952.4 

139088.1 

139949.4 

1 39950.6 

Taa 

- 77.26 

- 77.65 

- 77.50 

- 77.53 

Tbb 

- 43.38 

- 77.65 

- 43.44 

- 43.44 

'T 

•cc 

- 1.934 

40.47 

- 1.705 

- 1.809 

Tl 

13.99 

33.9 

14.22 

14.31 

(‘■'>12 

1.703 X 10 ^' 

4.418 X 10 <’ 

1.686 X 10 <' 

1.727 X 10 ^ 

(a)This table i 

s identical to Table V- 

14 . except that the 

' constants in this table 

are given in 


MHz. 

ib)Notc that for the AS-llF form the quartic constant 5j = 0. 
(c)Units for are in (MHz)-, all others are in MHz. 
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Table V-16 


GF/ENF Factors Best Ground State Fits 


AS Form 
Constant 

AS, 1*^ 

AS. IlF 

NS. 1^ 

NS. HI'' 

NS Form 
Constant 

XAS 

1.0 

1.0 

1.0 

1.0 

XNS 

yAS 

0.38 

1.0 

0.38 

0.93 

yNS 

2AS 

0.32 

0.022 

0.32 

0.13 

ZNS 


0.36 

0.36 

0.35 

0.36 

Dj 

^JK 

0.49 

0.050 

0.49 

0.041 

^JK 


0.20 

0.059 

0.20 

0.059 

*^R 


0.060 

0.021 

0.060 

0,39 

5i 


0.23 

0.021 

0.67 

0.48 

h 

Hj 

0.057 

0.052 

0.058 

0.057 

Hi 

Hjk 

0.25 

0.010 

0.25 

0.0067 

HiR 

Hkj 

0.62 

0.0031 

0.062 

0.0032 

Hkj 

Hr 

0.070 

0.0016 

0.070 

0.0016 

Hk 


0.011 

0.30 

0.011 

0.22 


^JK 

0.013 

0.013 

0.041 

0.13 

h2 


0.016 

0.0030 

0.015 

0.31 


Lj 


0.0040 

0.0026 


L) 

•-JJK 


0.0017 

0.055 


‘-JJR 

'-JR 

0,070 


0.015 


^R 

Lrrj 

0.011 


0.0085 


^kRJ 

Lr 



0.025 


l'k 


0.0036 

0.099 




■^JR 

0.0077 




8. 

«RJ 

0.0080 

0.0015 




«R 

0.075 

0.00082 



^4 

GMDT 

1.1 X 10 '"^ 

6.6 X 10 "^^ 

1.1 X 10 -^'^ 

1.2 X lO--"*' 
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Table V-17 


GF/ENF Factors - Best 1^2 Band Fits 


AS Form 
Constant 

AS, I*" 

AS, 111*^ 

NS, I*" 

NS, III^ 

NS Form 
Constant 

‘^0 

1.0 

1.0 

1.0 

1.0 

*^0 

XAS 

0.66 

0.62 

0.66 

0.63 

XNS 

yAS 

0.38 

0.38 

0.38 

0.36 

yNS 

^AS 

0.34 

0.37 

0.34 

0,37 

ZNS 


0.064 

0.064 

0.064 

0.064 

Dj 


0.13 

0.027 

0.13 

0.027 

Djk 

Ak 

0.24 

0.034 

0.24 

0.034 

Dr 


0.031 


0.031 

0.24 


% 

0.14 

0.096 

0.10 

0.66 

6. 

Hj 

0.027 

0.027 

0.027 

0.027 

Hi 

Hjk 

0.031 


0.030 

0.0072 

HiR 

Hrj 

0.032 

0.0076 

0.032 

0.0034 

Hkl 

Hr 

O.OoQ 

0.0065 

0.069 

0.0066 

Hk 


0.0088 


0.0088 

0.063 

^1 

‘’JK 

0.027 

0.053 

0.019 

0.16 



0.12 

0.042 

0.23 

0.52 

•'3 

CMDT 

1.5 X 10 -^-'’ 

2.1 X \ 0 ~^^ 

1.3 X 10 ’ 

3.2 X lO --^'^ 





orders of magnitude smaller for the AS-IIF form ground state than for the other Hamil- 
tonian forms.* The upper state GF/ENF factors are also smaller for the AS-llF form, 
though its Gram Determinant is not smaller since the AS-llF form fit contains three less 
constants thn*! ffee o^resrs. ftee 4ap©n4encies amo^rag the AS-HF form Hamil- 
tonian terms are consistent with it l^ng a m©re ste®% tei^e 

individual coefficients. It is well known, for example, that a power series approximation 
u.sing orthogonal p«.>lynomials, will have smaller aiefficients and will converge more rapidly 
than one using non-orthogonal polynomials. 

Tlie convergence properties of the AS-IIF form can also be studied by examining the 
unitary' transformation that leads to this form. The possible poor convergence of some re- 
duced Hamiltonians has been previously recognized only for nearly symmetric asymmetrre 
rotors. For the almost symmetric rotors C F, (k = 0.97d8> and SO, (k = -0.d42). Carpente4-‘'‘^^ 
found that the AS form Hamiltonian converged very slowly in the IIF representation for CFi 
and in the F representation for SOi. He t lien calculated sj | j. given in equation 11-22, for these 
different representations and found that the slow convergence of the Hamiltonian was related 
to a relatively large value of sm . This is expected for a nearly symmetric prolate or oblate 
top in the F and IIF rep-esentations. respectively, because sj j ] in the .AS form is proportional 
to l /lY-\), which becomes very large as Y -* X in these particular symmetric lop limits. 

Since the transformed Hamiltonian contains terms like S] j i (\-Z) in its quartic coefficients 
and like sjj | (Z-Y) and Sj 1 1 ( - 1\.,) in its sextic coefficients wc see that, in general, a 
bigger sj j [ will transform a larger part of the original lower order coefficients into the new 
higher order eoefficients. This will result in larger transformed coefficients and make the 
Hamiltonian power series converge more slowly. One would expect this type of behavior, 
though, to be noticeable only in nearly .symmetric asymmetric top molecules. 


The ground state tiram ITctcrininants compared here onlv included tH /l NF faelors for eoiustants up to 
sixtii order so they would be directly comparable to each other 



A similar calculation of sm was performed here using equations 11-22 and 11-25 for 
the AS and NS forms respectively, to see if it is much larger for the AS-IIF form than for the 
other three Hamiltonian forms. Unfortunately, sjj j cannot be calculated directly from the 
coefficients in a reduced Hamiltonian unless the planarity conditions for the quartic constants 
are impasad. the quartic pter^rity relations are not strictly correct, the errors in- 


troduced by using them in the calculation of sj ji were found to be small compared to the 
variation of sj j j with Hamiltonian form. All the untransformed quartic coefficients, and 
therefore Sj j j , can also be approximately calculated from the harmonic molecular force field 
if it is known. This calculation gives equilibrium distortion constants and not the vibrationally 
averaged ones observed in practice. For H 2 S, a molecular force field obtained solely from 
vibrational band centers is available.^^) Therefore, the S| 1 j values calculated with this 
method will contain no input from H^S rotational frequencies or from the fitting routines 
used in this work and can serve as a check on the other calculations of Sj u . 

To calculate sj j j from the coefficients obtained in the least-squares fits, the two planar- 
ity relations _ 


T 


ac 


1 


A^C^ 








(V-17) 


and 


Tbc = 


2r2 


B-'C 


A4 34 (-4 


(V-18) 


were used along with the equations for the determinable coefficients Tj and Tj given in equa- 
tion 11-18. The planarity relations given above are written in terms of the T^^^ constants since 
W = ’’Lice " '•'ac W = ’’bbcc " '•'be- Note that the third planarity relation for r^^b 

given in Appendix A cannot be written using the T^^ notation because - 2 * 

Tjj,. With the values given for the determinable constants Tjj, Tj,b- and T(.(. in Tables V-13 
and V-14, T^j. and T^^ can be calculated using the planarity relations. If only Tg^ or Tjj^. is 
calculated using these relations, the other constant as well as Tgj, can be obtained by solving 


the relations for T| and T 2 simultaneously. The procedure used here was to calculate the 


required coefficients twice, each time using a -tifferent planarity relation along with the 
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equations for Tj and T 2 . There are a number of additional ways to calculate sjjj using these 
equations, but they should all give roughly the same results, as was tested by using our two 
different methods for calculating Sjj j . The determinable constants A^, B^, and were 
used in equations V-1 7 and V-1 8 as (very good) approximations for A, B, and C. The results 
of these calculations are given in Table V-1 8 where sj ^ - 1 denotes the Sj j j value calculated 
using equation V-1 7 and the relations for Tj and Tt, and sj j j - 2 denotes the sj j j value cal- 
culated using equation V-1 8 and the Tj and T 2 relations. The best fit values for the quartic 
distortion constants were used in all cases. 

To calculate sj | j from the harmonic force field the following relations for the quartic 
distortion constants derived by Kivelson^^'^ for a planar triatomic molecule were 
used. 


^aaaa ^ 


A> 


= - + tan^0 tan 6 Fj 


-1 


(V-19) 


’’bbbb ^ 


B2 


= 2 F"j + cot- 0 F;\ + 2 v^2 cot 0 F‘^ 


(V-20) 


^aabb 

AB 


= 2 F"j - F"^ + (cot 0 - tan 0) F'* 


(V-21) 


where 


T.U.K R 2 M 


* abab 

”aF 


lUv 


(1+2 sin' 0)-' F:1 
lUx nix 


R = 


r- X 10-^0 
2h 


(V-22) 


0 is '/z the bond angle for H 2 S = 46.07°, M = atomic weight of H 2 S = 34 amu. lUy = rnhydfoge^ 

1 amu, m^ = = 32 amu, and r = H - S bond length = 1.3518 A. The inverse force 

field matrix F”' was calculated from a force field matrix given by Cook et. al.1^0) which 
was derived solely from the vibrational band centers of H^S. The non-zero elements of these 


matrices are 


Table V-18 


Sj j j Parameter Calculations 


Type of 
Data 

Hamiltonian 

Form 

Mll-1 

dll-2 

(X 10-^) 

s j j j_pp (c)Standard Deviation 

of Fit 

(MHz) 

MW 

ASd") 

62 

70 

56 

1.96 

MW 

AS(lir) 

837 

849 

729 

136 

MW 

NSd'") 

90 

82 

82 

2.91 

MW 

NS(111'’) 

-24 

-17 

-17 

2.09 

(a)MW 

AS(lllS) 

864 

861 

729 


(b)MW 

AS(11F) 

770 

812 

729 


IR 

AS(l^) 

69 

91 

53 

54 

IR 

AS (111^) 

1554 

962 

661 

5100 

IR 

NS(I'’) 

93 

75 

75 

54 

IR 

NS (111^) 

-37 

-17 

-16 

54 


^^^From MW fit of Helminger et al/ 

^^^From MW fit ofGiilisct al/*^^ 

^‘-'Standard deviation obtained when Hamiltonians include terms only up to sixth power in 
angular momentum. 


(V-23) 


and 


F,, =4.2731 md/A 
Fj 2 = 0.4250 md/A 
F33 = 4.2960 md/A 
Fi 2 = F2 j = 0.064 md/A 

F‘[ = 0.23455 (md/A)“‘ 

F*^ = 2.3582 (md/A)-’ 

F-’ = 0.23278 (md/A)-’. 

^12 " ^ -0 03532 (md/A)-’ 


(V-24) 


The rest of the required distortion constants were calculated using the planarity relations and 
the equations given in 11-8 relating the constants to the = 4 constants. The 
values for S| j j calculated u.sing the molecular force field are denoted by S| j j - FF in Table V-1 8. 

Another general check on the accuracy of the quartic planarity relations when applied 
to H2S was made by evaluating the planarity relations derived by Watson. These relations 
are stated directly in terms of the AS and NS Hamiltonian constants, and are given in Appendix 
A. equati m A-4. They involve four terms, composed of products of quartic and quadratic 
rotational constants, whose sum should equal zero if true equilibrium constants are used. 

These relations were evaluated for each of the four Hamiltonian forms in both the ground and 
upper vibrational states. A measure of the accuracy of these relations is the deviation from 
zero of ;he sum of the four terms divided by their average magnitude. Using this measure, the 
planarity relations were found to be accurate to within \ 7 < to 57< for the gic and state constants 
and 57 to 407 for the upper state constants. These results indicate that any errors in the cal- 
culation of sj 1 1 due to the planarity relations should be small compared to the variation of 
S| j| with Hamiltonian form observed in Table V-1 8. 

From Table V-18, we .see that sj ] | is about 1 2 times larger for the AS-llF form in the 
ground state and about 25 times larger for this form in the upper vibrational state than for 
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the other three Hamiltonian forms. The rather good agreement among the three calculated 
values for S|n for each form gives us confidence in the methods used for these calculations. 
Also, since the values for sjj| calculated from the distortion constants are always close to 
those calculated from the molecular force field, it would have been possible to predict the 
poor performance of the AS-lIH Hamiltonian in fitting the transition frequencies without 
performing the least-squares fits. The similarity of these two calculations is also anothe; 
strong assurance that the AS-lIU form fitting routine does not contain any program errors 
or numerical problems. It should also be noted that the values calculated for S| | j from the 
constants of Helminger et. al. and Gillis and Edwards are very close to our value for Sj j ] for 
the AS-lIF form. Considering the fact that molecular rotational data is usually only fit to a 
single, somewhat arbitrarily chosen reduced Hamiltonian, our results show the usefulness of 
performing a calculation of sj u before choosing a Hamiltonian form for the least-squares 
fit of the. data, sj | j - FF was also calculated for the AS and NS forms in the IF representa- 
tion and was found to equal - 55 X 10“^ for the AS form and - 31 X 10”^ for the NS form. 
Both of these values for S| jj are close to those obtained for the AS-F, NS-F, and NS-IIF 
Hamiltonians, and therefore fits made using the IF representation will probably be no better 
than those already performed. 

To illustrate how much the unitary transformation for the AS-IIF form can alter the 
rotational constants, we can evaluate the quartic coefficient Tq 22 , which corresponds to -6)^, 
and part of the sextic coefficient corresponds to Hj^. From equation 11-21 and 

the ground state AS-llF form constants we have 

-5^ = 257 MHz = To22=To22-2Bo2o ^u \ = Tq22 238 MHz (V-25) 
where the AS-llF form values for Bq 20 ^nd sj n have been inserted. Using the formula for 
0Q6O Table 2 of Ref. (3) and the value found for Hj^ in the AS-llF ground state fit we have 

Hk = 0.42 MHz = 00^0 = <> 060 ^^^ Bq 20 *11 1 

= 0QgQ - 1.3 MHz + other terms 
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(V-26) 



where the AS-lII* form values for Bq 20 and Sj j j have again been inserted. We see that the 
AS-lIU unitary transformation contributes 238 MHz to which can be compared to a con- 
tribution of only - 13.2 MHz to calculated for the AS-F form unitary transformation. 
Similarly, the AS-IIF transformation contributes - 1.3 MHz to Hj^ from the Bq 20 | 
term while the AS-F transformation only contributes 0.0058 MHz to from this term. 

We can also calculate Tq 22 froni equation V-25 to obtain Tq 22 = 19 MHz. Therefore the 
AS-llF unitary transformation changed Tq 22 by more than an order of magnitude, illustrat- 
ing the poor characteristics of this particular Hamiltonian form. The large contributions to the 
higher order transformed constants from the original lower order ones for the AS-llH form, 
combitied with the already large centrifugal distortion corrections required to accurately cal- 
culate the H 2 S rotational energy levels, results in a power series for this reduced Hamiltonian 
with extrem^'ly slow convergence. 

It is possible to understand the poor behavior of the AS-llF Hamiltonian in a more 
qualitative manner. Because of the geometry of H 2 S, the rotational constants A and B are 
sensitive to the value of the bond angle, while C is not. As the bond angle goes to 180°. 
for example, the A rotational constant diverges. Also, the force constant F 22 . which is 
associated with changes in the bond angle, is very small. Consequently, rotation of the 
molecule about a and b axes can significantly alter the bond angle and cause large changes 
in the instantaneous A and B rotational constants. Rotation about the c axis, though, 
only causes a small amount of bond stretching, which affects A, B, and C very little. 
Therefore, the centrifugal distortion constants that go with terms in the Hamiltonian con- 
taining the operators and will be large relative to those constants that go with terms 
containing only J^.. For example, from Table V-13, it can be seen that the quartic distor- 
tion coefficients T^^ and T^b about 25 times larger than T^.^.. 

If the operators present in the four reduced Hamiltonians used here are ex- 
amined (see equations 11-24 and 11-27). a striking difference between the AS-llH Hamil- 
tonian and the other three Hamiltonian forms is apparent. To any order, the highest 
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powers of and J 5 that ever appear in the AS-IIK Hamiltonian are and The other 
three Hamiltonian forms, though, when expanded to the n'^ power of angular momentum, 
contain at least one of the operators or or the product Since the opera- 

tors with the largest distortion coefficients have been dropped from the AS-lIH Hamilton- 
ian, it is not surprising that it fits the transition frequencies so poorly compared to the 
other three Hamiltonian forms. As has been shown, the AS-llH Hamiltonian must be 
taken to higher orders of angular momentum than the other three Hamiltonians to make 
up for its lack of terms with high powers of Jg and J 5 . 

In conclusion, a very satisfactory fit of the measured V 2 band frequencies has been 
made to the AS-1*^ form reduced Hamiltonian for the upper vibrational state. It should 
be mentioned that tic* use of the AS-1'' constants to predict transition frequencies involv- 
ing states with J or Kg greater than about 9 (roughly the maximum value of J or Kg used 
in the least-squares fits) will result in errors higher than those indicated by the statistical 
uncertainties in these constants. The AS-lIl*’ reduced Hamiltonian fit the measured fre- 
quencies very poorly due to its very slow convergence. I'his unexpected slow conv^^gence 
could have been predicted by the large value of the Hamiltonian transformation parameter, 
siji , calculated using the molecular force field constants. In most previous studies of 
asymmetric rotor spectra, a single and somewhat arbitrary choice for the Hamiltonian 
form has usually been made. <''>ur re.->ults. though, illustrate that in some cases, it can be 
very important to choose the proper Hamiltonian form and representation, even when 
fitting rotational data for very asymmetric tops. 
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CHAPTER VI 


LINE STRENGTHS IN THE 1^2 BAND OF H 2 S - COMPARISON 
TO RIGID ROTOR LINE STRENGTHS 

In this chapter the results of the H 2 S line strength measurements will be discussed 
and compared to calculated rigid rotor line strengths. Since a number of doublets in the 
t '2 band of H 2 S havo lines with an intensity of 3 to 1 due to nuclear spin statistics, a par- 
tial check of the internal consistency of the measured line strengths will be made by com- 
paring the observed line strength ratios to the known ratio of 3 to 1. Some general obser- 
vations will also be made on the dependence of the line strength perturbations on J. Ka 
and K(;. 

The line strength ratios in twelve doublets, all measured at low pressures, were com- 
puted for comparison to the expected ratios of 3 to I . The expected deviation of any of these 
strength ratios from 3 to 1 due to centrifugal distortion or other causes should be small 
since the doublets were split by only 0.01 enr* to 0.1 cm*’. The standard deviation of the 
observed ratios from a 3 to 1 ratio was found to be l.H7r. Since both lines of each doub- 
let were always in the same diode laser scan, any errors in their strengths due to uncertain- 
ties in the gas pressure or due to distortion of the lines by the laser lineshape should be 
the same. Therefore, the experimental errors in the observed line strength ratios were 
determined solely from the standard deviation of the measurements used to obtain an 
average line .strength for each transition. The average experimental error calculated in this 
manner for the twelve ratios was 1.89?. the same as the standard deviation oi the observed 
ratios from a 3 to 1 ratio. This indicates that the strengths measured within one laser 
mode were internally consistent. 

The line strengths measured here were compared to calculated rigid rotor line strengths 
in order to determine the magnitude of the centrifugal distortion perturbations. To calcu- 
late the rigid rotor line strengths the total band strength S^. which is the sum of alt the 


127 



line strengths in the band, is needed. Since Sy cannot be determined directly from our 
measured line strengths, the band strength value of S° = 2.17cm'Vatm-cm measured by 
Emerson and Egger^*) was used for the rigid rotor calculation. In terms of the total 
band strength, the rigid rotor line strengths are given b> 


S = - (1 -e-h^/kT) 

»^Qr 

This equation can be derived from equation 11-32 and the definition S°. V is the band 
center frequency. In this approximation, S° is proportional to the vibrational part of the 
dipole moment matrix element, which is given in equation 11-40. 

The usual approximation was made in equation VI- 1 that the partition function Q 
can be separated into a rotational part, Qr and a vibrational part, Qy, which is included 
in Sy. Qr is given approximately by(^3) 


Or == 



(VI-2) 


wb‘,e A, B, C are the ground state quadratic rotational constants. The accuracy of equa- 
tion VI-2 was checked by comparing the value it gives for Qr for H 2 O to the more accurate 
value of Qr calculated by Camy-Peret and Flaud'^^) using equation 11-33. These two 
values of Qr for H 2 O were different by only 2.3%. This uncertainty of about 2% in Or 
and the quoted uncertainty of 5% in the total band strength were the main contributors to 
any errors in the rigid rotor line strength calculation. 

The lower state energy levels, Ej, and the direction cosine matrix elements ^ 
were calculated using a computer program written by Maki.(^) This program uses a 
Watson AS-F Hamiltonian to calculate the rotational energy levels and eigenfunctions required 
to determine these matrix elements. Since only terms up to the sixth power of angular 
momentum were included in the ground and upper state AS-F Hamiltonians for these cal- 
culations. our best fit rotational constants for the same ground and upper state Hamilton- 
ians were used here. It should be noted that by using quartic and sextic terms in the rotational 
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Hamiltonian in this calculation, the average centrifugal distortion within each vibrational 
state is included in the rotational part of the total wavefunction. The vibrational part of the 
dipole moment matrix element, though, is kept constant for all the indhidual line strengths 
calculated for the band. The final rigid rotor line strengths calculated using equation VI-1 
were divided by 0.9502 to correct for the naturally occurring abundancef^^^ of the sulphur 
32 isotope. Since Hmerson and Eggers made no mention of isotopic impurities in their 
measurements, we have assumed that their quoted value for includes contributions 
from all the sulphur isotopes. 

The calculated rigid rotor line strengths and the observed line strengths are plotted as 
a function of frequency in Figure Vl-1. The experimental F factors (the ratios of the ob- 
served to calculated line strengths) for each line are listed in Table IV-5. Also, the F 
factors are plotted as a function of frequency in Figure Vl-2. The 0 branch plots have 
been separated from the R and P branch plots in this figure for clarity. It should be 
remembered that only the H 2 S transitions measured with the diode laser spectrometer are 
shown in these figures. There are a number of other H^S absorption lines in the 1080cnr* 
to 1260cnr* region th.;t were not measured and do not appear in Figures VI- 1 and VI-2. 

As expe' ted fa>m the low resolution measurements of this band, the R and 
branch F factors are greater than unity, and the P and branch F factors are almost all 
lesj than unity. The F factors range from as low as 0.057 on the P branch side of the band 
center to as high as 4.71 on the R branch side. The maximum values of J, K;,. and K^. 
observed for the measured transitions were 12. 7. and 10. respectively. It is also evident 
in Figure Vl-1 that in all the branches it is the weaker lines whose stpuigths differ the 
most relative to their rigid rotor strengths. 

The strong lines in Figure Vl-1 . whose strengths differ from the rigid rotor line 
strengths by only 1091 to 209'. belong to the four sub-bands I J)qj -(J-Di j-i.lJ-D] j.i - 
(J)oj . - U-1 ^o.J-l ' rotational energy levels involved in 

these transitions always have one slate with Kg = I and the other with = 0. The values 
are always equal to J lor each energy level, the maximum jxjssible value allowed. This 
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Figure VI- 1 . Stick Spectrum of Observed Line Strengths and Rigid Rotor Line Strengths 
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Figure VI-2 Ratio ol Observed Line Strengths to C alculated Rigid Rotor Line Strengths CF Factor) 




corresponds classically to rotation almost exclusively about the c axis of the molecule in both 
the upper and lower states of the transition. To a large degree, rotations about the c axis 
stretch the H-S bonds but do not change the bond angle. Centrifugal forces arising from ro- 
tations about the a and b axes, though, will alter the bond angle. Since the H2S force con- 
stantst^®^ involving bond an^le changes are so much smaller than those involving bond stretch- 
ing, rotations about the a and b axes distort the molecule more severely and therefore affect 
the dipole moment the most. Thus it is not surprising that the transitions with K^. at a maxi- 
mum in both the upper and lower rotational states have line strengths close to rigid rotor 
values. 

It is also expected that lines with a large change in Kjj and would have F factors 
very different from unity. Unfortunately lines with and AKj. greater then one are al- 
ways very weak and hard to observe. Our only example of a line with both AKg and 
AK^. greater than one, the 54]-5j4 transition, had a F factor of 4.71. the largest one 
observed. For this transition, the rotational angular momentum about the a axis is much 
greater in the upper state than in the lower state. Therefore, the bond angle is smaller in the 
upper state. This difference in the molecular geometry between the upper and lower states 
results in a much larger change in the dipole moment for this transition relative to a rigid rotor 
transition. 

The largest and smallest F factors were usually for the and 0***^ lines respectively. 
This is partially because the observed Q branch lines had higher values of J and than the 
observed R and P branch transitions. Since the distortion of the molecule is proportional 
to the square of the components of the total angular momentum, the difference in the 
distortion between two states with high J and Ky values is greater than the difference in the 
distortion between states with low values of J and Ky. In general, this leads to more 
strongly perturbed line strengths at high J and Ky values. 

As mentioned in (’h. 11. it is hoped that the line strength measurements presented here 
can serve as a severe test of the presently existing theories of centrifugal distortion in asym- 
metric rotor line strengths. The combination of the large strength perturbations present in 



the vj accuracy of these measurements should make these line 

strengths very useful for that purpose. 
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CHAFTKR VII 


NUCLKAR QUADRUPOLF HYPERFINK STRUCTURE IN THE INFRARED 
SPECTRUM OF HYDROC.EN IODIDE 


VII. A. Introduction 

NuJear quadrupolc hyporfinc structure in the 1-0 vibration-rotation band of hydrogen 
iodide near 4.5 iim was observed using a tunable semiconductor diode laser.* Hypertlne 
structure arising from nuclear quadru^xile interactions is commonly seen in micarwave aita- 
tion sjx'ctra with splittings ranging from less than 1 MHz to several hundred MHz for mole- 
cules amtaining iodine. In vibration-aitation transitions, nuclear quadrupole splittings have 
been observed only in the past several years using the sub-Doppler resolution techn ques of 
laser saturation spectarsaipy in studies of methyl halide spectra.* In the case of hydrogen 
iodide however, its extremely large nuclear quadrupole moment and molecular weight cx>m- 
bine to give relatively large hyperfine splittings with small Doppler widths. This allows reso- 
lution of some In perfine structure in the low J transitions of hydrogen iodide with Doppler 
limited spectrosarpy. Several studies of the 1-0 band of the HI molecule have been made 
with grating spectrometers with resolution too low to observe any hyperfine structure. The 
grating sjvctra were used in locating the transitions studied in the present work with the 
diode laser.* 

VII. B. Nuclear Quadrupsde Interaction Theory 

Nuclear quadrupole splittings arise through the interaction of the nuclear electric quad- 
rvqxile moment with the electric field gradient of the molecular electron cloud. This inter- 
action couples the nuclear spin. I (if 1 > I) to the rotational angular momentum J to form 
the total angular momentum F = 1 + J and splits each rotational level into 21 + 1 levels for 
J ^ 1 or 2J + 1 levels for 1 > J. These new levels each have different values of total angular 
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momentum denoted by the quantum number F, which varies from I J - 1 I to I J + 1 I in 
steps of one. 

The interaction energy is given by^^^ 

= (Vll-l) 

where 

Qij = f p(r)l3rir3-£-6ijl dr, (VIl-2) 

•'nucleus 

is the nuclear quadrupole moment tensor and 

a-v aU; 

V = 1 (Vll-3) 

^ arjdrj 3rj 

is the molecular electric field gradient at the nucleus. p(r> represents the nuclear charge den- 
sity and the rj are space-fixed cartesian coordinates. The dependence of the interaction 
energy on the rotational state of the molecule was originally derived by Casimit^^^^ using 
first order perturbation theory, and can be found in several texts.^-'*-^'^^ Casimir obtained 


<IJFIHq 11JF> = - eQtij 


3/4C(C+ n-UI U JU 


U21- n JCJ - 1) 


+ \) 


(Vll^) 


where 


C F(F+ n- 1(1 + 1)-J(J + 1). (Vll-5) 

eQ is called the “quadrupole moment of the nucleus” and is defined in the limiting case of 
M| = I as 

eQ = <1,M| = 1 !plr)(3r^ -i") 11, M, = l>. tVll-b) 

Similarly, qj is defined as the value of the molecular electric field gradient in the z direction 
averaged over the state 1 J. Mj = J >. or 

qj = < J, Mj = J IV„1J. Mj = J >. (Vll-7) 

For a diatomic molecule 


- J 

" CJ + 3) 


(Vll-8) 
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where q is the field gradient along the molecular bond axis. The rotational level splittings 
for a diatomic molecule can then be written as 

AEq = <IJFIHq 11JF> = - (eQq) Y(I, J, F) (VII-9) 

where 


3/4C(C-H)-l(I+ 1) J(J-H) 
21(21- 1)(2J - 1)(2J + 3) 


(Vll-10) 


The combined constant eQq is called the nuclear quadrupole coupling constant. The ob- 
served splittings in a vibration-rotation transition will be given by 

A;; = (eQq)” Y(lJ'',F")-(eQq)' Y(U',F') (Vll-11) 

where a double prime indicates the upper state and a single prime the lower state of the trans- 
iiion. Though Y(l, J, F) in general varies for each hyperfine component and rotational level, 
eQq will only vary with vibrational state. 

The relative intensities of the various hyperfine components within a given vibration- 
rotation transition are as follows/ ^0) 


^Q(F)Q(F- 1) AF = ^ 1 
F 

P(F)Q(F) AF = 0 (Vll-12) 

F(F + 1) 

I P(F)P(F- 1) AF = ± 1 
F 

where 


P(F) = (F + J)(F + J+ 1)-1(1+ 1) 
Q(F) = 1(1+ 1)-(F-J)(F-J + 1). 


(Vll-13) 


The top sign on AF is to be used for AJ = + 1 (R Branch) transitions and the bottom sign for 
AJ = - 1 (P Branch) transitions. 

The '^^1 nuclear spin of 5/2 and microwave determined ground vibrational state eQq 
of - 1828.4 for hydrogen iodide give line splittings in the pure rotation spectrum as 

targe as 800 MHz. Splittings in the hydrogen iodide vibration-rotation spectrum will be very 
similar to the rotational splittings, modified slightly due to any dependence of eQq on the 
vibrational slate of the molecule. The selection rules of At = 0, AF = 0, ± 1 gives 3 lines for 
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R(0) and P(l), 9 lines for P(2), and 14 lines for P(3V Stick spectra for P(l), P(2), and P(3) 
calculated using the known ground state eQq and the first excited vibrational state eQq value 
measured in this work are given in Figure VII-1 . Only 1 2 lines are shown for P(3) because 
two pairs of lines are degenerate to first Oider, although they are split by about 1 .5 MHz, due 
to nuclear magnetic hyperfine effects,^^* ) which is too small to be shown in Figure VIl-1 (c). 

VII. C. Experimental Considerations 

The R(0) through P(3) lines, which fall in the region from 2190 to 2242 cm"*, were re- 
corded using a PbSSe diode laser. The experimental apparatus was almost identical to that 
described in Ch. HI except for the use of a one meter single pass cell instead of the White 
cell. Also the germanium etalon was not thermally stabilized during this work, which should 
not adversely affect the splitting measurements since the splittings are so small. 

The gas cell pressures used ranged from 3 to 7 torr. The small strength of the HI transi- 
tion^^2) and the somewhat impure gas samples used required these relatively high pressures in 
the 1 -meter cell. This resulted in linewidths generally 307r greater than the Doppler width 
of HI. It is unlikely that pressure shifts could introduce any systematic errors since all hyper- 
fine levels in a given transition should have very similar pressure-shift coefficients. 

The line centers, or minima within the hyperfine structure of a line, and etalon peaks 
were found digitally by least squares fitting the data to a cubic polynomial as described in 
Ch. IV. The splittings were determined using the calculated free spectral range for the etalon. 
Errors in the calculated free spectral range are negligible in this experiment since the split- 
tings measured were so small. Random laser frequency fluctuations between fringe peaks was 
the main source of error, giving standard deviations of 3 to 10 MHz in the measured splittings. 

Vll. D. Results 

Figure VII-1 (a) gives the observed and calculated profiles of the P(l) transition. The 
triplet structure due to splitting in the J = 1 rotational level is fully resolved. Since only the 
ground vibrational state is split in this transition, the infrared data can be directly compared 
to microwave measurements. The microwave value for eQq was used to generate the 
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Figure VII-1 : Observed and Calculated Profiles of P( 1 ), P{2), 
and P(3) of HI. The solid line is the observed profile and the 
dashed line the calculated profile. The splittings measured are 
labeled A and some individual transitions are identified by 
the quantum number F. Above each profile is a stick spec- 
trum of the transition with a line length proportional to the 
relative strength. 
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calculated profile in Figure VlI-1 (a). A Voigt line shape (sec Ch. IV, equation IV-6) was 
used, arbitrarily scaling the line strength to closely approximate the observed absorption. 

The overall agreement, both for splittings and relative strengths, is very good. The splittings 
were corrected for slight blending by determining the amount of line shift due to blending in 
the synthetic Voigt profile. The calculated splittings also included small corrections for mag- 
netic hyperfine contributions, averaging about l.SMHz, which is somewhat below the stand- 
ard deviations of the measurements. Second order nuclear quadrupole hyperfine contribu- 
tions were negligible Seven runs of the P( 1 ) line were made and their scatter about the 
mean is the main contributor to the stated errors. Table VlI-1 gives results for the measured 
splittings. Both splittings were measured relative to the strongest hyperfine component, 

F = 7/2 -* 5/2. They are compared in Table VII-l to those calculated using the well known 
microwave determined ground state eQq. Agreement is within experimental error for both 
splittings. 

The R(0), P(2), and P(3) transitions all involve rotational level splittings in the first ex- 
cited vibrational state, allowing a determination of eQq in that state. As J increases, the 
intensity of AF = AJ transitions become dominant over AF = 0 and AF = - AJ transitions, 
and the spectra condense to the unperturbed line center. This causes an extreme amount of 
line blending, especially for P(2) and P(3), so that quadmpole coupling constants could not 
be directly obtained from the observed absorption features. In order to convert measured 
frequency intervals between minima in the spectrum of each line to an eQq value for the 
upper vibrational state, simulated Voigt profiles were made for each lOtational line, keeping 
the ground state eQq constant while slightly varying the upper state cQq about the ground 
state value. From these profiles the change in observed splitting per unit change in upper 
state eQq was generated. Splittings analogous to those given by Aj and A2 in Figure VII-1 
(a) for P(l) were measured in R(0) to obtain a value for the upper state cQq, while the 
splittings identified by A’s in Figure VII- 1 (b) and VII-1 (c) were used for the P(2) and P(3) 
measurements. These spectral features were a combination of lines exhibiting the least 
amount of blending and greatest sensitivity to upper state eQq value. The theoretical spectra 
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Tabic VII-1 


Observed Splittinp in P(l) of HI 



Calculated - 

Splitting (MHz)* 

Observed (MHz) 

A, » 159.9 ±3.2 

2.6 

A2 = 3S5.5 ± 6.5 

-0.3 


•7 Tritl Average 


shown in Figure VIl-l (b) and VIl-1 (c) were generated in the same manner as for P(l), using 
both the microwave derived ground state eQq value and the measured eQq in the first excited 
vibrational state. As with P( I ), the overall agreement is very good for these transitions. Ivote, 
for example, the close agreement between the calculated and measured curves for the widely 
split but weak lines in P(3). Tlie R(0) profile, which is not shown, appears identical in form 
to P(l) except reversed in frequency. 

Table VIl-2 gives the measured differences in eQq between the ground and fu st excited 
vibrational state using the method described above along with the first excited state cQq values 
themselves. The average value for the difference was found to be very small, as expected, not 
much larger than its standard deviation. Estimates for the standard deviations of the two 
weighted averages in Table VII-2 were derived from the standard deviations of the four indi* 
vidual measurements and not from the spread of these values about the weighted average. 

This was done because the uncertainties in the individual eQq values were larger than their 
spread about the mean. The last column in Table VII-2 illustrates the sensitivity of the meas- 
ured splittings to the upper state eQq value, showing that P(2), for example, has the most 
sensitive splitting measured. In order to reveal any systematic errors in the splittings the 
measured difference in the ground state eQq from the much more accurate microwave value 
is also given in Table VII-2. This observed difference was determined from the P( I ) transi- 
tion. Its value, which lies less tlian one standard deviation away from 0%, indicates that there 
are no large systematic errors contributing to the upper state eQq determination. The meas- 
ured difference between the ground and first exdted vibrational state eQq of 1.2 i 0.7% is 
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Table VII-2 


Observed Difference in Quadrupoie Coupling 
Constant from Ground to 1st Excited Vibrational State 


Transition 

No. of Trials 

Excited State eQq 
(MHz) 

eQq Difference 
(%) 

Change in 
Splitting for 
1% Difference 
in e(^ (MHz) 

R(0)A, 

5 

-1854 ± 27 

1.4 ± 1.5 

1.8 

A 2 


-1853 ±48 

1 .4 ± 2.6 

3.8 

P(2) 

4 

-1851 ± 16 

1.2 ± 0.9 

5.8 

P(3) 

3 

-1830 ±44 

0.1 ± 2.4 

1.6 


Weighted Averages: 

-1850 ± 13 

1.2 ±0.7 (-22 ± 13 MHz) 

P(l) 

7 


-0.6 ± 1.3* 



•Measured DiffeierKC in Ground State eQq from Microwave Value 


comparable in magnitude to the vibrational dependences found for eQq in methyl iodide^^^*^^) 
using laser saturated absorption experiments, where eQq was found to change for various vi- 
brational modes from -.3% to +1.5% per unit change in vibration quantum number. 

In summary, nuclear quadrupoie hyperfine structure has been observed in the vibration- 
rotation spectrum of hydrogen iodide. The measured splittings are in good agreement with 
microwave measurements of HI. The relative intensities observed for the hyperfine compo- 
nents were shown to be reproduced very well using the theoretical expressions for the relative 
strengths. Also, some evidence for a slight increase in the quadrupoie coupling constant from 
the ground to first excited vibrational state was found. 
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APPENDIX A 


PLANARITY RELATIONS AND DETERMINABLE CONSTANTS 


A number of approximate relations can be derived for the rotational constants of a 
planar asymmetric rotor. It can be easily shown that the equilibrium quadratic constants 
obey the relation 


AeBe 
Ae + Be 


(A-1) 


since 


1/ = V + I^'. (A-2) 

These -elations, which are only strictly valid for the equilibrium rotational constants (denoted 
by the superscript e), are approximately correct for the vibrationally averaged constants 
that are observed in practice. For example, equation A-1 is accurate to about 2% when the 
observed ground state rotational constants of H 2 S are used. 

Similar planarity relations can be derived for the equilibrium quartic constants by util- 
izing the properties of the inertial derivatives of a planar molecule. These relations, originally 
derived by Dowling(^2)^ are 


^ I ^**23 ’’bbbb ^cccc' 

r~~bb = ~ A^ I 7 — + — ~ 

2 V a4 B4 C4 j 


'aacc 


a2c2| 


’ 3233 


’’bbbb 


'cccci 


a4 B^ 


^bbec 


= - B2c2 (- 


^aaaa ^ ^bbbb 
A* B^ 


^ccc c 

j 


(A-3) 


^acac ^bebe 


Therefore for a planar molecule, the nine quartic coefficients are reduced to four independent 
coefficients. In the above relations the superscript denoting equilibrium coefficients has been 
dropped. 
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These quartic planarity relations :an also be applied to the quartic distortion constants 
in Watson’s AS and NS reduced Hamiltonians to obtain a single relation among the fine re- 
duced quartic constants. The relations for AS and NS Hamiltonians in the P and IIF repre- 
sentations are^^‘^ 

AS-F: 4CAj - (B - OAj- - 2(2A + B + C)6j + 2(B - C)6|c = 0 

AS-IIF: 4CAj + (A + B + 2 C)Ajk + 2(A + B)Ak + 2(A - B) {Sj + 6^) = 0 

(A-4) 

NS-F; 4CDj - (B - ODjk + 2(2A + B + C)6i - 4(4A + B - 306, = 0 
NS-IIF; 4CDj + (A + B + 2 C)Djk + 2(A + BID^ - 2(A - B)6i = 0. 

Again, these relations are only strictly true for equilibrium quadratic and quartic constants, 
but should be approximately correct when the observed rotational constants are used instead. 

A set of quadratic and quartic determinable constants in terms of the AS and NS re- 
duced form constants have been derived by Watson^^), and arc as follows. 

AS Form 

XD = xAS + 2 Aj + Ajk - - 28k 

yD = yAS + 2Aj + Ajk + 25j + 25 k 
zD = zAS + 2 Aj 


T'xx = -Aj - 26 j 
Tyy = - Aj + 25j 
Tzz = - Aj - Ajk - Ak 


(A-5) 


T| - -3Aj - Ajk 

T2 = - (X + Y + Z)Aj - \{\ + Y)Ajk + (X - Y) (6j + 6k) 
NS Form 

X» = XNS + 2 Dj + Djk -r 26, + 482 
yD ^ yNS + 2Dj + Djk - 26, + 450 
7P = ZNS + 2Dj + 662 
Txx = -Dj+ 26, + 26. 


!A-6) 


Tyy = -Dj - 26 , + 262 
^73. = - Dj - Djk - Dk 
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T, = - 3 Dj - Djk - 662 

T2 * -(X + Y + Z)Dj - -7(X + Y)Djk -(X-Y)6i -6Z62- 


(A-6) 

cont. 


The X, Y, Z constants in the above relations for Tj and T2 correspond to the X, Y, Z con- 
stants in equation 11 - 10 . In this work X^, Y*^, and Z^ were used in place of X. Y, and Z. 
This should not introduce any errors in these relations larger than any already present. 
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APPENDIX B 


CORRELATION TABLES FOR MOLECULAR CONSTANTS OF 
GROUND AND 1^2 = 1 EXCITED STATES OF H2S 
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Table B-1 

Correlation Table, AS Form, V Representation, pj = ’ State Fit, Ground State Fixed at Best Fit Value 
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